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The BaBAR Detector 

The BABAR Collaboration 



BABAR, the detector for the SLAC PEP-II asymmetric e"'"e~ B Factory operating at the T(4S') resonance, was 
designed to allow comprehensive studies of CP- violation in _B-meson decays. Charged particle tracks are measured 
in a multi-layer silicon vertex tracker surrounded by a cylindrical wire drift chamber. Electromagnetic showers 
from electrons and photons are detected in an array of Csl crystals located just inside the solenoidal coil of a 
superconducting magnet. Muons and neutral hadrons are identified by arrays of resistive plate chambers inserted 
into gaps in the steel flux return of the magnet. Charged hadrons are identified by dE/dx measurements in the 
tracking detectors and in a ring-imaging Cherenkov detector surrounding the drift chamber. The trigger, data 
acquisition and data-monitoring systems , VME- and network-based, are controlled by custom-designed online 
software. Details of the layout and performance of the detector components and their associated electronics and 
software are presented. 



1. Introduction 

The primary physics goal of the BABAR ex- 
periment is the systematic study of CP-violating 
asymmetries in the decay of neutral B mesons 
to CP eigenstates. Secondary goals are precision 
measurements of decays of bottom and charm 
mesons and of r leptons, and searches for rare 
processes that become accessible with the high 
luminosity of the PEP-II B Factory Q. The de- 
sign of the detector is optimized for CP violation 
studies, but it is also well suited for these other 
physics topics. The scientific goals of the BABMi 
experiment were first presented in the Letter of 
Intent |^ and the Technical Design Report [||; 
detailed physics studies have been documented 
in the BABMi Physics Book [|j and earlier work- 
shops [||. 

The PEP-II B Factory is an asymmetric e+e^ 
collider designed to operate at a luminosity of 
3 X 10'^^ cm~^s~^ and above, at a center-of- 
mass energy of 10.58 GeV, the mass of the T{AS) 
resonance. This resonance decays exclusively to 
B^B^ and B^B~ pairs and thus provides an ideal 
laboratory for the study of B mesons. In PEP- 
II, the electron beam of 9.0 GeV collides head- 
on with the positron beam of 3.1 GeV resulting 
in a Lorentz boost to the T(4S') resonance of 
/37 = 0.56. This boost makes it possible to recon- 



struct the decay vertices of the two B mesons, to 
determine their relative decay times, and thus to 
measure the time dependence of their decay rates. 
The crucial test of CP invariance is a comparison 
of the time-dependent decay rates for B° and B° 
to a self-conjugate state. For the cleanest exper- 
imental test, this requires events in which one B 
meson decays to a CP eigenstate that is fully re- 
constructed and the other B meson is tagged as 
a B" or a B" by its decay products: a charged 
lepton, a charged kaon, or other flavor sensitive 
features such as a low momentum charged pion 
from a D* decay. 

The very small branching ratios of B mesons to 
CP eigenstates, typically 10~^, the need for full 
reconstruction of final states with two or more 
charged particles and several 7r°s, plus the need 
to tag the second neutral B, place stringent re- 
quirements on the detector, which should have 

• a large and uniform acceptance down to 
small polar angles relative to the boost di- 
rection; 

• excellent reconstruction efficiency for 
charged particles down to 60 MeV/ c and for 
photons to 20 MeV; 

• very good momentum resolution to separate 
small signals from background; 
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Figure 1 . BABAR detector longitudinal section. 



• excellent energy and angular resolution for 
the detection of photons from 7r°and 77° de- 
cays, and from radiative decays in the range 
from 20 MeV to 4 GeV; 

• very good vertex resolution, both transverse 
and parallel to the beam direction; 

• efficient electron and niuon identification, 
with low misidentification probablities for 
hadrons. This feature is crucial for tagging 
the B flavor, for the reconstruction of char- 
monium states, and is also important for 
the study of decays involving leptons; 

• efficient and accurate identification of 
hadrons over a wide range of momenta for 



B flavor-tagging, and for the reconstruction 
of exclusive states; modes such as 
K^TT^ or B^ 7r+7r~, as well as in charm 
meson and r decays; 

• a flexible, redundant, and selective trigger 
system; 

• low-noise electronics and a reliable, high 
bandwidth data-acquisition and control sys- 
tem; 

• detailed monitoring and automated calibra- 
tion; 

• an online computing and network system 
that can control, process, and store the ex- 
pected high volume of data; and 
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Figure 2. BABAR detector end view. 



• detector components that can tolerate sig- 
nificant radiation doses and operate reliably 
under high-background conditions. 

To reach the desired sensitivity for the most in- 
teresting measurements, data sets of order 10^ B 
mesons will be needed. For the peak cross section 
at the T(4S') of about 1.1 nb, this will require an 
integrated luminosity of order 100 fb~^ or three 
years of reliable and highly efficient operation of 
a detector with state-of-the art capabilities. 

In the following, a brief overview of the princi- 
pal components of the detector, the trigger, the 
data-acquisition, and the online computing and 
control system is given. This overview is followed 
by brief descriptions of the PEP-II interaction re- 
gion, the beam characteristics, and of the impact 
of the beam generated background on the design 



and operation of the detector. Finally, a detailed 
presentation of the design, construction, and per- 
formance of all BABAR detector systems is pro- 
vided. 

2. Detector Overview 

The BABAR detector was designed and built by 
a large international team of scientists and en- 
gineers. Details of its original design are docu- 
mented in the Technical Design Report , issued 
in 1995. 

Figure ^ shows a longitudinal section through 
the detector center, and Figure || shows an end 
view with the principal dimensions. The detector 
surrounds the PEP-II interaction region. To max- 
imize the geometric acceptance for the boosted 
T{4:S) decays, the whole detector is offset rela- 
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tive to the beam-beam interaction point (IP) by 
0.37 m in the direction of the lower energy beam. 

The inner detector consists of a sihcon ver- 
tex tracker, a drift chamber, a ring-imaging 
Cherenkov detector, and a Csl calorimeter. These 
detector systems are surrounded by a supercon- 
ducting solenoid that is designed for a field of 
1.5 T. The steel flux return is instrumented for 
muon and neutral hadron detection. The polar 
angle coverage extends to 350 mrad in the forward 
direction and 400 mrad in the backward direction, 
defined relative to the high energy beam. As in- 
dicated in the two drawings, the right handed co- 
ordinate system is anchored on the main tracking 
system, the drift chamber, with the z-axis coin- 
ciding with its principal axis. This axis is offset 
relative to the beam axis by about 20 mrad in the 
horizontal plane. The positive y-axis points up- 
ward and the positive a;-axis points away from the 
center of the PEP-II storage rings. 

The detector is of compact design, its trans- 
verse dimension being constrained by the 3.5 m el- 
evation of the beam above the floor. The solenoid 
radius was chosen by balancing the physics re- 
quirements and performance of the drift chamber 
and calorimeter against the total detector cost. 
As in many similar systems, the calorimeter was 
the most expensive single system and thus consid- 
erable effort was made to minimize its total vol- 
ume without undue impact on the performance 
of either the tracking system or the calorimeter 
itself. The forward and backward acceptance of 
the tracking system are constrained by compo- 
nents of PEP-II, a pair of dipole magnets (Bl) 
followed by a pair of quadrupole magnets (Ql). 
The vertex detector and these magnets are placed 
inside a support tube (4.5 m long and 0.217 m in- 
ner diameter) that is cantilevered from beamline 
supports. The central section of this tube is fab- 
ricated from a carbon-fiber composite. 

Since the average momentum of charged par- 
ticles produced in B-meson decay is less than 
1 GeV/ c, the precision of the measured track 
parameters is heavily influenced by multiple 
Coulomb scattering. Similarly, the detection effi- 
ciency and energy resolution of low energy pho- 
tons are severely impacted by material in front 
of the calorimeter. Thus, special care has been 
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Figure 3. Amount of material (in units of radia- 
tion lengths) which a high energy particle, origi- 
nating from the center of the coordinate system 
at a polar angle 9, traverses before it reaches the 
first active element of a specific detector system. 

taken to keep material in the active volume of 
the detector to a minimum. Figure ^ shows the 
distribution of material in the various detector 
systems in units of radiation lengths. Each curve 
indicates the material that a high energy particle 
traverses before it reaches the first active element 
of a specific detector system. 

2.1. Detector Components 

An overview of the coverage, the segmentation, 
and performance of the BABAR detector systems 
is presented in Table |l|. 

The charged particle tracking system is made of 
two components, the silicon vertex tracker (SVT) 
and the drift chamber (DCH). 

The SVT has been designed to measure angles 
and positions of charged particles just outside the 
beam pipe. The SVT is composed of five layers 
of double-sided silicon strip detectors that are as- 
sembled from modules with readout at each end, 
thus reducing the inactive material in the accep- 
tance volume. The inner three layers primarily 
provide position and angle information for the 
measurement of the vertex position. They are 
mounted as close to the water-cooled beryllium 
beam pipe as practical, thus minimizing the im- 
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pact of multiple scattering in the beam pipe on 
the extrapolation to the vertex. The outer two 
layers are at much larger radii, providing the co- 
ordinate and angle measurements needed for link- 
ing SVT and DCH tracks. 

The principal purpose of the DCH is the mo- 
mentum measurement for charged particles. It 
also supplies information for the charged particle 
trigger and a measurement of dE/dx for particle 
identification. The DCH is of compact design, 
with 40 layers of small, approximately hexagonal 
cells. Longitudinal information is derived from 
wires placed at small angles to the principal axis. 
By choosing low-mass wires, and a helium-based 
gas mixture the multiple scattering inside the 
DCH is minimized. The readout electronics are 
mounted on the backward cndplatc of the c;liam- 
ber, minimizing the amount of material in front 
of the calorimeter endcap. 

The DIRC, the detector of internally reflected 
Cherenkov light, is a novel device providing sep- 
aration of pions and kaons from about 500 MeV/ c 
to the kinematic limit of 4.5GeV/c. Cherenkov 
light is produced in 4.9 m long bars of synthetic 
fused silica of rectangular cross section, 1.7 cm x 
3.5 cm, and transported by total internal reflec- 
tion, preserving the angle of emission, to an ar- 
ray of photomultiplier tubes. This array forms 
the backward wall of a toroidal water tank that is 
located beyond the backward end of the magnet. 
Images of the Cherenkov rings are reconstructed 
from the position and time of arrival of the signals 
in the photomultiplier tubes. 

The electromagnetic calorimeter (EMC) is de- 
signed to detect electromagnetic showers with ex- 
cellent energy and angular resolution over the en- 
ergy range from 20 MeV to 4 GeV. This coverage 
allows the detection of low energy tt^s and rj'^s 
from B decays and higher energy photons and 
electrons from electromagnetic, weak, and radia- 
tive processes. The EMC is a finely segmented 
array of projective geometry made of thallium 
doped cesium iodide (CsI(Tl)) crystals. The crys- 
tals are arranged in modules that are supported 
individually from an external support structure. 
This structure is built in two sections, a barrel 
and a forward endcap. To obtain the desired res- 
olution, the amount of material in front of and in- 



between the crystals is held to a minimum. The 
individual crystals are read out by pairs of silicon 
PIN diodes. Low noise analog circuits and fre- 
quent, precise calibration of the electronics and 
energy response over the full dynamic range are 
crucial for maintaining the desired performance. 

The instrumented flux return (IFR) is designed 
to identify muons and to detect neutral liadrons. 
For this purpose, the magnet flux return steel in 
the barrel and the two end doors is segmented 
into layers, increasing in thickness from 2 cm on 
the inside to 10 cm at the outside. Between these 
steel absorbers, single gap resistive plate cham- 
bers (RPCs) are inserted which detect stream- 
ers from ionizing particles via external capacitive 
readout strips. There are 19 layers of RPCs in 
the barrel sectors and 18 layers in the end doors. 
Two additional cylindrical layers of RPCs with 
four readout planes are placed at a radius just 
inside the magnet cryostat to detect particles ex- 
iting the EMC. 

2.2. Electronics, Trigger, Data Acquisition 
and Online Computing 

The electronics, trigger, data acquisition, and 
online computing represent a collection of tightly 
coupled hardware and software systems. These 
systems were designed to maximize the physics 
data acceptance, maintainability, and reliabil- 
ity while managing complexity, and minimizing 
deadtime, and cost. 

• Front-End Electronics (FEE) assemblies, 
located on the detector, consist of signal 
processing and digitization electronics along 
with the data transfer via optical fiber to 
the data acquisition system. 

• A robust and flexible two-level trigger copes 
with the full beam-beam interaction rate. 

The first level. Level 1 (LI), is implemented 
in hardware, the other. Level 3 (L3), in soft- 
ware. Provision is made for an intermediate 
trigger (Level 2) should severe conditions 
require additional sophistication. 

• The Online Dataflow (ODF), handles digi- 
tized data from the FEE through the event 
building. ODF includes the fast control and 
timing system (FCTS). 
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Table 1 

Overview of the coverage, segmentation, and performance of the BABAR detector systems. The nota- 
tion (C), (F), and (B) refers to the central barrel, forward and backward components of the system, 
respectively. The detector coverage in the laboratory frame is specified in terms of the polar angles 9i 
(forward) and 62 (backward). The number of readout channels is listed. The dynamic range (resolution) 
of the FEE circuits is specified for pulse height (time) measurements by an ADC (TDC) in terms of the 
number of bits (nsec). Performance numbers are quoted for 1 GeV/c particles, except where noted. The 
performances for the SVT and DCH are quoted for a combined Kalman fit (for the definition of the track 
parameters, see Section ^) 
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• A farm of commercial Unix processors and 
associated software. Online Event Process- 
ing (OEP), provides the realtime environ- 
ment within which complete events are pro- 
cessed by the L3 trigger algorithms, partial 
event reconstruction is performed for mon- 
itoring, and event data are logged to an in- 
termediate storage. 

• Software running on a second farm of 
processors. Online Prompt Reconstruction 
(OPR), completely reconstructs all physics 
events, and performs monitoring and con- 



stants generation in near realtime. Physics 
event data are transferred to an object 
database ^ and are made available for fur- 
ther analyses. 

• An Online Run Control (ORC) system im- 
plements the logic for managing the state of 
the detector systems, starting and stopping 
runs, and performing calibrations as well as 
providing a user control interface. 

• A system to control and monitor the detec- 
tor and its support systems, Online Detec- 
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tor Control (ODC), is based upon the Ex- 
perimental Physics Industrial Control Sys- 
tem (EPICS) toolbox 0. This system in- 
cludes communication links with PEP-II. 

2.2.1. Electronics 

All BABAR detector systems share a common 
electronics architecture. Event data from the de- 
tector flows through the FEE, while monitoring 
and control signals are handled by a separate, 
parallel system. All FEE systems are mounted 
directly on the detector to optimize performance 
and to minimize the cable plant, thereby avoiding 
noise pickup and ground loops in long signal ca- 
bles. All detector systems utilize standard BABAR 
interfaces to the data acquisition electronics and 
software. 

Each FEE chain consists of an amplifier, a digi- 
tizer, a trigger latency buffer for storing data dur- 
ing the LI trigger processing, and an event buffer 
for storing data between the LI Accept and sub- 
sequent transfer to the data acquisition system 
(see Figure ^. Custom integrated circuits (ICs) 
have been developed to perform the signal pro- 
cessing. The digital LI latency buffers function 
as fixed length data pipelines managed by com- 
mon protocol signals generated by the FCTS. All 
de-randomizing event buffers function as FIFOs 
(first-in-first-out) capable of storing a fixed num- 
ber of events. During normal operation, analog 
signal processing, digitization, and data readout 
occur continuously and simultaneously. 




Digitizer 



to LI trigger 




LI trigger accept 



Figure 4. Schematic diagram of the Front-End 
Electronics (FEE) . Analog signals arrive from the 
left, proceed conditionally through the indicated 
steps and are injected into the remainder of the 
data acquisition system. 

Since many of the front-end circuits are inac- 
cessible or require significant downtime for ac- 
cess, stringent requirements were placed on reli- 



ability. Most components underwent comprehen- 
sive mean time between failure (MTBF) studies. 
All circuits underwent a burn-in procedure prior 
to installation with the goal of minimizing initial 
failure rates. 

2.2.2. Trigger 

The trigger system operates as a sequence of 
two independent stages, the second conditional 
upon the first. The LI trigger is responsible for 
interpreting incoming detector signals, recogniz- 
ing and removing beam-induced background to a 
level acceptable for the L3 software trigger which 
runs on a farm of commercial processors. 

LI consists of pipelined hardware processors 
designed to provide an output trigger rate of 
< 2 kHz. The LI trigger selection is based on data 
from DCH, EMC, and IFR. The maximum LI re- 
sponse latency for a given collision is 12 /is. Based 
on both the complete event and LI trigger infor- 
mation, the L3 software algorithms select events 
of interest which are then stored for processing. 
The L3 output rate is administratively limited to 
120 Hz so as not to overload the downstream stor- 
age and processing capacity. 

BABAR has no fast counters for triggering pur- 
poses, and bunch crossings are nearly continuous 
at a 4.2 ns spacing. Dedicated LI trigger proces- 
sors receive data continuously clocked in from the 
DCH, EMC, and IFR detector systems. These 
processors produce clocked outputs to the fast 
control system at 30 MHz, the time granularity 
of resultant LI Accept signals. The arrival of 
an LI Accept signal by the data acquisition sys- 
tem causes a portion of each system's LI latency 
buffer to be read out, ranging from about 500 ns 
for the SVT to 4-16 /is for the EMC. Absolute 
timing information for the event, i.e., associating 
an event with a particular beam crossing, is deter- 
mined offline, using DCH track segment timing, 
waveforms from the EMC, and accelerator timing 
fiducials. 

2.2.3. Data Acquisition 

and Online Computing 

The data acquisition and computing systems, 
responsible for the transport of event data from 
the detector FEE to mass storage with a min- 
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imum of dead time are shown schematically in 
Figure ||. These systems also interface with the 
trigger to enable calibrations and testing. Other 
parts of these systems provide for the control and 
monitoring of the detector and supporting facili- 
ties. 

Hardware 

The data acquisition system hardware consists 
of VME crates, specialized VME-based proces- 
sors called readout modules (ROMs), the FCTS, a 
Unix processor farm, various server machines and 
an Ethernet network. A ROM consists of a Mo- 
torola MVME2306 PowerPC single board com- 
puter, event buffers, an interface to the FCTS, 
and a custom Personality Card that connects 
with the FEE circuits via 1.2Gbits/s fiber optic 
cables. The ROM provides the standard interface 
between the detector specific FEE, the FCTS, 
and the event builder. There are 157 ROMs in 
the system located in 19 physical VME crates 
divided into 24 logical crates by virtue of seg- 
mented backplanes. The FCTS system consists 
of a VME crate plus individual Fast Control Dis- 
tribution Modules in each of the data acquisition 
VME crates. The Unix processor farm consists of 
32 Sun workstations. 

The detector monitoring and control system 
consists of a standard set of components, includ- 
ing Motorola MVME177 single-board computers, 
and other VME modules. With the exception of 
the solenoid magnet, which has its own control 
and monitoring, all BABAR detector components 
use this system. 

The online computing system relies on a com- 
plex of workstation consoles and servers with 
0.8 Tbytes of attached storage, all interconnected 
with switched 100 Mbits/s and 1 Gbits/s Ethernet 
networks. Multiple 1 Gbits/s Ethernet links con- 
nect the experimental hall with the SLAC com- 
puting center. 

Online Dataflow (ODF) 

The ODF software connects, controls and man- 
ages the flow of data in the acquisition hard- 
ware with little dead time. This code is di- 
vided between embedded processors in the ROMs 



running the VxWorks ^ realtime operating sys- 
tem and Unix processors running the Solaris 
operating system. ODF provides configuration 
and readout of the FEE over fiber links to 
the ROMs; data transport, buffering, and event 
building from the ROMs to the Unix farm over 
a switched 100 Mbits/s Ethernet network; mask- 
ing and prescaling of LI triggers; and logical par- 
titioning of DAQ hardware into multiple, inde- 
pendent data acquisition systems for parallel cal- 
ibrations and diagnostics. Additional feature ex- 
traction (FEX) code in the ROMs extracts phys- 
ical signals from the raw data, performs gain and 
pedestal correction, data sparsification, and data 
formatting. Data from electronics calibrations 
are accumulated in the ROMs, channel response 
functions are evaluated, results are compared to 
reference data and subsequently applied in fea- 
ture extraction. Calibration data are stored in a 
conditions database. 

Online Event Processing (OEP) 

OEP receives and processes data from the ODF 
event builder on each of the Unix processors. 
OEP orchestrates the following tasks: L3 trigger 
algorithms; fast monitoring to assure data qual- 
ity; and logging the selected events to disk while 
merging the multiple data output streams to a 
single file. 

Online Prompt Reconstruction (OPR) 

OPR bridges the online and offline systems |^ . 
This system reads raw data recorded to disk by 
OEP and, operating on a farm of 150 Unix pro- 
cessors, selects physics events, performs complete 
reconstruction, performs rolling calibrations, col- 
lects extensive monitoring data for quality assur- 
ance, and writes the result into an event store. 
A rolling calibration is the generation of recon- 
struction constants during normal event process- 
ing, which are then applied to the processing of 
subsequent data. 

Online Detector Control (ODC) 
and Run Control (ORC) 

The ODC system controls and extensively 
monitors the electronics, the environment, and 
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Figure 5. Schematic diagram of the data acquisition. 



assures the safety of the detector. Its implementa- 
tion is based on EPICS, providing detector-wide 
standardization for control and monitoring, diag- 
nostics and alarm handling. ODC also provides 
communication with PEP-II and the magnet con- 
trol system. Monitoring data are archived in an 
ambient database. 

The ORG system ties together the various com- 
ponents of the online system and provides the op- 
erator with a single graphical interface to control 
detector operation. Complex configurations are 
stored in a configuration database; keys to the 
configuration used for any run are stored along 
with the data. The event store, conditions, ambi- 
ent, and configuration databases are implemented 
in an object database while other data are 
stored in a relational database. 

3. The PEP II Storage Rings and Their 
Impact on the JBABAR Detector 

3.1. PEP-II Storage Rings 

PEP-II is an e+e^ storage ring system designed 
to operate at a center of mass (cm.) energy 
of 10.58 GeV, corresponding to the mass of the 
T{AS) resonance. The parameters of these en- 
ergy asymmetric storage rings are presented in 
Table ||. PEP-II has surpassed its design goals, 
both in terms of the instantaneous and the inte- 
grated daily luminosity, with significantly fewer 
bunches than anticipated. A detailed description 



of the design and operational experience of PEP- 
II can be found in references and [pl . 

Table 2 

PEP-II beam parameters. Values are given both 
for the design and for typical colliding beam oper- 
ation in the first year. HER and LER refer to the 
high energy e~ and low energy ring, respec- 
tively. aLx, (^Ly, and glz refer to the horizontal, 
vertical, and longitudinal rms size of the luminous 
region. 



Parameters 


Design 


Typical 


Energy HER/LER (GeV) 


9.0/3.1 


9.0/3.1 


Current HER/LER (A) 


0.75/2.15 


0.7/1.3 


# of bunches 


1658 


553-829 


Bunch spacing (ns) 


4.2 


6.3-10.5 


(TLx (^ni) 


110 


120 


(JLy {^im) 


3.3 


5.6 


<JLz (mm) 


9 


9 


Luminosity (10^^ cm~^s~^) 


3 


2.5 


Luminosity (pb~^/d) 


135 


120 



PEP-II typically operates on a 40-50 minute 
fill cycle. At the end of each fill, it takes about 
three minutes to replenish the beams. After a 
loss of the stored beams, the beams are refilled in 
approximately 10-15 minutes. BABAR divides the 
data into runs, defined as periods of three hour 
duration or less during which beam and detector 
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conditions are judged to be stable. While most 
of the data are recorded at the peak of the T{AS) 
resonance, about 12% are taken at a cm. energy 
40 MeV lower to allow for studies of non-resonant 
background. 

3.2. Impact of PEP-II on BABAR Layout 

The high beam currents and the large num- 
ber of closely-spaced bunches required to produce 
the high luminosity of PEP-II tightly couple the 
issues of detector design, interaction region lay- 
out, and remediation of machine- induced back- 
ground. The bunches collide head-on and are 
separated magnetically in the horizontal plane 
by a pair of dipole magnets (Bl), followed by 
a series of offset quadrupoles. The tapered Bl 
dipoles, located at ± 21 cm on cither side of the 
IP, and the Ql quadrupoles are permanent mag- 
nets made of samarium-cobalt placed inside the 
field of the BABAR solenoid, while the Q2, Q4, and 
Q5 quadrupoles, located outside or in the fringe 
field of the solenoid, are standard iron magnets. 
The collision axis is off-set from the z-axis of the 
BABAR detector by about 20 mrad in the horizon- 
tal plane Jl2| to minimize the perturbation of the 
beams by the solenoidal field. 

The interaction region is enclosed by a water- 
cooled beam pipe of 27.9 mm outer radius, com- 
posed of two layers of beryllium (0.83 mm and 
0.53 mm thick) with a 1.48 mm water channel be- 
tween them. To attenuate synchrotron radiation, 
the inner surface of the pipe is coated with a 4 fim 
thin layer of gold. In addition, the beam pipe 
is wrapped with 150 /im of tantalum foil on ei- 
ther side of the IP, beyond z — -f-lO.lcm and 
z = —7.9cm. The total thickness of the cen- 
tral beam pipe section at normal incidence cor- 
responds to 1.06% of a radiation length. 

The beam pipe, the permanent magnets, and 
the SVT were assembled and aligned, and then 
enclosed in a 4.5 m-long support tube which spans 
the IP. The central section of this tube was fabri- 
cated from a carbon-fiber epoxy composite with 
a thickness of 0.79% of a radiation length. 

3.3. Monitoring of Beam Parameters 

The beam parameters most critical for BABAR 
performance are the luminosity, the energies of 



the two beams, and the position, angles, and size 
of the hmrinous region. 

3.3.1. Luminosity 

While PEP-II measures radiative Bhabha scat- 
tering to provide a fast monitor of the relative lu- 
minosity for operations, BABAR derives the abso- 
lute luminosity offline from other QED processes, 
primarily e^e~, and fJ.~^fJ.~ pairs. The measured 
rates are consistent and stable as a function of 
time. For a data sample of Ifb^^, the statistical 
error is less than 1%. The systematic uncertainty 
on the relative changes of the luminosity is less 
than 0.5%, while the systematic error on the ab- 
solute value of the luminosity is estimated to be 
about 1.5%. This error is currently dominated by 
uncertainties in the Monte Carlo generator and 
the simulation of the detector. It is expected that 
with a better understanding of the efhciency, the 
overall systematic error on the absolute value of 
the luminosity will be significantly reduced. 

3.3.2. Beam Energies 

During operation, the mean energies of the 
two beams are calculated from the total mag- 
netic bending strength (including the effects of 
off- axis quadrupole fields, steering magnets, and 
wigglers) and the average deviations of the ac- 
celerating frequencies from their central values. 
While the systematic uncertainty in the PEP-II 
calculation of the absolute beam energies is es- 
timated to be 5-10 MeV, the relative energy set- 
ting for each beam is accurate and stable to about 
1 MeV. The rms energy spreads of the LER and 
HER beams are 2.3 MeV and 5.5 MeV, respec- 
tively. 

To ensure that data are recorded close to the 
peak of the T(4S') resonance, the observed ratio of 
BB enriched hadronic events to lepton pair pro- 
duction is monitored online. Near the peak of the 
resonance, a 2.5% change in the BB production 
rate corresponds to a 2 MeV change in the cm. 
energy, a value that is close to the tolerance to 
which the energy of PEP-II can be held. How- 
ever, a drop in the BB rate does not distinguish 
between energy settings below or above the T(45) 
peak. The sign of the energy change must be de- 
termined from other indicators. The best mon- 
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itor and absolute calibration of the cm. energy 
is derived from the measured cm. momentum 
of fully reconstructed B mesons combined with 
the known _B-meson mass. An absolute error of 
1.1 MeV is obtained for an integrated luminosity 
of lfb~^. This error is presently limited by the 
uncertainty in the _B-meson mass and by the 
detector resolution. 

The beam energies are necessary input for the 
calculation of two kinematic variables that are 
commonly used to separate signal from back- 
ground in the analysis of exclusive i?-meson de- 
cays. These variables, which make optimum use 
of the measured quantities and are largely uncor- 
related, arc Lorentz-invariants which can be eval- 
uated both in the laboratory and cm. frames. 

The first variable, AE, can be expressed in 
Lorentz invariant form as 

AE ^ (2^90 - s)/2yi, (1) 

where — 2E^^^^^ is the total energy of the 
e+e" system in the cm. frame, and qs and qo = 
{Eq,Pq) are the Lorentz vectors representing the 
momentum of the B candidate and of the e+e" 
system, go = <le+ + qe- ■ In the cm. frame, AE 
takes the familiar form 

AE ^ Eg - El^^„^, (2) 

here Eg is the reconstructed energy of the B me- 
son. The AE distribution receives a sizable con- 
tribution from the beam energy spread, but is 
generally dominated by detector resolution. 

The second variable is the energy-substituted 
mass, mES, defined as tties^ — Qb- ^^'^ lab- 
oratory frame, ttt-es can be determined from the 
measured three-momentum pB of the B candi- 
date without explicit knowledge of the masses of 
the decay products: 

mES = ^Jis/2+pB-Po)yEi~pl. (3) 

In the cm. frame {pq = 0), this variable takes the 
familiar form 

™es = ^Ell^~p*i, (4) 

where p*g is the cm. momentum of the B meson, 
derived from the momenta of its decay products. 




Figure 6. The energy-substituted mass for a sam- 
ple of 6,700 neutral B mesons reconstructed in 
the final states I?(*)-7r+, L>(*)-p+, L>(*)-af , and 
J/ijjK*^. The background is extrapolated from 
events outside the signal region. 

and the i?- meson energy is substituted by E^^^^. 
Figure ^ shows the toes distribution for a sam- 
ple of fully reconstructed B mesons. The resolu- 
tion in Toes is dominated by the spread in -Ebeam' 
^^b*ca,„ = 2.6 MeV. 

3.3.3. Beam Direction 

The direction of the beams relative to BABAR 
is measured iteratively run- by-run using e~^e~ 
e^e~ and e'^e^ iJ-^fJ-^ events. The resultant 
uncertainty in the direction of the boost from the 
laboratory to the cm. frame, /3, is about 1 mrad, 
dominated by alignment errors. This translates 
into an uncertainty of about 0.3 MeV in rngs- P 
is consistent to within 1 mrad with the orienta- 
tion of the elongated beam spot (see below) . It is 
stable to better than 1 mrad from one run to the 
next. 

3.3.4. Beam Size and Position 

The size and position of the luminous region are 
critical parameters for the decay-time-dependent 
analyses and their values are monitored contin- 
uously online and offline The design values for 
the size of the luminous region are presented in 
Table |[ The vertical size is too small to be mea- 
sured directly. It is inferred from the measured 
luminosity, the horizontal size, and the beam cur- 
rents; it varies typically by 1-2 /im. 
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The transverse position, size, and angles of the 
luminous region relative to the BABAR coordinate 
system are determined by analyzing the distribu- 
tion of the distance of closest approach to the 
z-axis of the tracks in well measured two-track 
events as a function of the azimuth cf). The lon- 
gitudinal parameters are derived from the lon- 
gitudinal vertex distribution of the two tracks. 
A combined fit to nine parameters (three aver- 
age coordinates, three widths, and three small 
angles) converges readily, even after significant 
changes in the beam position. The uncertain- 
ties in the average beam position are of the order 
of a few /im in the transverse plane and 100 /im 
along the collision axis. Run-by-run variations in 
the beam position are comparable to these mea- 
surement uncertainties, indicating that the beams 
are stable over the period of a typical run. The 
fit parameters are stored run-by-run in the con- 
ditions database. These measurements are also 
checked offline by measuring the primary vertices 
in multi-hadron events. The measured horizontal 
and longitudinal beam sizes, corrected for track- 
ing resolution, are consistent with those measured 
by PEP-II. 

3.4. Beam Background Sources 

The primary sources of steady-state accelera- 
tor backgrounds are, in order of increasing im- 
portance: synchrotron radiation in the vicinity of 
the interaction region; interactions between the 
beam particles and the residual gas in either ring; 
and electromagnetic showers generated by beam- 
beam collisions p"^-p^. In addition, there are 
other background sources that fluctuate widely 
and can lead to very large instantaneous rates, 
thereby disrupting stable operation. 

3.4.1. Synchrotron Radiation 

Synchrotron radiation in nearby dipoles, the 
interaction-region quadrupole doublets and the 
Bl separation dipoles generates many kW of 
power and is potentially a severe background. 
The beam orbits, vacuum-pipe apertures and 
synchrotron-radiation masks have been designed 
such that most of these photons are channeled 
to a distant dump; the remainder are forced to 
undergo multiple scatters before they can en- 



ter the BABAR acceptance. The remaining syn- 
chrotron radiation background is dominated by 
x-rays (scattered from tungsten tips of a mask) 
generated by beam tails in the high-field region of 
the HER low-/? quadrupoles. This residual back- 
ground is relatively low and has not presented 
significant problems. 

3.4.2. Beam-Gas Scattering 

Beam-gas bremsstrahlung and Coulomb scat- 
tering off residual gas molecules cause beam par- 
ticles to escape the acceptance of the ring if their 
energy loss or scattering angle are sufficiently 
large. The intrinsic rate of these processes is pro- 
portional to the product of the beam current and 
the residual pressure (which itself increases with 
current). Their relative importance, as well as the 
resulting spatial distribution and absolute rate of 
lost particles impinging the vacuum pipe in the 
vicinity of the detector, depend on the beam op- 
tical functions, the limiting apertures, and the 
entire residual-pressure profile around the rings. 
The separation dipoles bend the energy-degraded 
particles from the two beams in opposite direc- 
tions and consequently most BABAR detector sys- 
tems exhibit occupancy peaks in the horizontal 
plane, i.e., the LER background near cj) ~ and 
HER background near = 180 . 

During the first few months of operation and 
during the first month after a local venting of 
the machine, the higher pressures lead to signifi- 
cantly enhanced background from beam-gas scat- 
tering. The situation has improved significantly 
with time due to scrubbing of the vacuum pipe 
by synchrotron radiation. Towards the end of the 
first year of data-taking, the dynamic pressure in 
both rings had dropped below the design goal, 
and the corresponding background contributions 
were much reduced. Nevertheless, beam-gas scat- 
tering remains the primary source of radiation 
damage in the SVT and the dominant source of 
background in all detectors systems, except for 
the DIRC. 

3.4.3. Luminosity Background 

Radiative Bhabha scattering results in energy- 
degraded electrons or positrons hitting aperture 
limitations within a few meters of the IP and 
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spraying BABAR with electromagnetic shower de- 
bris. This background is directly proportional 
to the instantaneous luminosity and thus is ex- 
pected to contribute an increasing fraction of the 
total background in the future. Already this is 
the dominant background in the DIRC. 

3.4.4. Background Fluctuations 

In addition to these steady-state background 
sources, there are instantaneous sources of radia- 
tion that fluctuate on diverse time scales: 

• beam losses during injection, 

• intense bursts of radiation, varying in du- 
ration from a few ms to several minutes, 
currently attributed to very small dust par- 
ticles, which become trapped in the beam, 
and 

• non-Gaussian tails from beam-beam inter- 
actions (especially of the e+ beam) that are 
highly sensitive to adjustments in collima- 
tor settings and ring tunes. 

These effects typically lead to short periods of 
high background and have resulted in a large 
number of i3ABAR-initiated beam aborts (see be- 
low). 

3.5. Radiation Protection 
and Monitoring 

A system has been developed to monitor the in- 
stantaneous and integrated radiation doses, and 
to either abort the beams or to halt or limit the 
rate of injection, if conditions become critical. In 
addition, DCH and IFR currents, as well as DIRC 
and IFR counting rates, are monitored; abnor- 
mally high rates signal critical conditions. 

Radiation monitoring and protection systems 
are installed for the SVT, the DCH electronics, 
and the EMC. The radiation doses are measured 
with silicon photodiodes. For the SVT, 12 diodes 
are arranged in three horizontal planes, at, above, 
and below the beam level, with four diodes in each 
plane, placed at z — +12.1 cm and z = —8.5 cm 
and at a radial distance of 3 cm from the beam 
line [0. The diode leakage current, after cor- 
rection for temperature and radiation damage ef- 
fects, is proportional to the dose rate. The four 
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Figure 7. The integrated radiation dose as mea- 
sured by PIN diodes located at three different 
positions, showing contributions from the HER 
{(j) = 180°), and the LER (0 = 0°) in the hori- 
zontal plane, and from both beams combined else- 
where. Also shown is the SVT radiation budget 
for the first 500 days of operation. 

diodes in the middle are exposed to about ten 
times more radiation than the others. These mid- 
plane diodes are connected to the beam abort 
system, while the remaining eight diodes at the 
top and bottom are used to monitor the radia- 
tion dose delivered to the SVT. The accuracy of 
the measured average dose rate is better than 0.5 
mRad/s. The integrated dose, as measured by 
the SVT diodes, is presented in Figure 0. 

The radiation level at the DCH and the EMC 
is more than two orders of magnitude lower than 
at the SVT. To amphfy the signal, the PIN diodes 
for the DCH and EMC are mounted on small 
CsI(Tl) crystals (with a volume of about 10 cm^). 
These silicon diodes are installed in sets of four. 
Three sets are mounted on the front face of the 
endcap calorimeter and one set on the backward 
endplate of the DCH, close to the readout elec- 
tronics. The signals of the four diodes in each 
set are summed, amplified, and fed into the ra- 
diation protection electronics. Only one of the 
three diode sets of the EMC is used at any given 
time. The DCH and the EMC use identical hard- 
ware and decision algorithms. They limit injec- 
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tion rates whenever an instantaneous dose equiv- 
alent to about 1 Rad/day is exceeded. 

The SVT employs a different strategy and cir- 
cuitry to assess whether the measured radiation 
levels merit a beam abort or a reduction in single- 
beam injection rate. Every beam dump initi- 
ated by BABAR is followed by a 10-15 minute 
period of injection with significant radiation ex- 
posure. Thus, to minimize the ratio of the in- 
tegrated radiation dose to the integrated lumi- 
nosity, it has been beneficial to tolerate transient 
high-dose events as long as the integrated dose 
remains less than the typical dose accumulated 
during injection. To differentiate between very 
high instantaneous radiation and sustained high 
dose rates, trip thresholds are enforced on two 
different time scales: an instantaneous dose of 
1 Rad accumulated over 1 ms, and an average 
of 50 mRad/s measured over a 5-minute period. 
During injection, higher thresholds are imposed, 
since an aborted injection will delay the return to 
taking data. 

Figure || shows the daily rate of beam aborts 
initiated by the SVT protection diodes during the 
year 2000. Initially, as many as 80 beam aborts 
were triggered per day, while the average for sta- 
ble operation was significantly below ten at the 
end of the run. The measures described above, 
combined with a significant reduction in large 
background fluctuations, have been very effective 
in protecting the detector against radiation dam- 
age, as well increasing the combined live time of 
the machine and detector to greater than 75%. 

3.6. Impact of Beam-Generated 
Background on BABAR 

Beam-generated backgrounds affect the detec- 
tor in multiple ways. They cause radiation dam- 
age to the detector components and the electron- 
ics and thus may limit the lifetime of the ex- 
periment. They may also cause electrical break- 
down and damage or generate large numbers of 
extraneous signals leading to problems with band- 
width limitations of the data acquisition system 
and with event reconstruction. Backgrounds can 
degrade resolution and decrease efficiency. 

The impact of the beam-generated background 
on the lifetime and on the operation of the dif- 
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Figure 8. Daily rates of beam aborts initiated 
by the SVT radiation protection diodes, summed 
over regular data-taking and PEP-II injection. 

ferent detector systems varies significantly. Ta- 
ble ^ lists the limits on the instantaneous and 
integrated background levels in terms of the to- 
tal dose and instantaneous observables. These 
limits are estimates derived from beam tests and 
experience of earlier experiments. For each de- 
tector system, an annual radiation allowance has 
been established taking into account the total es- 
timated lifetime of the components and the ex- 
pected annual operating conditions. The typical 
values accumulated for the first year of operation 
are also presented in the table. 

Systematic studies of background rates were 
performed with stable stored beams. Measure- 
ments of the current-dependence of the back- 
grounds were carried out for single beams, two 
beams not colliding, and colliding beams with the 
goal to identify the principal background sources, 
to develop schemes of reducing these sources, and 
to extrapolate to operation at higher luminosity 
||l6|] . These experimental studies were comple- 
mented by Monte Carlo simulations of beam-gas 
scattering and of the propagation of showers in 
the detector. The studies show that the rela- 
tive importance of the single-beam and luminos- 
ity background contributions varies, as illustrated 
in Figure ^ Data for the IFR are not shown be- 
cause this system is largely insensitive to beam- 
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Table 3 

BABAR background tolerance. Operational limits are expressed either as lifetime limits (radiation-damage 
and aging-related quantities), or in terms of instantaneous observables (DCH current, DIRC and Ll- 
trigger rates). 





Limiting factor 


Operational 


First-year 


Detector system 


and impact 


limit 


typical 


SVT sensors 


Integrated dose: 


2 MRad 


0.33 MRad 


and electronics 


radiation damage 




(hor. -plane modules) 








u.uo ivixiaci 








(other modules) 


SVT sensors 


Instantaneous dose: 


1 Rad/ms 


N/A 




diode shorts 






DCH: electronics 


Integrated dose: 


20 kRad 


< 100 Rad 




radiation damage 






DCH: wire current 


Accumulated charge: 


100 mC/cm 


8 mC/cm 




wire agmg 






DCH: total current 


HV system limitations 


1000 AiA 


250 fiA 








(steady-state) 


DIRC PMTs 


Counting rate: 


200 kHz 


110 kHz (steady-state. 




TDC deadtimc 




well-shielded sector) 


EMC crystals 


Integrated dose: 


10 kRad 


0.25 kRad 




radiation damage 




(worst case) 


LI trigger 


Counting rate: 


2 kHz 


0.7 kHz 




DAQ dead time 




(steady-state) 



generated backgrounds, except for the outer layer 
of the forward endcap, due to insufficient shield- 
ing of the external beam line components. 

The experience of the first year of operation 
and the concern for future operation for each of 
the detectors are summarized as follows. 

SVT: The most significant concern for the 
SVT with regard to machine background is the 
integrated radiation dose. The instantaneous and 
integrated dose rates in the radiation protection 
diodes are representative, to within about a fac- 
tor of two, of the radiation doses absorbed by 
the SVT modules. The exposure in the horizon- 
tal planes is an order of magnitude larger than 
elsewhere, averaging 15-25 mRad/s during sta- 
ble beam operation. The highest integrated dose 
is 450 kRad, roughly 1 kRad/day. This dose is 
about 30% below the allowance, giving some con- 



fidence that the SVT can sustain operation for 
several more years (see Figure |^) . 

DCH: For the DCH, the currents on the wires 
are the main concern, both because of the lim- 
ited capacity of the HV power supplies and the 
effect of wire aging. The currents drawn are ap- 
proximately uniformly distributed among the 44 
HV supplies, one for each quadrant of superlay- 
ers 2-10, and two per quadrant for superlayer 1. 
Consequently, the total current limit is close to 
the sum of the limits of the individual supplies. 
During stable operation the total chamber cur- 
rent is 200-300 /iA. However, radiation spikes can 
lead to currents that occasionally exceed the limit 
of 1000 /iA, causing HV supplies to trip. Other 
background effects are measured to be well be- 
low the estimated lifetime limits and thus are not 
a serious issue at this time. The average wire 
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Figure 9. Fractional steady-state background 
contributions in BABAR detector systems, mea- 
sured for single beams and colliding beams under 
stable conditions (/+ 1.25 A, I~ = 0.75 A, 
L = 2.3 X 10^3 cm-2 s'^) in July 2000. The con- 
tributions are derived from the measured doses in 
the horizontal plane for the SVT, the total cur- 
rents in the DCH, the rates in the DIRC photo- 
multipliers, the occupancy and number of pho- 
tons above 10 MeV in the EMC, and the LI trig- 
ger rates. 

occupancy has not exceeded 1-2% during stable 
operation, but the extrapolation to future opera- 
tion at higher luminosity and currents remains a 
major concern. 

DIRC: The DIRC radiators, made of syn- 
thetic fused sihca, were tested up to doses of 
100 kRad without showing any measurable effects 
and thus radiation damage is not a concern. The 
present operational limit of the DIRC is set by 
the TDC electronics which induce significant dead 
time at rates above 250 kHz, well above the sta- 
ble beam rate of 110 kHz in well shielded areas. 
Roughly half of the present rate is luminosity- 
related and can be attributed to radiative Bhabha 
scattering. The counting rate is due to debris 
from electromagnetic showers entering the water- 
filled stand- ojf box. Efforts are underway to im- 
prove the shielding of the beam pipe nearby. 

EMC: The hfetime of the EMC is set by the 
reduction in light collection in the Csl crystals 
due to radiation damage. The cumulative dose 
absorbed by the EMC is measured by a large set 
of RadFETs placed in front of the barrel and end- 



cap crystals. RadFETs [Q are realtime integrat- 
ing dosimeters based on solid-state Metal Oxide 
Semiconductor (MOS) technology. The absorbed 
dose increases approximately linearly with the in- 
tegrated luminosity. The highest dose to date 
was observed in the innermost ring of the end- 
cap, close to 250 Rad, while the barrel crystals 
accumulated about 80 Rad. The observed reduc- 
tion in light collection of 10-15% in the worst 
place, and 4-7% in the barrel, is consistent with 
expectation (see Section ^ . 

The energy resolution is dependent on the sin- 
gle crystal readout threshold, currently set at 
1 MeV. During stable beam conditions the aver- 
age crystal occupancy for random triggers is 16%, 
with 10% originating from electronics noise in the 
absence of any energy deposition. The spectrum 
of photons observed in the EMC from the LER 
and HER is presented in Figure |l^. The HER 
produces a somewhat harder spectrum. The av- 
erage occupancy for a threshold of 1 MeV and the 
average number of crystals with a deposited en- 
ergy of more than 10 MeV are shown in Figure 
pj] as a function of beam currents for both single 
and colliding beams. The occupancy increases 
significantly at smaller polar angles, in the for- 
ward endcap and the backward barrel sections, 
and in the horizontal plane. The rate increase is 
approximately linear with the single beam cur- 
rents. Background rates recorded with separated 
beams are consistent with those produced by sin- 
gle beams. For colliding beams, there is an addi- 
tional flux of photons originating from small angle 
radiative Bhabha scattering. This effect is larger 
for low energy photons and thus it is expected 
that at higher luminosities the low energy back- 
ground will raise the occupancy and thereby limit 
the EMC energy resolution. 

LI Trigger: During stable beam operation, 
the typical LI trigger rate is below 1 kHz, more 
than a factor of two below the data acquisition 
bandwidth limit of about 2.0-2.5 kHz. Experi- 
ence shows that background bursts and other rate 
spikes can raise the data volume by as much as a 
factor of two and thus it is necessary to aim for 
steady state rates significantly below the stated 
limit. For the LI trigger, the dominant sources 
of DCH triggers are particles generated by inter- 
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Figure 10. Tlie energy spectrum of pliotons 
recorded in tlie EMC by random triggers witli 
single beams at typical operating currents, LER 
at l.lA and HER at 0.7A. The electronic noise 
has been subtracted. 

actions in vacuum flanges and the Bl magnets 
(see Figure ^ in Section |ll|). This effect is most 
pronounced in the horizontal plane. At present, 
the HER background is twice as high as that of 
the LER, and the colliding beams contribute less 
than half of the combined LER and HER single 
beam triggers. 

3.7. Summary and Outlook 

Towards the end of the first year of data-taking, 
PEP-II routinely delivered beams close to de- 
sign luminosity. Due to the very close coop- 
eration with the PEP-II operations team, the 
machine-induced backgrounds have not been a 
major problem once stable conditions were estab- 
lished. The background monitoring and protec- 
tion system has become a reliable and useful tool 
to safeguard the detector operation. 

Currently planned upgrades are expected to 
raise the luminosity to 1.5 x lO'^'^cm"^ s~^ within 
a few years. The single beam backgrounds will 
increase with beam currents and the luminosity 
background is projected to exceed, or at best 
remain comparable to, the beam-gas contribu- 
tion. Measures are being prepared to reduce 
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Figure 11. Average rates in the EMC for ran- 
dom triggers as a function of the HER current for 
a fixed LER current of I.IA, both for separated 
and colliding beams; a) the single crystal occu- 
pancy for thresholds of 1 MeV and b) the number 
of crystals with a deposited energy greater than 
10 MeV. The solid curves represent a fit to the 
colliding beam data, the dashed curves indicate 
the sum of rates recorded for single beams. 

the sources and the impact of machine-related 
background on BABAR, among them upgrades to 
the DCH power supply system and to the DIRC 
TDC electronics, the addition of localized shield- 
ing against shower debris (especially for the DIRC 
stand-off box), new vacuum chambers, adjustable 
collimators, and additional pumping capacity in 
critical regions upbeam of the interaction point. 
With the expected increase in LER current 
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Figure 12. The LI trigger rate as a function of 
the HER current for single beam only, for both 
beams, separated and colliding (with a LER cur- 
rent of 1.1 A). 

and in luminosity, both the single-beam and 
the luminosity-generated LI trigger rates will in- 
crease and are projected to exceed 2 kHz (see Fig- 
ure |2|). Therefore, the DCH trigger is being up- 
graded to improve the rejection of background 
tracks originating from outside the luminous re- 
gion. In addition, the data acquisition and data 
processing capacity will need to be expanded to 
meet the demands of higher luminosity. 

Overall, the occupancy in all systems, except 
the IFR, will probably reach levels that are likely 
to impact the resolution and reconstruction efh- 
ciency. For instance, the occupancy in the EMC 
is expected to more than double. Thus, be- 
yond the relatively straight forward measures cur- 
rently planned for BABAR system upgrades, de- 
tailed studies of the impact of higher occupancy 
will be necessary for all systems. 

4. The Solenoid Magnet and Flux Return 

4.1. Overview^ 

The BABAR magnet system consists of a super- 
conducting solenoid [ p^ , a segmented flux return 
and a field compensating or bucking coil. This 
system provides the magnetic field which enables 
charged particle momentum measurement, serves 
as the hadron absorber for hadron/muon separa- 
tion, and provides the overall structure and sup- 



port for the detector components. Figures^ and ^ 
show key components of the BABAR magnet sys- 
tem and some of the nearby PEP-II magnets. 

The magnet coil cryostat is mounted inside the 
hexagonal barrel flux return by four brackets on 
each end. The flux return end doors are each 
split vertically and mounted on skids that rest 
on the floor. To permit access to the inner detec- 
tor, the doors can be raised and moved on rollers. 
At the interface between the barrel and the end 
doors, approximately 60% of the area is occupied 
by structural steel and filler plates; the remaining 
space is reserved for cables and utilities from the 
inner detectors. A vertical, triangular chase cut 
into the backward end doors contains the cryostat 
chimney. Table ^ lists the principal parameters of 
the magnet system. The total weight of the flux 
return is approximately 870 metric tons. 

To optimize the detector acceptance for un- 
equal beam energies, the center of the BABAR de- 
tector is offset by 370 mm in the electron beam 
direction. The principal component of the mag- 
netic field, Bz, lies along the +z axis; this is also 
the approximate direction of the electron beam. 
The backward end door is tailored to accommo- 
date the DIRC bar boxes and to allow access 
to the drift chamber electronics. Both ends al- 
low space and adequate shielding for the PEP-II 
quadrupoles. 

4.2. Magnetic Field Requirements 

and Design 
4.2.1. Field Requirements 

A solenoid magnetic field of 1.5 T was specified 
in order to achieve the desired momentum res- 
olution for charged particles. To simplify track 
finding and fast and accurate track fitting, the 
magnitude of the magnetic field within the track- 
ing volume was required to be constant within a 
few percent. 

The magnet was designed to minimize dis- 
turbance of operation of the PEP-II beam el- 
ements. The samarium-cobalt Bl dipole and 
Ql quadrupole magnets are located inside the 
solenoid as shown in Figure |l|. Although these 
magnets can sustain the high longitudinal field of 
1.5 T, they cannot tolerate a large radial compo- 
nent. Specifically, the field cannot exceed 0.25 T 
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Table 4 

Magnet Parameters 



Field Parameters 



Central Field 


1.5 


T 


Max. Radial Field 


<0.25 


T 


at Ql and r ~ 200 mm 






Leakage into PEP-II 


<0.01 


T 


Stored Energy 


27 


MJ 


Steel Parameters 


Overall Barrel Length 


4050 


mm 


Overall Door Thickness 


1149 


mm 


(incl. gaps for RPCs) 






Overall Height 


6545 


mm 


Plates in Barrel 


18 




9 


20 


mm 


4 


30 


mm 


3 


50 


mm 


2 


100 


mm 


Plates in Each Door 


18 




9 


20 


mm 


4 


30 


mm 


4 


50 


mm 


1 


100 


mm 


Main Coil Parameters 


Mean Diameter of 


3060 


mm 


Current Sheet 






Current Sheet Length 


3513 


mm 


Number of layers 


2 




Operating Current 


4596 


A 


Conductor Current 


1.2 


kA/ mm-^ 


Density 






Inductance 


2.57 


H 


Bucking Coil Parameters 


Inner Diameter 


1906 


mm 


Operating Current 


200 


A 


Number of Turns 


140 




Cryostat Parameters 


Inner Diameter 


1420 


mm 


Radial Thickness 


350 


mm 


Total Length 


3850 


mm 


Total Material (Al) 


- 126 


mm 



at a radius r = 200 mm (assuming a linear depen- 
dence of Br on r) without degrading their field 
properties due to partial demagnetization. The 
conventional iron quadrupoles Q2, Q4, and Q5 
are exposed to the solenoid stray fields. To avoid 
excessive induced skew octupole components, the 
stray field leaking into these beam elements is re- 
quired to be less than 0.01 T averaged over their 
apertures. 

4.2.2. Field Design Considerations 

Saturation of the steel near the coil and the gap 
between the coil and the steel leads to field non- 
uniformities. To control these non-uniformities, 
the current density of the coil is increased at the 
ends relative to the center by reducing the thick- 
ness of the aluminum stabilizer. While the re- 
quirements on the radial field component at Ql 
inside the solenoid can be satisfied easily at the 
forward end, the shape of the backward plug had 
to be specifically designed to simultaneously con- 
trol field uniformity and unwanted radial compo- 
nents. 

Leakage of magnetic flux is a problem, in par- 
ticular at the backward end. A bucking coil, 
mounted at the face of the backward door and 
surrounding the DIRC strong support tube, is 
designed to reduce the stray field to an accept- 
able level for the DIRC photomultipliers and the 
PEP-II quadrupoles. 

4.2.3. Magnetic Modeling 

Extensive calculations of the magnetic field 
were performed to develop the detailed design of 
the flux return, the solenoid coil, and the bucking 
coil. To crosscheck the results of these calcula- 
tions the fields were modeled in detail in two and 
three dimensions using commercial software [ pO[ . 
The shape of the hole in each end door was de- 
signed by optimizing various parameters, such as 
the minimum steel thickness in areas of satura- 
tion. The design of the hole in the forward door 
was particularly delicate because the highly sat- 
urated steel is very close to the Q2 quadrupole. 
The multiple finger design of the hole was chosen 
to control the saturation of the steel. 

Most of the design work was performed in two 
dimensions, but some three dimensional calcula- 
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tions were necessary to assure the accuracy of 
modeling the transitions between the end doors 
and the barrel, the leakage of field into the PEP- 
II magnets, and the impact of that leakage on 
the multipole purity ||2l|,|2^. The computations 
of the leakage field were done for central field of 
1.7 T instead of 1.5 T to provide some insurance 
against uncertainties in the modeling of complex 
steel shapes and the possible variations of the 
magnetic properties of the steel. 

4.3. Steel Flux Return 

4.3.1. Mechanical 

and Magnetic Forces 

The magnet flux return supports the detector 
components on the inside, but this load is not 
a major issue. Far greater demands are placed 
on the structural design by the magnetic forces 
and the mechanical forces from a potential earth- 
quake. 

Magnetic forces are of three kinds. First, there 
is a symmetric magnetic force on the end doors 
which was taken into consideration in their de- 
sign and construction. Second, there is an ax- 
ial force on the solenoid due to the forward- 
backward asymmetry of the steel. Because the 
steel is highly saturated in places, the magnitude 
of the field asymmetry changes when the current 
is raised from zero, and there is no position of 
the solenoid at which the force remains zero at all 
currents. Because it is important that this axial 
force should not change sign, which could cause 
a quench, the superconducting solenoid was de- 
liberately offset by 30 mm towards the forward 
door. This offset was chosen to accommodate a 
worst case scenario, including uncertainties in the 
calculation. Third, during a quench of the super- 
conducting coil, eddy currents in conducting com- 
ponents inside the magnetic volume could gen- 
erate sizable forces. These forces were analyzed 
for components such as the endplates of the drift 
chamber and the electromagnetic calorimeter and 
were found not to be a problem. 

4.3.2. Earthquake Considerations 

Because SLAG is located in an earthquake 
zone, considerable attention has been given to 
protecting the detector against severe damage 



from such an event. The entire detector is sup- 
ported on four earthquake isolators, one at each 
corner, which limit the component acceleration in 
the horizontal plane to 0.4 g. However, these iso- 
lators offer no protection in the vertical direction. 
Vertical ground accelerations of 0.6 g are consid- 
ered credible and actual component accelerations 
may be considerably larger due to resonances. By 
taking into account resonant frequencies and the 
expected frequency spectra of earthquakes, the 
magnet and all detector components have been 
designed to survive these accelerations without 
serious damage. Because the magnet is isolated 
from the ground moving beneath it, worst case 
clearances to external components, e.(7.,PEP-II 
components, are provided. It is expected that 
even during a major earthquake, damage would 
be modest. 

4.3.3. Fabrication 

The flux return was fabricated |p3| from draw- 
ings prepared by the BABAR engineering team. A 
primary concern was the magnetic properties of 
the steel. The need for a high saturation field 
dictated the choice of a low carbon steel, speci- 
fied by its chemical content (close to AISI 1006). 
The manufacturer supplied sample steel for crit- 
ical magnetic measurements and approval. The 
availability of very large tools at the factory made 
it possible to machine the entire face of each end 
of the assembled barrel, thus assuring a good fit 
of the end doors. The entire fiux return was as- 
sembled at the factory, measured mechanically, 
and inspected before disassembly for shipment. 

4.4. Magnet Coils 

The design of the superconducting solenoid 
is conservative and follows previous experience. 
The superconducting material is composed of 
niobium-titanium (46.5% by weight Nb) fila- 
ments, each less than 40 /im in diameter. These 
filaments are then wound into 0.8 mm strands, 
16 of which are then formed into Rutherford 
cable measuring 1.4x6. 4 mm. The final con- 
ductor consists of Rutherford cable co- 
extruded with pure aluminum stabilizer measur- 
ing 4.93x20.0 mm for use on the outer, high 
current density portion of the solenoid, and 
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Figure 13. A portion of cryostat assembly. The 
forward end is shown. Legend: (A) evacuated 
spaces fihed with IR-reflective insulator; (B) su- 
perconducting coil (2-layers); (C) aluminum sup- 
port cylinder; (D) aluminmii heat shield; (E) alu- 
minum cryostat housing. 

8.49x20.0 mm for the central, lower current den- 
sity portion. The conductor is covered in an insu- 
lating dry wrap fiberglass cloth which is vacuum 
impregnated with cpoxy. The conductor has a 
total length of 10.3 km. 

The solenoid is indirectly cooled to an oper- 
ating point of 4.5K using a thermo-syphon tech- 
nique. Liquid helium P5| is circulated in chan- 
nels welded to the solenoid support cylinder. Liq- 
uid helium and cold gas circulate between the 
solenoid, its shields, the liquefier-refrigerator and 
a 4000 £ storage dewar via 60 m of coaxial, gas- 
screened, flexible transfer line. The solenoid coil 
and its cryostat were fabricated |2^ to drawings 
prepared by the BABAR engineering team. Before 
shipment the fully assembled solenoid was 
cooled to operating temperature and tested with 
currents of up to 1000 A, limited by coil forces in 
the absence of the iron flux return. 

A portion of the cryostat assembly, containing 
the solenoid coil, its support cylinder and heat 
shield, is shown in Figure 

To reduce the leakage fields into the PEP-II 
components and the DIRC photomultiplicrs, an 
additional external bucking coil is installed [ p8[ . 
This is a conventional water cooled copper coil 
consisting of ten layers. Although the nominal 
operating current is 200 A, a current of up to 
575 A is attainable, if needed, to demagnetize the 
DIRC shield. 



To optimally control the stray fields and avoid 
a magnetization of the DIRC magnetic shield, the 
currents in the solenoid and the bucking coil are 
ramped together under computer control. High 
precision transducers are used to measure the 
currents and provide the feedback signals to the 
power supplies. The values of the currents are 
recorded in the BABAR database. 

4.5. Magnetic Field Map 

The goal of the magnetic field mapping and 
subsequent corrections was to determine the mag- 
netic field in the tracking volume to a precision 
of 0.2mT. 

4.5.1. Mapping Procedure 

A field mapping device was built specifically for 
the BABAR magnet based on a design concept de- 
veloped at Fermilab j29j . The magnetic field sen- 
sors were mounted on a propeller at the end of a 
long cantilevered spindle which reached through 
the hole in the forward end door. The spindle in 
turn rode on a carriage which moved on precision- 
aligned rails. The propeller rotated to sample the 
magnetic volume in 0, and the carriage moved 
along its axis to cover z. Measurements were ob- 
tained from five sets of Br and Bz and two B^ 
Hall probes, all of which were mounted on a plate 
at different radial positions. This plate was at- 
tached to the propeller and its position could be 
changed to cover the desired range in the radial 
distance r from the axis. Precision optical align- 
ment tools were used to determine the position of 
the sensors transverse to the z-axis. 

The Br and Bz probes were two-element de- 
vices with a short-term (few month) precision of 
0.01%, the B^ probes were single element devices 
with a precision of 0.1% |3^. In addition to the 
Hall probes, an NMR probe |^l[] was mounted at 
a radius of 89 mm on the propeller to provide a 
very precise field reference near the 2-axis as a 
function of z for \z\ < 1000mm, where z = at 
the magnet center. The NMR measurements set 
the absolute scale of the magnetic field. 

The magnetic field was mapped at the nomi- 
nal central operating field of 1.5 T, as well as at 
1.0 T. Measurements were recorded in 100 mm 
intervals from —1800 to +1800 mm in z, and in 
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24 azimuthal positions spaced by 15° for each of 
three different radial positions of the HaU probe 
plate. Thus for each z and position, the com- 
ponents and were measured at 13 distinct 
radii from 130 mm to 1255 mm and B^ at six radii 
between 505 mm and 1180 mm. 

The field map was parameterized in terms of a 
polynomial of degree up to 40 in r and z plus ad- 
ditional terms to account for expected perturba- 
tions The fit reproduced the measurements 
to within an average deviation of 0.2 mT through- 
out the tracking volume. The fitting procedure 
also served as a means of detecting and removing 
questionable measurements. 

4.5.2. Perturbations to the Field Map 

During the mapping process, the permanent 
magnet dipoles (Bl) and quadrupoles (Ql) were 
not yet installed. Their presence inside the 
solenoid results in field perturbations of two 
kinds. The first is due to the fringe fields of the 
Bl and Ql permanent magnets, and of the dipole 
and quadrupole trim coils mounted on Ql. The 
Bl field strength reaches a maximum of ~20mT 
close to the surface of the Bl casing and decreases 
rapidly with increasing radius. The fields associ- 
ated with the trim coils were measured and pa- 
rameterized prior to installation; they are essen- 
tially dipole in character. 

The second field perturbation is due to the per- 
meability of the permanent magnet material. Sin- 
tered samarium-cobalt has a relative permeabil- 
ity of 1.11 to 1.13 in the z direction, and as a 
result the solenoid field is modified significantly. 
Probes between the Bl and Ql magnets at a ra- 
dius of about 190 mm measure the effect of the 
permeability. The field perturbation is obtained 
from a two-dimensional, finite element analysis 
which reproduces the r and z dependence of B^ 
and -B^ . The induced magnetization increases B^ 
by about 9mT at the interaction point; the effect 
decreases slowly with increasing radius. 

4.5.3. Field Quality 

To illustrate the quality of magnetic field. Fig- 
ure ^ shows the field components i?^ and B^ as 
a function of z for various radial distances r. In 
the tracking volume the field is very uniform, the 
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Figure 14. The magnetic field components B^ 
and Br as a function of z for various radial dis- 
tances r (in m). The extent of the DCH and 
the location of the interaction point (IP) are in- 
dicated. 

B^ component does not exceed ImT. The vari- 
ation of the bend field, i.e., the field transverse 
to the trajectory, along the path of a high mo- 
mentum track is at most 2.5% from maximum to 
minimum within the tracking region, as shown in 
Figure |l5|. 

4.5.4. Field Computation 

In order to reduce the computation of the mag- 
netic field for track reconstruction and momen- 
tum determination, the field values averaged over 
azimuth are stored in a grid of r~z space points 
spanning the volume interior to the cryostat. Lo- 
cal values are obtained by interpolation. Within 
the volume of the SVT, a linear interpolation is 
performed in a 20 mm grid; elsewhere the interpo- 
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Figure 15. Relative magnitude of magnetic field 
transverse to a high momentum track as a func- 
tion of track length from the IP for various polar 
angles (in degrees). The data are normalized to 
the field at the origin. 

lation is quadratic in a 50 mm grid. Azimuthal de- 
pendence is parameterized by means of a Fourier 
expansion at each r-z point. The Fourier coef- 
ficients at the point of interest are obtained by 
interpolation on the r-z grid, and the average 
field value is corrected using the resulting Fourier 
series. 

4.6. Summary 

Since its successful commissioning, the magnet 
system has performed without problems. There 
have been no spontaneous quenches of the su- 
perconducting solenoid. In the tracking region, 
the magnetic field meets specifications, both in 
magnitude and uniformity. The field compen- 
sation and magnetic shielding work well for the 
DIRC photomultiplier array and the external 
quadrupoles. Measurements indicate that the 
bucking coil reduces the field at the face of Q2 
from '--^SOmT to ^ImT ||2^, in agreement with 
calculations. 

5. Silicon Vertex Tracker 

5.1. Charged Particle Tracking 

The principal purpose of the BABAR charged 
particle tracking systems, the SVT and the DCH, 



is the efficient detection of charged particles and 
the measurement of their momentum and angles 
with high precision. Among many applications, 
these precision measurements allow for the re- 
construction of exclusive B- and D-meson de- 
cays with high resolution and thus minimal back- 
ground. The reconstruction of multiple decay ver- 
tices of weakly decaying B and D mesons is of 
prime importance to the physics goals. 

Track measurements are also important for the 
extrapolation to the DIRC, EMC, and IFR. At 
lower momenta, the DCH measurements are more 
important, while at higher momenta the SVT 
measurements dominate. Most critical are the 
angles at the DIRC, because the uncertainties in 
the charged particle track parameters add to the 
uncertainty in the measurement of the Cherenkov 
angle. Thus, the track errors from the combined 
SVT and DCH measurements should be small 
compared to the average DIRC Cherenkov angle 
measurements, i.e., of order of Imrad, particu- 
larly at the highest momenta. 

5.2. SVT Goals and Design Requirements 

The SVT has been designed to provide pre- 
cise reconstruction of charged particle trajecto- 
ries and decay vertices near the interaction region. 
The design choices were driven primarily by di- 
rect requirements from physics measurements and 
constraints imposed by the PEP-II interaction re- 
gion and BABAR experiment. In this section the 
mechanical and electronic design of the SVT are 
discussed, with some discussion of the point res- 
olution per layer and dE/dx performance. The 
tracking performance and efhciency of the SVT 
alone and in combination with the DCH are de- 
scribed in Section 0. 

5.2.1. SVT Requirements and Constraints 

The SVT is critical for the measurement of the 
time-dependent CP asymmetry. To avoid signif- 
icant impact of the resolution on the CP asym- 
metry measurement the mean vertex resolution 
along the z-axis for a fully reconstructed B de- 
cay must be better than 80 /im . The required 
resolution in the x-y plane arises from the need 
to reconstruct final states in B decays as well as 
in T and charm decays. For example, in decays 
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of the type i?" D^D^ , separating the two D 
vertices is important. The distance between the 
two -D's in the x-y plane for this decay is typi- 
cally ~ 275 /im. Hence, the SVT needs to provide 
resolution of order ^100 /im in the plane perpen- 
dicular to the beam line. 

Many of the decay products of B mesons have 
low pt ■ The SVT must provide standalone track- 
ing for particles with transverse momentum less 
than 120MeV/c, the minimum that can be mea- 
sured reliably in the DCH alone. This feature is 
fundamental for the identification of slow pions 
from Z?*-meson decays: a tracking efficiency of 
70% or more is desirable for tracks with a trans- 
verse momentum in the range 50-120 MeV/c. The 
standalone tracking capability and the need to 
link SVT tracks to the DCH were crucial in choos- 
ing the number of layers. 

Beyond the standalone tracking capability, the 
SVT provides the best measurement of track an- 
gles, which is required to achieve design resolu- 
tion for the Cherenkov angle for high momentum 
tracks. 

Additional constraints are imposed by the stor- 
age ring components. The SVT is located inside 
the ^--^4.5 m-long support tube, that extends all 
the way through the detector. To maximize the 
angular coverage, the SVT must extend down to 
350mrad (20°) in polar angle from the beam line 
in the forward direction. The region at smaller 
polar angles is occupied by the Bl permanent 
magnets. In the backward direction, it is suffi- 
cient to extend the SVT sensitive area down to 
30°. 

The SVT must withstand 2 MRad of ionizing 
radiation. A radiation monitoring system capa- 
ble of aborting the beams is required. The ex- 
pected radiation dose is 1 Rad/day in the hori- 
zontal plane immediately outside the beam pipe 
(where the highest radiation is concentrated), and 
0.1 Rad/day on average otherwise. 

The SVT is inaccessible during normal detec- 
tor operations. Hence, reliability and robustness 
are essential: all components of the SVT inside 
the support tube should have long mean-time-to- 
failure, because the time needed for any replace- 
ment is estimated to be 4-5 months. Redundan- 
cies are built in whenever possible and practical. 




Figure 16. Fully assembled SVT. The silicon sen- 
sors of the outer layer are visible, as is the carbon- 
fiber space frame (black structure) that surrounds 
the silicon. 

The SVT is cooled to remove the heat gener- 
ated by the electronics. In addition, it operates 
in the 1.5 T magnetic field. 

To achieve the position resolution necessary to 
carry out physics analyses, the relative position 
of the individual silicon sensors should be sta- 
ble over long time periods. The assembly allows 
for relative motion of the support structures with 
respect to the Bl magnets. 

These requirements and constraints have led to 
the choice of a SVT made of five layers of double- 
sided silicon strip sensors. To fulfill the physics 
requirements, the spatial resolution, for perpen- 
dicular tracks, must be 10-15 /im in the three in- 
ner layers and about 40 //m in the two outer lay- 
ers. The inner three layers perform the impact 
parameter measurements, while the outer layers 
are necessary for pattern recognition and low pt 
tracking. 

5.3. SVT Layout 

The five layers of double-sided silicon strip sen- 
sors, which form the SVT detector, are organized 
in 6, 6, 6, 16, and 18 modules, respectively; a pho- 



tograph is shown in Figure 16, The strips on the 



opposite sides of each sensor are oriented orthog- 
onally to each other: the (j) measuring strips {(p 
strips) run parallel to the beam and the z mea- 
suring strips {z strips) are oriented transversely 
to the beam axis. The modules of the inner three 
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Figure 17. Schematic view of SVT: longitudinal section. The roman numerals label the six different types 
of sensors. 



layers are straight, while the modules of layers 4 
and 5 are arc/i-shaped (Figures |lj and p^ . 

This arch design was chosen to minimize the 
amount of silicon required to cover the solid angle, 
while increasing the crossing angle for particles 
near the edges of acceptance. A photograph of 
an outer layer arch module is shown in Figure ^ 
The modules are divided electrically into two half- 
modules, which are read out at the ends. 




Figure 18. Schematic view of SVT: tranverse sec- 
tion. 



To satisfy the different geometrical require- 
ments of the five SVT layers, five different sen- 
sor shapes are required to assemble the planar 
sections of the layers. The smallest detectors 
are 43 x 42 mm^ (z x (f>), and the largest are 
68 X 53 mm^. Two identical trapezoidal sensors 
are added (one each at the forward and back- 
ward ends) to form the arch modules. The half- 
modules are given mechanical stiffness by means 
of two carbon fiber/kevlar ribs, which are visible 
in Figure |9|. The strips of sensors in the same 
half-module are electrically connected with wire 
bonds to form a single readout strip. This results 
in a total strip length up to 140 mm (240 mm) in 
the inner (outer) layers. 

The signals from the z strips are brought to the 
readout electronics using fanout circuits consist- 
ing of conducting traces on a thin (50 fim) insu- 
lating Upilex |33|] substrate. For the innermost 
three layers, each z strip is connected to its own 
preamplifier channel, while in layers 4 and 5 the 
number of z strips on a half-module exceeds the 
number of electronics channels available, requir- 
ing that two z strips on different sensors be elec- 
trically connected (ganged) to a single electronics 
channel. The length of a z strip is about 50 mm 
(no ganging) or 100mm (two strips connected). 
The ganging introduces an ambiguity on the z 
coordinate measurement, which must be resolved 
by the pattern recognition algorithms. The to- 
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Figure 19. Photograph of an SVT arch module 
in an assembly jig. 

tal number of readout channels is approximately 
150,000. 

The inner modules are tilted in (j) by 5°, allow- 
ing an overlap region between adjacent modules, 
a feature that provides full azimuthal coverage 
and is advantageous for alignment. The outer 
modules cannot be tilted, because of the arch ge- 
ometry. To avoid gaps and to have a suitable 
overlap in the (j) coordinate, layers 4 and 5 are 
divided into two sub-layers (4a, 4b, 5a, 5b) and 
placed at slightly different radii (see Figure |l^) . 
The relevant geometrical parameters of each layer 
are summarized in Table 1|. 

In order to minimize the material in the 
acceptance region, the readout electronics are 
mounted entirely outside the active detector vol- 
ume. The forward electronics must be mounted 
in the 10 mm space between the 350mrad stay- 
clear space and Bl magnet. This implies that 
the hybrids carrying the front-end chip must be 



Table 5 

Geometric parameters for each layer and readout 
plane of the SVT. Floating strips refers to the 
number of strips between readout (R-0) strips. 
Note: parts of the (j) sides of layers 1 and 2 are 
bonded at 100 ^m and 110 /.tm pitch, respec- 
tively, with one floating strip. Strip length of z- 
strips for layers 4 and 5 includes ganging. The 
radial range for layers 4 and 5 includes the radial 
extent of the arched sections. 







R-0 




Strip 


Layer/ 


Radius 


pitch 


Floating 


length 


view 


(mm) 


(^m) 


strips 


(mm) 


1 z 


32 


100 


1 


40 
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positioned at an angle of 350mrad relative to the 
sensor for the layers 3, 4, and 5 (Figure |l7|) . In the 
backward direction, the available space is larger 
and the inner layer electronics can be placed in 
the sensor plane, allowing a simplified assembly. 

The module assembly and the mechanics are 
quite complicated, especially for the arch mod- 
ules, and are described in detail elsewhere pi[ |. 
The SVT support structure (Figure |l^) is a rigid 
body made from two carbon- fiber cones, con- 
nected by a space frame, also made of carbon- fiber 
epoxy laminate. 

An optical survey of the SVT on its assembly 
jig indicated that the global error in placement of 
the sensors with respect to design was ~200 ^m, 
FWHM. Subsequently, the detector was disas- 
sembled and shipped to SLAG, where it was re- 
assembled on the IR magnets. The SVT is at- 
tached to the Bl magnets by a set of gimbal rings 
in such a way as to allow for relative motion of 
the two Bl magnets while fixing the position of 
the SVT relative to the forward Bl and the orien- 
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tation relative to the axis of both Bl dipoles. The 
support tube structure is mounted on the PEP- 
II accelerator supports, independently of BABAR, 
allowing for movement between the SVT and the 
rest of BABAR. Precise monitoring of the beam 
interact ion point is necessary, as is described in 
Section 5.5 . 

The total active silicon area is 0.96 m^ and the 
material traversed by particles is ^ 4% of a radi- 
ation length (see Section^). The geometrical ac- 
ceptance of SVT is 90% of the solid angle in the 
cm. system, typically 86% are used in charged 
particle tracking. 

5.4. SVT Components 

A block diagram of SVT components is shown 
in Figure The basic components of the de- 
tector are the silicon sensors, the fanout circuits, 
the Front End Electronics (FEE) and the data 
transmission system. Each of these components 
is discussed below. 
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Figure 20. Schematic block diagram showing the 
different components of the SVT. 



5.4.1. Silicon Sensors 



The SVT sensors 1 35 are 300 /im thick double- 
sided silicon strip devices. They were designed 
at INFN Pisa and Trieste (Italy) and fabri- 
cated commercially [ p6[ . They are built on high- 
resistivity (6-15kf2-cm) n-type substrates with 
p^ strips and n"*" strips on the two opposite sides. 
The insulation of the n"*" strips is provided by in- 
dividual p-stops, so as to achieve an inter-strip 



resistance greater than 100 MJ7 at operating bias 
voltage, normally about 10 V above the depletion 
voltage. 

Typical depletion voltages are in the range 25- 
35 V. The strips are biased on both sides with 
polysilicon resistors (4-20 MJ7) to ensure the re- 
quired radiation hardness, keeping the voltage 
drop across resistors and the parallel noise as 
low as possible. Strips are AC-coupled to the 
electronics via integrated decoupling capacitors, 
the capacitance of which depends on the sen- 
sor shape, but is always greater than 14pF/cm. 
The sensors were designed to maximize the ac- 
tive area, which extends to within 0.7 mm of the 
physical edges. Another design goal was to con- 
trol the inter-strip capacitance: values between 
0.7pF/cm and l.lpF/cm were obtained for the 
various sensor shapes. To achieve the required 
spatial resolution, while keeping the number of 
readout channels as low as possible, most of the 
modules have a floating strip {i.e., not read out) 
between two readout strips. 

The leakage currents, because of the excellent 
performance of the manufacturing process, were 
as low as 50nA/cm'^ on average, measured at 
10 V above depletion voltage. The silicon sensor 
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parameters have been measured after irradiation 
with ^"^Co sources. Apart from an increase in the 
inter-strip capacitance of about 12% during the 
first 100 krad, the main effect was an increase of 
the leakage current by 0.7//A/cm^/MRad. How- 
ever, in a radiation test performed in a 1 GeV/c 
electron beam, an increase in leakage current of 
about 2 /iA/cm^/MRad and a significant shift in 
the depletion voltage, dependent on the initial 
dopant concentration, were observed. A shift 
of about 8-10 V was seen for irradiation corre- 
sponding to a dose of approximately 1 MRad. 
These observations indicate significant bulk dam- 
age caused by energetic electrons. As indicated 
by the change in depletion voltage, the SVT sen- 
sors could undergo type inversion after about 1~ 
3 MRad. Preliminary tests show that the sen- 
sors continue to function after type inversion ]37| . 
Studies of the behavior of SVT modules as a func- 
tion of radiation dose continue. 

5.4.2. Fanout Circuits 

The fanout circuits, which route the signals 
from the strips to the electronics, have been de- 
signed to minimize the series resistance and the 
inter-strip capacitance. As described in ref. ]38| , 
a trace on the fanout has a series resistance about 
1.6 f?/cm, an inter-strip resistance > 20 MX?, and 
an inter-strip capacitance < 0.5pF/cm. The elec- 
trical parameters of the final assembly of sensors 
and fanouts (referred to as Detector Fanout As- 
semblies or DFAs) are summarized in Table |^. 
Due to the different strip lengths, there are large 
differences between the inner and the outer lay- 
ers. Smaller differences are also present between 
the forward and backward halves of the module, 
that are of different lengths. 

5.4.3. Front End Electronics 

The electrical parameters of a DFA and the 
general BABAR requirements are the basic inputs 
that drove the design of the SVT front-end cus- 
tom IC; the ATOM {A Time- Over- Threshold Ma- 
chine). In particular, the front-end IC had to 
satisfy the following requirements: 

• signal to noise ratio greater than 15 for min- 
imum ionizing particle (MIP) signals for all 
modules; 



• signals from all strips must be retained, 
in order to improve the spatial resolution 
through interpolation, while keeping the 
number of transmitted hits as low as pos- 
sible. A hit refers to a deposited charge 
greater than 0.95 fC, corresponding to 0.25 
MIP; 

• the amplifier must be sensitive to both neg- 
ative and positive charge; 

• the peaking time must be programmable, 
with a minimum of 100 ns (in layers 1 and 
2, because of the high occupancy), up to 
400 ns (outer layers, with high capacitance); 

• capability to accept random triggers with a 
latency up to 11.5 /iS and a programmable 
jitter up to ±1 /is, without dead time; 

• radiation hardness greater than 2.5 MRad; 

• small dimensions; 128 channels in a 6.2 mm- 
wide chip. 

These requirements are fully satisfied by the 
ATOM IC [39], which is depicted schematically 



in Figure |21 




From Silicon 



Figure 21. Schematic diagram of the ATOM front 
end IC. 

The linear analog section consists of a charge- 
sensitive preamplifier followed by a shaper. Gains 
of 200mV/fC (low) or 300mV/fC (high) may 
be selected. The channel gains on a IC are 
uniform to 5mV/fC. Signals are presented to a 
programmable-threshold comparator, designed so 
that the output width of the pulse ( Time over 
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Table 7 

ATOM chip ENC parameters at different peaking 
times 



Peaking 


ENC 


Noise 


time 


(OpF) 


slope 


100 ns 


380 e- 


40.9 e-/pF 


200 ns 


280 e- 


33.9 e-/pF 


400 ns 


220 e- 


25.4e-/pF 



Threshold or ToT) is a quasi-logarithmic function 
of the collected charge. This output is sampled at 
30 MHz and stored in a 193 location buffer. Upon 
receipt of a Level 1 (LI) trigger, the time and 
ToT is retrieved from this latency buffer, sparci- 
fied, and stored in a four event buffer. Upon the 
receipt of an LI Accept command from the data 
acquisition system, the output data (the 4 bits for 
the ToT, 5 bits for the time stamp, and 7 bits for 
the strip address) are formatted, serialized, and 
delivered to the ROM. The IC also contains a test 
charge injection circuit. The typical noise behav- 
ior of the ATOM, as described by the Equivalent 
Noise Charge (ENC) of the linear analog section 
IS given m Table 0. 

The average noise for the various module types 
is shown in Table Given that shot noise due to 
sensor leakage current is negligible, the expected 
noise may be calculated from the parameters of 
Tables ^ and ^ The results of such a calculation 
are also shown in Table |^. The maximum average 
noise is 1,600 electrons, leading to a signal-to- 
noise ratio greater than 15. 

The power consumption of the IC is about 
4.5mW/channel. Radiation hardness was stud- 
ied up to 2.4 MRad with a 6°Co source. At that 
dose, the gain decreased 20%, and the noise in- 
creased less than 15%. 

The ATOM ICs are mounted on thick-film 
double-sided hybrid circuits (known as High 
Density Interconnects or HDIs) based on an 
aluminum-nitride substrate with high thermal 
conductivity. The electronics are powered 
through a floating power supply system, in such a 
way as to guarantee a small voltage drop (< 1 V) 
across the detector decoupling capacitors. 



5.4.4. Data Transmission 

The digitized signals are transmitted from the 
ATOM ICs through a thin kapton tail or ca- 
ble to the matching cards, from where they are 
routed to more conventional cables. Just outside 
the detector, signals are multiplexed by the MUX 
modules, converted into optical signals and trans- 
mitted to the Readout Modules (ROMs). The 
MUX modules also receive digital signals from the 
DAQ via a fiber optical connection. The SVT is 
connected to the BABAR online detector control 
and monitoring system via the industry standard 
CAN bus. Details on SVT data transmission sys- 
tem and DAQ can be found in references |^,^ . 
Power to SVT modules (silicon sensor bias volt- 
age and ATOM low voltages) is provided by a 
CAEN A522 power supply system |4|]. 

5.5. Monitoring and Calibration 

To identify immediately any operational prob- 
lems, the SVT is integrated in the control and 
monitoring system (see Section |l^). Major 
concerns for SVT monitoring are temperature 
and humidity, mechanical position, and radiation 
dose. 

5.5.1. Temperature 

and Humidity Monitors 

The total power dissipation of the SVT mod- 
ules is about 350 W, mainly dissipated in the 
ATOM ICs. External cooling is provided by 
chilled water at 8°C. In addition, humidity is re- 
duced by a stream of dry air in the support tube. 

Since condensation or excessive temperature 
can permanently damage the FEE, temperature 
and humidity monitoring are very important to 
the safe operation of the SVT. Thermistors are 
located on the HDIs (for the measurements of 
FEE temperature), around the SVT, along the 
cooling systems, and in the electronics (MUX) 
crates. The absolute temperatures are monitored 
to 0.2°C and relative changes of 0.1°C. Addition- 
ally, a series of humidity sensors are employed 
to monitor the performance of the dry air sys- 
tem. The temperature and humidity monitors 
also serve as an interlock to the HDI power sup- 
phes. 
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Figure 22. Horizontal motion between the DCH 
and the support tube measured with the capaci- 
tive sensors (curve) compared to the mean x coor- 
dinate of the interaction point (circles) measured 
with e'^e~ and events over a three-day pe- 

riod in July 1999. An arbitrary offset and scale 
has been applied to the beam position data. 

5.5.2. Position Monitors 

A system of capacitive sensors was installed to 
identify and track changes in the position of the 
SVT with respect to the PEP-II Bl magnets and 
the position of the support tube with respect to 
the DCH. An example of the understanding that 
can be achieved by this system is given in Fig- 
ure |2^, where the measured changes in the hor- 
izontal position of the SVT relative to the DCH 
are shown for a period of six day in the summer of 
1999. These position changes can be attributed 
to local temperature variations. The sensor data 
are compared to measurements of the mean po- 
sition of the interaction point (in the horizontal 
plane) determined with e^e~ and events 
recorded over this period. While the amplitude 
of motion at the time was uncharacteristically 
large, the strong correlation between these inde- 
pendent measurements is quite evident. Align- 
ment with charged particle tracks is now per- 
formed routinely to correct for relative motion 
of the tracking systems, as described in Section 
5.6.2. 

5.5.3. Radiation Monitors 

Radiation monitoring is extremely important 
to ensure the SVT does not exceed its radiation 
budget, which could cause permanent damage to 



the device. To date, the measured radiation ab- 
sorbed by the SVT is well within the allowed bud- 
get. 

The monitoring of radiation dose to the SVT 
is discussed in detail in Section |^. 

5.5.4. Calibrations 

Once a day, and each time the SVT configu- 
ration has changed, calibrations are performed in 
absence of circulating beams. All electronic chan- 
nels are tested with pulses through test capaci- 
tors, for different values of the injected charge. 
Gains, thresholds, and electronic noise are mea- 
sured, and defective channels are identified. The 
calibration results have proven very stable and re- 
peatable. The main variation in time is the occa- 
sional discovery of a new defective channel. The 
calibration procedures have also been very useful 
for monitoring noise sources external to the SVT. 

5.5.5. Defects 

Due to a series of minor mishaps incurred dur- 
ing the installation of the SVT, nine out of 208 
readout sections (each corresponding to one side 
of a half-module) were damaged and are cur- 
rently not functioning. There is no single fail- 
ure mode, but several causes: defective connec- 
tors, mishandling during installation, and not- 
fully-understood problems on the FEE hybrid. 
There has been no module failure due to radi- 
ation damage. It should be noted that due to 
the redundancy afforded by the five layers of the 
SVT, the presence of the defective modules has 
minimal impact on physics analyses. 

In addition, there are individual channel de- 
fects, of various types, at a level of about 1%. Cal- 
ibrations have revealed an increase in the num- 
ber of defective channels at a rate of less than 
0.2%/year. 

5.6. Data Analysis and Performance 

This section describes the reconstruction of 
space points from signals in adjacent strips on 
both sides of the sensors, the SVT internal and 
global alignment, single hit efficiency, and resolu- 
tion and dE/dx performance of the SVT. 
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5.6.1. Cluster and Hit Reconstruction 

Under normal running conditions, the average 
occupancy of the SVT in a time window of 1 /xs 
is about 3% for the inner layers, with a signif- 
icant azimuthal variation due to beam-induced 
backgrounds, and less than 1% for the outer lay- 
ers, where noise hits dominate. Figure shows 
the typical occupancy as a function of IC index 
(equivalent to azimuthal angle, in this case) for 
layer 1, cf) side. In the inner layers, the occupancy 
is dominated by machine backgrounds, which are 
significantly higher in the horizontal plane, seen 
in the plot as the peaks near IC indices 3 and 25. 




10 20 30 

IC Index 



Figure 23. Typical occupancy in percent as 
a function of IC index in layer 1, side for 
a) forward half-modules and b) backward half- 
modules. The IC index increases with azimuthal 
angle and the higher occupancy in the horizontal 
plane is visible near chip indices 3 and 25. 

The first step of the reconstruction program 
consists in discarding out-of-time channels. A 
time correction, i.e., the time between the pas- 
sage of the particle and the time the shaper ex- 
ceeds threshold, is performed, after which hits 
with times more than 200 ns from the event time 
(determined by the DCH) are discarded. The 
loss of real hits from this procedure is negligible. 
The resulting in-time hits are then passed to the 



cluster finding algorithm. First, the charge pulse 
height (Q) of a single pulse is calculated from 
the ToT value, and clusters are formed group- 
ing adjacent strips with consistent times. In a 
second pass, clusters separated by just one strip 
are merged into one cluster. The two original 
clusters plus the merged cluster are made avail- 
able to the pattern recognition algorithm, which 
chooses among the three. 

The position x oi a cluster formed by n strips 
is determined, with the "head-to-tail" algorithm: 

_ {Xi + Xn) p {Qn - Ql) 

" 2 +2 (g„ + Ql) ' 

where Xi and Qi are the position and collected 
charge of i-th. strip, respectively, and p is the read- 
out pitch. This formula results in a cluster posi- 
tion that is always within p/2 of the geometrical 
center of the cluster. The cluster pulse height 
is simply the sum of the strip charges, while the 
cluster time is the average of the signal times. 

5.6.2. Alignment 

The alignment of the SVT is performed in two 
steps. The first step consists of determining the 
relative positions of the 340 silicon sensors. Once 
this is accomplished, the next step is to align 
the SVT as a whole within the global coordi- 
nate system defined by the DCH. The primary 
reason for breaking the alignment procedure into 
these two steps is that the local positions are rela- 
tively stable in time compared to the global posi- 
tion. Also, the local alignment procedure is con- 
siderably more complex than the global alignment 
procedure. Thus, the global alignment can be up- 
dated on a run- by- run basis, while the local align- 
ment constants are changed as needed, typically 
after magnet quenches or detector access. 

The local alignment procedure is performed 
with tracks from e^e~ — s- fi^fJ-~ events and cos- 
mic rays. Well isolated, high momentum tracks 
from hadronic events are also used to supplement 
di-muon and cosmic data. Data samples suffi- 
cient to perform the local alignment are collected 
in one to two days of typical running conditions. 

In fJ-'^ fi~ events, the two tracks are simulta- 
neously fit using a Kalman filter technique and 
the known beam momentum. The use of tracks 
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from cosmic rays reduces any systematic distor- 
tion that may be introduced due to imprecise 
knowledge of the beam momenta. No informa- 
tion from the DCH is used, effectively decoupling 
the SVT and DCH alignment. 

In addition to the information from tracks, 
data from an optical survey performed during the 
assembly of the SVT are included in the align- 
ment procedure. The typical precision of these 
optical measurements is 4 /im. This survey infor- 
mation is only used to constrain sensors relative 
to other sensors in the same module, but not one 
module to another or one layer to another. Fur- 
thermore, only degrees of freedom in the plane of 
the sensor are constrained as they are expected to 
be the most stable, given the assembly procedure. 

Using the hit residuals from the aforementioned 
set of tracks and the optical survey information, a 
is formed for each sensor and minimized with 
respect to the sensor's six local parameters. The 
constraints coming from the overlapping regions 
of the silicon sensors, the di-muon fit, the cosmic 
rays, and the optical survey result in internally 
consistent local alignment constants. 

Figure ^ shows a comparison between the op- 
tical alignment made during the SVT assembly in 
February 1999 and a local alignment using data 
taken during January 2000. The alignment from 
tracking data was made without using cosmics or 
constraints from the optical survey. The width of 
the distributions in these plots has four contribu- 
tions: 1) displacement during the transfer of the 
SVT from the assembly jig to the IR magnets, 2) 
time dependent motion of the SVT after mount- 
ing, 3) statistical errors, and 4) systematic errors. 
The second set of plots shows the difference in two 
alignment sets for data taken in January 2000 as 
compared to March 2000. In general, the sta- 
bility of the inner three layers is excellent, with 
slightly larger tails in the outer two layers. The 
radial coordinate is less tightly constrained in all 
measurements because the radial location of the 
charge deposition is not well known, and most of 
the information about the radial locations comes 
only from constraints in the overlap region of the 
sensors. 

After the internal alignment, the SVT is con- 
sidered as a rigid body. The second alignment 
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Figure 24. Comparison of a local alignment of all 
the sensors in the SVT using data from January 
2000 with the optical survey of the SVT made 
during assembly in February 1999 in the (a) rA(f>, 
(b) Az and (c) Ar coordinates. Plots (d), (e), and 
(f) show the difference between two local align- 
ments using data from January 15-19 and March 
6-7, 2000 for the rAcf), Az, and Ar coordinates, 
respectively. In all the plots, the shaded regions 
correspond to the sensors in the first three layers. 
In comparing the different alignments and optical 
survey, a six parameter fit (three global transla- 
tions and three global rotations) has been applied 
between the data sets. 

step consists in determining the position of the 
SVT with respect to the DCH. Tracks with suf- 
ficient numbers of SVT and DCH hits are fit 
two times: once using only the DCH informa- 
tion and again using only the SVT hits. The six 
global alignment parameters, three translations 
and three rotations, are determined by minimiz- 
ing the difference between track parameters ob- 
tained with the SVT-only and the DCH-only fits. 
As reported above, because of the diurnal move- 
ment of the SVT with respect to the DCH, this 
global alignment needs to be performed once per 
run every 2-3 hours). The alignment con- 
stants obtained in a given run are then used to 
reconstruct the data in the subsequent run. This 
procedure, known as rolling calibration, ensures 
that track reconstruction is always performed 
with up-to-date global alignment constants. 
A record of the changes in the relative posi- 
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on front-end chips. Actually, since most of the 
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Figure 25. Global alignment of the SVT relative 
to the DCH based on e^e~ and /i events: 
changes in the relative vertical placement mea- 
sured a) over the entire ten-month run in the year 
2000, and b) a ten-day period, illustrating diurnal 
variations. 

tion of the SVT as determined by rolling calibra- 
tions is shown in Figure The position is stable 
to better than ±100 fim over several weeks, but 
changes abruptly from time to time, for instance, 
during access to the detector. The calibrations 
track diurnal variations of typically ±50 /im that 
have been correlated with local changes in tem- 
perature of about ±2°C. Movements within a 
single run are small compared to the size of the 
beam. 

5.6.3. Performance 

The SVT efficiency can be calculated for each 
half-module by comparing the number of asso- 
ciated hits to the number of tracks crossing the 
active area of the module. As can be seen in 
Figure a combined hardware and software ef- 
ficiencies of 97% is measured, excluding defective 
readout sections (9 out of 208), but employing 
no special treatment for other defects, such as 
broken AC coupling capacitors or dead channels 



Figure 26. SVT hit reconstruction efficiency, as 
measured on /i'^/.t" events for a) forward half- 
modules and b) backward half-modules. The 
plots show the probability of associating both a 
(f) and z hit to a track passing through the ac- 
tive part of the detector. The horizontal axis 
corresponds to the different modules, with the 
vertical lines separating the different layers as 
numbered. Missing values correspond to non- 
functioning half-modules. 

defects affect a single channel, they do not con- 
tribute to the inefficiency, because most tracks 
deposit charge in two or more strips due to track 
crossing angles, and charge diffusion. 

The spatial resolution of SVT hits is deter- 
mined by measuring the distance (in the plane of 
the sensor) between the track trajectory and the 
hit, using high- momentum tracks in two prong 
events. The uncertainty due to the track trajec- 
tory is subtracted from the width of the resid- 
ual distribution to obtain the hit resolution. Fig- 
ure |2^ shows the SVT hit resolution for z and (j> 
side hits as a function of track incident angle, for 
each of the five layers. The measured resolutions 
are in excellent agreement with expectations from 
Monte Carlo simulations. 

Initial studies have shown that hit reconstruc- 
tion efficiency and spatial resolution are effec- 
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Figure 27. SVT hit resolution in the a) z and b) 
(j) coordinate in microns, plotted as a function of 
track incident angle in degrees. Each plot shows 
a different layer of the SVT. The plots in the cj) 
coordinate for layers 1-3 are asymmetric around 
= because of the "pinwheel" design of the 
inner layers. There are fewer points in the </> res- 
olution plots for the outer layers as they subtend 
smaller angles than the inner layers. 



tively independent of occupancy for the occu- 
pancy levels observed so far. 

Measurement of the ToT value by the ATOM 
ICs enables one to obtain the pulse height, and 
hence the ionization dE/dx in the SVT sensor. 
The values of ToT are converted to pulse height 
using a lookup table computed from the pulse 
shapes obtained in the bench measurements. The 
pulse height is corrected for track length vari- 
ation. The double-sided sensors provide up to 
ten measurements of dE/dx per track. For every 
track with signals from at least four sensors in 
the SVT, a 60% truncated mean dE/dx is cal- 
culated. The cluster with the smallest dE/dx 
energy is also removed to reduce sensitivity to 
electronics noise. For MIPs, the resolution on the 
truncated mean dE/dx is approximately 14%. A 
2(7 separation between the kaons and pions can 
be achieved up to momentum of 500 MeV/ c, and 
between kaons and protons beyond 1 GeV/c. 

5.7. Summary and Outlook 

The SVT has been operating efficiently since 
its installation in the BABAR experiment in May 
1999. The five layer device, based on double- 
sided silicon sensors, has satisfied the original de- 
sign goals, in particular the targets specified for 
efficiency, hit resolution, and low transverse mo- 
mentum track reconstruction. The radiation dose 
during the first 25 fb~^ of integrated luminosity is 
within the planned budget, and no modules have 
failed due to radiation damage. The performance 
of the SVT modules at high radiation dose is cur- 
rently being studied. Early results indicate that 
the sensors will continue to function after type 
inversion (at 1-3 MRad), but further tests with 
irradiated sensors and ATOM ICs need to be per- 
formed. A program of spare module production 
has commenced, with the goal of replacing mod- 
ules that are expected to fail due to radiation 
damage. Beam-generated backgrounds are ex- 
pected to rise with increasing luminosity. Physics 
studies at five times the current backgrounds lev- 
els indicate no change in mass or vertex resolution 
for the mode B° J/ipK° and a ~ 20% loss of 
resolution in the D*+ — i?" mass difference. In 
this study the detector efficiency for both decay 
modes was lower by 15-20%. 
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6. Drift Chamber 

6.1. Purpose and Design Requirements 

The principal purpose of the drift chamber 
(DCH) is the efficient detection of charged parti- 
cles and the measurement of their momenta and 
angles with high precision. These high preci- 
sion measurements enable the reconstruction of 
exclusive B- and £)-meson decays with minimal 
background. The DCH complements the mea- 
surements of the impact parameter and the di- 
rections of charged tracks provided by the SVT 
near the IP. At lower momenta, the DCH mea- 
surements dominate the errors on the extrapola- 
tion of charged tracks to the DIRC, EMC, and 
IFR. 

The reconstruction of decay and interaction 
vertices outside of the SVT volume, for instance 
the K'^ decays, relies solely on the DCH. For this 
purpose, the chamber should be able to measure 
not only the transverse momenta and positions, 
but also the longitudinal position of tracks, with 
a resolution of ~1 mm. 

The DCH also needs to supply information for 
the charged particle trigger with a maximum time 
jitter of 0.5 ^is (Section ^Tj). 

For low momentum particles, the DCH is re- 
quired to provide particle identification by mea- 
surement of ionization loss (dE/dx). A resolu- 
tion of about 7% will allow t:/K separation up 
to 700MeV/c. This capability is complementary 
to that of the DIRC in the barrel region, while 
in the extreme backward and forward directions, 
the DCH is the only device providing some dis- 
crimination of particles of different mass. 

Since the average momentum of charged parti- 
cles produced in B- and Z?-meson decays is less 
than 1 GeV/c, multiple scattering is a significant, 
if not the dominant limitation on the track pa- 
rameter resolution. In order to reduce this contri- 
bution, material in front of and inside the cham- 
ber volume has to be minimized. 

Finally, the DCH must be operational in 
the presence of large beam-generated back- 
grounds, which were predicted to generate rates 
of ~5 kHz/cell in the innermost layers. 

6.2. Mechanical Design and Assembly 



6.2.1. Overview 

The DCH is relatively small in diameter, but al- 
most 3 m long, with 40 layers of small hexagonal 
cells providing up to 40 spatial and ionization loss 
measurements for charged particles with trans- 
verse momentum greater than 180MeV/c. Longi- 
tudinal position information is obtained by plac- 
ing the wires in 24 of the 40 layers at small angles 
with respect to the z-axis. By choosing low-mass 
aluminum field wires and a helium-based gas mix- 
ture, the multiple scattering inside the DCH is 
held to a minimum, less than 0.2%A'o of mate- 
rial. The properties of the chosen gas, a 80:20 
mixture of helium:isobutane, are presented in Ta- 
ble ||. This mixture has a radiation length that 
is five times larger than commonly used argon- 
based gases. The smaller Lorentz angle results in 
a rather uniform time-distance relationship and 
thereby improved spatial resolution. 

Table 8 

Properties of helium-isobutane gas mixture at at- 
mospheric pressure and 20°C. The drift velocity is 
given for operation without magnetic field, while 
the Lorentz angle is stated for a 1.5 T magnetic 
field. 



Parameter 


Values 


Mixture He : C4H10 


80:20 


Radiation Length 


807 m 


Primary Ions 


21.2/cm 


Drift Velocity 


22 /im/ ns 


Lorentz Angle 


32° 


dE/dx Resolution 


6.9% 



The inner cylindrical wall of the DCH is kept 
thin to facilitate the matching of the SVT and 
DCH tracks, to improve the track resolution for 
high momentum tracks, and to minimize the 
background from photon conversions and inter- 
actions. Material in the outer wall and in the 
forward direction is also minimized so as not to 
degrade the performance of the DIRC and the 
EMC. For this reason, the HV distribution and 
all of the readout electronics are mounted on the 
backward endplate of the chamber. This choice 
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Figure 28. Longitudinal section of the DCH with principal dimensions; the chamber center is offset by 
370mm from the interaction point (IP). 



also eliminates the need for a massive, heavily 
shielded cable plant. 

A longitudinal cross section and dimensions of 
the DCH are shown in Figure |8[ The DCH is 
bounded radially by the support tube at its in- 
ner radius and the DIRC at its outer radius. The 
device is asymmetrically located with respect to 
the IP. The forward length of 1749 mm is chosen 
so that particles emitted at polar angles of 17.2° 
traverse at least half of the layers of the chamber 
before exiting through the front endplate. In the 
backward direction, the length of 1015 mm means 
that particles with polar angles down to 152.6° 
traverse at least half of the layers. This choice en- 
sures sufficient coverage for forward-going tracks, 
and thus avoids significant degradation of the in- 
variant mass resolution, while at the same time 
maintaining a good safety margin on the electri- 
cal stability of the chamber. The DCH extends 
beyond the endplate by 485 mm at the backward 
end to accommodate the readout electronics, ca- 
bles, and an rf shield. It extends beyond the for- 
ward endplate by 68 mm to provide space for wire 
feed-throughs and an rf shield. 

6.2.2. Structural Components 

Details of the DCH mechanical design are pre- 
sented in Figure |2^. The endplates, which carry 
an axial load of 31,800 kN, are made from alu- 
minum plates of 24 mm thickness. At the forward 



end, this thickness is reduced to 12 mm beyond 
a radius of 46.9 cm to minimize the material in 
front of the calorimeter endcap. For this thick- 
ness, the estimated safety margin on the plastic 
yield point for endplate material (6061T651 alu- 
minum) is not more than a factor of two. The 
maximum total deflection of the endplates under 
loading is small, about 2 mm or 28% of the 7 mm 
wire elongation under tension. During installa- 
tion of the wires, this small deflection was taken 
into account by over-tensioning the wires. 

The inner and outer cylinder cylindrical walls 
are load bearing to reduce the maximum stress 
and deflections of the endplates. The stepped 
forward endplate created a complication during 
the assembly, because the thinner forward end- 
plate would deflect more than the thicker back- 
ward endplate. The outside rim of the forward 
endplate had to be pre-loaded, i.e., displaced by 
2.17 mm in the forward direction, to maintain the 
inside and outside rims of the rear endplate at the 
same longitudinal position after the load of the 
wires was transfered from the stringing fixture to 
the outer cylinder. 

Prior to installation on the inner cylinder, the 
two endplates were inspected on a coordinate- 
measuring machine. All sense wire holes, as well 
as 5% of the field and clearing field wire holes, 
were measured to determine their absolute loca- 
tions. The achieved accuracy of the hole place- 
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Honeycomb Carbon Fiber 




Figure 29. Details of the structural elements of the DCH. All components are made of aluminum, except 
for the 1 mm-thick inner beryllium wall and the 9 mm-thick outer composite wall. 



ment was 38 /im for both sense and field wires, 
better than the specification by more than a fac- 
tor of two. In addition, the diameters of the same 
sample of endplate holes were checked with pre- 
cision gauge pins. All holes passed the diame- 
ter specification (4.500±q;qoo for sense wires and 
2.500±q;qq^ for the field and guard wires). 

The inner cylindrical wall of the DCH, which 
carries 40% of the wire load, was made from 
five sections, a central 1 mm-thick beryllium tube 
with two aluminum extensions which were in 
turn electron-beam welded to two aluminum end 
flanges to form a 3m-long cylindrical part. The 
central section was made from three 120° seg- 
ments of rolled and brazed beryllium. The end 
flanges have precision surfaces onto which the 
endplates were mounted. These surfaces set the 
angles of the two endplates with respect to the 
axis and significantly constrain the concentricity 
of the tube. The inner cylinder also provides a 
substantial rf shield down to the PEP-11 bunch- 
gap frequency of 136 kHz. 

The outer wall bears 60% of axial wire load 
between the endplates. To simplify its installa- 
tion, this external wall was constructed from two 



half-cylinders with longitudinal and circumferen- 
tial joints. The gas and electrical seals for these 
joints were made up in situ. The main structural 
element consists of two 1.6 mm-thick (O.OO6X0) 
carbon fiber skins laminated to a 6 mm-thick hon- 
eycomb core. The outer shell is capable of with- 
standing a differential pressure of 30mbar and 
temperature variations as large as ±20°C, con- 
ditions that could be encountered during ship- 
ping or installation. Aluminum foil, 25 /im-thick 
on the inside surface and 100 /im on the outside, 
are in good electrical contact with the endplates, 
thereby completing the rf shield for the chamber. 

The total thickness of the DCH at normal in- 
cidence is 1.08%Aro, of which the wires and gas 
mixture contribute 0.2%Xo, and the inner wall 
0.28%Xo. 

6.2.3. Wire Feed-Throughs 

A total of five different types of feed-throughs 
were required for the chamber to accommodate 
the sense, field, and clearing field wires, as well as 
two different endplate thicknesses. The five types 
are illustrated in Figure |3^. They incorporate 
crimp pins |^3|] of a simple design which fasten 
and precisely locate the wires. The choice of pin 
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material (gold-plated copper for the signal wires 
and gold-plated aluminum for the field wires) and 
wall thickness in the crimp region was optimized 
to provide an allowable range of almost 150 /im in 
crimp size, as a primary means for avoiding wire 
breakage. 

Crimp pins were either press-fit into an in- 
sulator made from a single piece of injection- 
molded thermoplastic reinforced with 30% silica 
glass fiber |Q, or swaged into a copper jacket for 
the field wires. The plastic insulates the sense, 
guard, and clearing field wires from the electri- 
cally grounded endplates, while the metal jackets 
provide good ground contact for field wires (< 
0.1!7) on the backward endplate. The outer diam- 
eter of the field and clearing field feed-throughs 
was maintained at 2.000^q q25 while the sense 
wire feed-through had a larger (4.5001o;q25 mm) 
outer diameter and a longer body (41.7mm). 
This choice provided both thicker insulating walls 
and a longer projection into the gas volume to 
better shield the HV from the grounded endplate. 

6.2.4. Assembly and Stringing 

Assembly of the chamber components and in- 
stallation of the wires was carried out in a large 
clean room (Class 10,000) at TRIUMF in Vancou- 
ver. The wires were strung horizontally without 
the outer cylindrical shell in place. The endplates 
were mounted and aligned onto the inner cylinder 
which in turn was supported by a central shaft in 
a mobile fixture. The endplates were mounted on 
the inner cylinder at the inside rim and attached 
to support rings at the outside. These rings were 
connected by radial spiders to the central shaft 
of the stringing frame. 

Two teams of two operators each worked in 
parallel as the wires were strung from the in- 
ner radius outward. The two teams were each 
assisted by an automated wire transporter ]45| . 
A wire was attached to a needle which was in- 
serted through one of the endplate hole, captured 
magnetically by one of the transporters, and then 
transported and inserted though the appropriate 
hole in the other endplate. The wire was then 
threaded through the feed-throughs, which were 
glued into the endplates, and the wire was ten- 
sioned and crimped. The automated wire trans- 




1-^12 mrn—^ 



Figure 30. Design of the five DCH wire feed- 
throughs for the 24mm-thick endplates and the 
12mm-thick endplate. The copper jacketed feed- 
through is for grounded field wires, the other four 
are for sense wires (4.5 mm diameter), and guard 
and clearing field wires (2.5mm diameter), all 
made from a Celenex insulator surrounding the 
crimp pins. 

porters were largely built from industrial compo- 
nents, employing commercial software and hard- 
ware. The semi-automatic stringing procedure 
ensured the correct hole selection, accelerated the 
stringing rate and greatly improved the cleanli- 
ness and quality of the stringing process. The in- 
stallation of a total of 28,768 wires was completed 
in less than 15 weeks. 

6.3. Drift Cells 

6.3.1. Layer Arrangement 

The DCH consists of a total of 7,104 small drift 
cells, arranged in 40 cylindrical layers. The layers 
are grouped by four into ten superlayers, with the 
same wire orientation and equal numbers of cells 
in each layer of a superlayer. Sequential layers 
are staggered by half a cell. This arrangement 
enables local segment finding and left-right ambi- 
guity resolution within a superlayer, even if one 
out of four signals is missing. The stereo angles of 
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the superlayers alternate between axial (A) and 
stereo (U,V) pairs, in the order AUVAUVAUVA, 
as shown in Figure The stereo angles vary be- 
tween ±45mrad and ±76mrad; they have been 
chosen such that the drilling patterns are identi- 
cal for the two endplates. The hole pattern has a 
16-fold azimuthal symmetry which is well suited 
to the modularity of the electronic readout and 
trigger system. Table ^ summarizes parameters 
for all superlayers. 

Table 9 

The DCH superlayer (SL) structure, specifying 
the number of cells per layer, radius of the inner- 
most sense wire layer, the cell widths, and wire 
stereo angles, which vary over the four layers in 
a superlayer as indicated. The radii and widths 
are specified at the mid-length of the chamber. 





#of 


Radius 


Width 


Angle 


SL 


Cells 


(mm) 


(mm) 


(mrad) 


1 


96 


260.4 


17.0-19.4 





2 


112 


312.4 


17.5-19.5 


45-50 


3 


128 


363.4 


17.8-19.6 


-(52-57) 


4 


144 


422.7 


18.4-20.0 





5 


176 


476.6 


16.9-18.2 


56-60 


6 


192 


526.1 


17.2-18.3 


-(63-57) 


7 


208 


585.4 


17.7-18.8 





8 


224 


636.7 


17.8-18.8 


65-69 


9 


240 


688.0 


18.0-18.9 


-(72-76) 


10 


256 


747.2 


18.3-19.2 






6.3.2. Cell Design and Wires 

The drift cells are hexagonal in shape, 11.9 mm 
by approximately 19.0 mm along the radial and 
azimuthal directions, respectively. The hexago- 
nal cell configuration is desirable because approx- 
imate circular symmetry can be achieved over a 
large portion of the cell. The choice of aspect 
ratio has the benefit of decreasing the number 
of wires and electronic channels, while allowing a 
40-layer chamber in a confined radial space. Each 
cell consists of one sense wire surrounded by six 
field wires, as shown in Figure The proper- 
ties of the different types of gold-coated wires that 
make up the drift cells are given in Table |l^. The 
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Figure 31. Schematic layout of drift cells for 
the four innermost superlayers. Lines have been 
added between field wires to aid in visualization 
of the cell boundaries. The numbers on the right 
side give the stereo angles (mrad) of sense wires in 
each layer. The 1 mm-thick beryllium inner wall 
is shown inside of the first layer. 

sense wires are made of tungsten- rhenium 
20 ^m in diameter and tensioned with a weight 
of 30 g. The deflection due to gravity is 200 /im 
at mid-length. Tungsten-rhenium has a substan- 
tially higher linear resistivity (290^?/m), com- 
pared to pure tungsten (160i7/m), but it is con- 
siderably stronger and has better surface quality. 
While the field wires are at ground potential, a 
positive high voltage is applied to the sense wires. 
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An avalanche gain of approximately 5 x 10^ is 
obtained at a typical operating voltage of 1960 V 
and a 80:20 helium:isobutane gas mixture. 

Table 10 

DCH wire specifications (all wires are gold 
plated) . 



Diameter Voltage Tension 



Type 


Material 


(/xm) 


(V) 


(g) 


Sense 


W-Re 


20 


1960 


30 


Field 


Al 


120 





155 


Guard 


Al 


80 


340 


74 


Clearing 


Al 


120 


825 


155 
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Figure 32. Drift cell isochrones, i.e., contours of 
equal drift times of ions in cells of layers 3 and 4 of 
an axial superlayer. The isochrones are spaced by 
100 ns. They are circular near the sense wires, but 
become irregular near the field wires, and extend 
into the gap between superlayers. 

to the dE/dx measurement. 

6.3.4. Cross Talk 

A signal on one sense wire produces oppositely- 
charged signals on neighboring wires due to ca- 
pacitive coupling. The cross talk is largest be- 
tween adjacent cells of adjacent layers, ranging 
from —0.5% at a superlayer boundary to —2.7% 
for internal layers within superlayers. For adja- 
cent cells in the same layer, the cross talk ranges 
from —0.8 to —1.8%, while for cells separated by 
two layers it is less than 0.5%. 

6.4. Electronics 

6.4.1. Design Requirements and Overview 

The DCH electronic system is designed to pro- 
vide a measurement of the drift time and the inte- 
grated charge, as well as a single bit to the trigger 
system [Q for every wire with a signal. In the 
80:20 hclium:isobutane gas mixture, there are on 
average some 22 primary and 44 total ionization 
clusters produced per cm. The position of the 
primary ionization clusters is derived from tim- 
ing of the leading edge of the amplified signal. 
The design goal was to achieve a position res- 



The relatively low tension on the approxi- 
mately 2.75 m-long sense wires was chosen so that 
the aluminum field wires have matching gravita- 
tional sag and are tensioned well below the elas- 
tic limit. A simulation of the electrostatic forces 
shows that the cell configuration has no instabil- 
ity problems. At the nominal operating voltage 
of 1960 V, the wires deflect by less then 60 fim. 

The field wires |4^ are tensioned with 155 g 
to match the gravitational sag of the sense wires 
to within 20 /im. This tension is less than one- 
half the tensile yield strength of the aluminum 
wire. For cells at the inner or outer boundary of a 
superlayer, two guard wires are added to improve 
the electrostatic performance of the cell and to 
match the gain of the boundary cells to those of 
the cells in the inner layers. At the innermost 
boundary of layer 1 and the outermost boundary 
of layer 40, two clearing wires have been added 
per cell to collect charges created through photon 
conversions in the material of the walls. 

6.3.3. Drift Isochrones 

The calculated isochrones and drift paths for 
ions in adjacent cells of layer 3 and 4 of an ax- 
ial superlayer are presented in Figure ^ The 
isochrones are circular near the sense wires, but 
deviate greatly from circles near the field wires. 
Ions originating in the gap between superlayers 
are collected by cells in the edge layers after a de- 
lay of several /is. These lagging ions do not affect 
the drift times measurements, but they contribute 
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Figure 33. Block diagram for a 1/16*'' wedge of the DCH readout system, showing logical organization 
of the three front-end assembhes and their connections to the trigger and data I/O modules 



olution of 140 /ini, averaged over the cells. To 
reduce the time jitter in the signal arrival and 
at the same time maintain a good signal-to-noise 
ratio, the signal threshold was set at about 2.5 
primary electrons. For the dE/dx measurement, 
a resolution of 7% was projected for a 40-layer 
chamber. 

The small cell size and the difficult access 
through the DIRC strong support tube require a 
very high density of electronics components. As a 
consequence, a compact and highly modular de- 
sign was chosen. The readout is installed in well 
shielded assemblies that are plugged into the end- 
plate and are easily removable for maintenance. 

A schematic overview of the DCH electronics 
is presented in Figure 



33 



49[. The 16-fold az- 
imuthal symmetry of the cell pattern is reflected 
in the readout segmentation. The DCH amplifier 
and digitizer electronics are installed in electron- 
ics front-end assemblies (FEAs) that are mounted 
directly onto the rear endplate. There are three 
FEAs in each of the 16 sectors. These sectors 
are separated by brass cooling bars that extend 
from the inner to the outer chamber walls. These 
bars provide mechanical support and water cool- 



ing channel for the FEAs. The assemblies connect 
to the sense wires through service boards, which 
route the signals and HV distribution. A readout 
interface board (ROIB) in each FEA organizes 
the readout of the digitized data. Data I/O and 
trigger I/O modules multiplex serial data from 
the FEAs to high-speed optical fibers for transfer 
to the readout modules that are located in the 
electronics building. 

6.4.2. Service Boards 

Service boards provide the electrostatic poten- 
tials for signal, guard, and clearing wires, and 
pass signals and ground to the front-end read- 
out electronics. A side view of a service board is 



shown in Figure 34. The HV board contains the 



HV buses and filtering, current limiting resistors, 
and blocking capacitors. Jumpers connect adja- 
cent boards. The stored energy is minimized by 
using a 220 pF HV blocking capacitors. 

The signals are connected via series resistors to 
the upper signal board which contains the pro- 
tection diodes and standard output connectors. 
Mounting posts, anchored into the rear endplate, 
also serve as ground connections. 
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Figure 34. Side view of service boards show- 
ing two-tiered structure for DCH HV distribution 
and signal collection. 

6.4.3. Pront-End Assemblies 

The FEAs plug into connectors on the back 
side of the service boards. These custom wedge- 
shaped crates are aluminum boxes that contain a 
ROIB and two, three, or four amplifier/digitizer 
boards (ADB) for superlayers 1-4, 5-7, and 8-10, 
respectively, as shown in Figure 35. The crates 



are mounted with good thermal contact to the 
water cooled radial support bars. The total heat 
load generated by the FEAs is 1.3 kW. 

The ADBs are built from basic building blocks 
consisting of two 4-channel amplifier ICs ||5^ feed- 
ing a single 8-channel digitizer custom ASIC jsi] . 
The number of channels serviced by an ADB is 
60, 48, or 45, for the inner, middle, and outer 
FEA modules, respectively. 

The custom amplifier IC receives the input sig- 
nal from the sense wire and produces a discrim- 
inator output signal for the drift time measure- 
ment and a shaped analog signal for the dE/dx 
measurement. Both outputs are fully differential. 
The discriminator has gain and bandwidth con- 
trol, and a voltage controlled threshold. The ana- 
log circuit has integrator and gain control. 

The custom digitizer IC incorporates a 4-bit 
TDC for time measurement and a 6-bit 15 MHz 
FADC to measure the total deposited charge. 
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Figure 35. Layout of 1/16*'' of the DCH rear end- 
plate, showing three FEA boxes between water 
cooled channels. 

The TDC is a phase-locked digital delay linear 
vernier on the sample clock of 15 MHz, which 
achieves a 1 ns precision for leading edge tim- 
ing. The FADC design is based on a resistor- 
divider comparator ladder that operates in bi- 
linear mode to cover the full dynamic range. The 
digitized output signals are stored in a trigger la- 
tency buffer for 12.9 /iS, after which a LI Accept 
initiates the transfer of a 2.2 /is block of data to 
the readout buffer. In addition, trigger informa- 
tion is supplied for every channel, based either 
on the presence of a TDC hit during the sam- 
ple period or FADC differential pulse height in- 
formation, should a higher discriminator level be 
desirable. 

The ROIB interprets FCTS commands to con- 
trol the flow of data and trigger information. 
Data are moved to FIFOs on the ROIBs, and then 
to data and trigger I/O modules via 59.5 MHz se- 
rial links. A total of four such links are required 
per 1/16*'' wedge, one for each of the outer two 
FEAs and two for the innermost of the FEA. Each 
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data I/O module services all FEAs one quad- 
rant and transmits the data to a single ROM 
via one optical fiber link. The trigger stream is 
first multiplexed onto a total of 30 serial lines per 
wedge for transmission to the trigger I/O mod- 
ule. Trigger data from two wedges of FEAs are 
then transmitted to the trigger system via three 
optical links. Thus, a total of 28 optical fibers, 
four for the data and 24 for the charged particle 
trigger, are required to transfer the DCH data to 
the readout. 

6.4.4. Data Acquisition 

The data stream is received and controlled by 
four BABAR standard readout modules. Drift 
chamber-specific feature extraction algorithms 
convert the raw FADC and TDC information into 
drift times, total charge, and a status word. The 
time and charge are corrected channel-by-channel 
for time offsets, pedestals, and gain constants. 
Based on measurements of the noise a thresh- 
old is typically 2-3 electrons is applied to dis- 
criminate signals. These algorithms take about 
1 us per channel, and reduce the data volume by 
roughly a factor of four. 

6.4.5. High Voltage System 

The HV bias lines on the chamber are daisy- 
chained together so that each superlayer requires 
only four power supplies, except for superlayer 
1 which has eight. The voltages are supplied to 
the sense, guard, and clearing wires by a CAEN 
SY527 HV mainframe Q, equipped with 24- 
channel plug-in modules. The sense wires are sup- 
plied by 44 HV channels providing up to 40 /zA of 
current each that can be monitored with a reso- 
lution of 0.1%. 

6.5. Gas System 

The gas system has been designed to provide a 
stable 80:20 heliumiisobutane mixture at a con- 
stant over pressure of 4mbar The chamber 
volume is about 5.2 m^. Gas mixing and recir- 
culation is controlled by precise mass flow con- 
trollers; the total flow is tuned to 15i!/min, of 
which 2.5 ^/min are fresh gas. During normal 
operation, the complete DCH gas volume is re- 
circulated in six hours, and one full volume of 
fresh gas is added every 36 hours. The pressure 



in the DCH is measured by two independent pres- 
sure gauges, one of which is connected to a regula- 
tor controlling the speed of the compressor. The 
relative pressure in the chamber is controlled to 
better than ±0.05 mbar. 

Oxygen is removed from the gas mixture us- 
ing a palladium catalytic filter. The water con- 
tent is maintained at 3500 ± 200 ppm by passing 
an adjustable fraction of the gas through a water 
bubbler. This relatively high level of water vapor 
is maintained to prevent electrical discharge. In 
addition to various sensors to monitor pressure, 
temperature, and flow at several points of the sys- 
tem, a small wire chamber with an ^^Fe source 
continuously monitors gain of the gas mixture. 

6.6. Calibrations and Monitoring 

6.7. Electronics Calibration 

The front-end electronics (FEEs) are calibrated 
daily to determine the channel-by-channel correc- 
tion constants and thresholds. Calibration pulses 
are produced internally and input to the pream- 
plifier at a rate of about 160 Hz. The calibration 
signals are processed in the ROM to minimize 
the data transfer and fully exploit the available 
processing power. The results are stored for sub- 
sequent feature extraction. The entire online cal- 
ibration procedure takes less than two minutes. 

6.7.1. Environmental Monitoring 

The operating conditions of the DCH are mon- 
itored in realtime by a variety of sensors and read 
out by the detector-wide CAN bus system. These 
sensors monitor the flow rate, pressure, and gas 
mixture; the voltages and currents applied to the 
wires in the chamber; the voltages and currents 
distributed to the electronics from power sup- 
plies and regulators; instantaneous and cumula- 
tive radiation doses; temperature and humidity 
around the chamber electronics and in the equip- 
ment racks. Additional sensors monitor the at- 
mosphere in and around the detector for excess 
isobutane, which could pose a flammability or ex- 
plosive hazard in the event of a leak. 

Many of the sensors are connected to hardware 
interlocks, which ensure that the chamber is au- 
tomatically put into a safe state in response to an 
unsafe condition. All of these systems have per- 
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formed reliably. In addition, automated software 
monitors raw data quality, chamber occupancies 
and efficiencies to sense variations in electronics 
performance that might indicate more subtle op- 
erational problems. 

6.7.2. Operational Experience 

The design of the DCH specifies a voltage of 
1960 V on the sense wires to achieve the desired 
gain and resolution. The chamber voltage was 
lowered for part of the run to 1900 V out of con- 
cern for a small region of the chamber that was 
damaged during the commissioning phase by in- 
advertently applying 2kV to the guard wires. 
Wires in this region (10.4% of superlayer 5, and 
4.2% of superlayer 6) were disconnected when 
continuous discharge was observed over extended 
periods of time. 

6.8. Performance 

The DCH was first operated with full mag- 
netic field immediately after the installation into 
BABAR. Cosmic ray data were recorded and ex- 
tensive studies of the basic cell performance were 
performed to develop calibration algorithms for 
the time-to-distance and dE/dx measurements. 
These algorithms were then implemented as de- 
scribed below for colliding beam data. Calibra- 
tions are monitored continuously to provide feed- 
back to the operation; some time varying parame- 
ters are updated continuously as part of OPR. For 
charge particle tracking the DCH and SVT infor- 
mation is combined; the performance of the com- 
bined tracking system is described in Section 0. 

6.8.1. Time-to-Distance Relation 

The precise relation between the measured drift 
time and drift distance is determined from sam- 
ples of e^e^ and /i^/i^ events. For each signal, 
the drift distance is estimated by computing the 
distance of closest approach between the track 
and the wire. To avoid bias, the fit does not use 
the hit on the wire under consideration. The es- 
timated drift distances and measured drift times 
are averaged over all wires in a layer, but the data 
are accumulated separately for tracks passing on 
the left of a sense wire and on the right. The time- 
distance relation is fit to a sixth-order Chebychev 
polynomial. An example of such a fit is shown in 
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Figure 36. The drift time versus distance rela- 
tion for left and right half of a cell. These func- 
tions are obtained from the data averaged over all 
cells in a single layer of the DCH. 

Figure 

An additional correction is made for tracks 
with varying entrance angle into the drift cell. 
This angle is defined relative to the radial vector 
from the IP to the sense wire. The correction is 
applied as a scale factor to the drift distance and 
was determined layer-by-layer from a Garfield js^ 
simulation. The entrance angle correction is im- 
plemented as a tenth-order Chebychev polyno- 
mial of the drift distance, with coefficients which 
are functions of the entrance angle. 

Figure |3^ shows the position resolution as a 
function of the drift distance, separately for the 
left and the right side of the sense wire. The 
resolution is taken from Gaussian fits to the dis- 
tributions of residuals obtained from unbiased 
track fits. The results are based on multi-hadron 
events, for data averaged over all cells in layer 18. 

6.8.2. Charge Measurement 

The specific energy loss, dE/dx^ for charged 
particles traversing the DCH is derived from mea- 
surement of total charge deposited in each drift 
cell. The charge collected per signal cell is mea- 
sured as part of the feature extraction algorithm 
in the ROM. Individual measurements are cor- 
rected for gain variations, pedestal-subtracted 
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Figure 37. DCH position resolution as a function 
of the drift distance in layer 18, for tracks on the 
left and right side of the sense wire. The data are 
averaged over all cells in the layer. 



and integrated over a period of approximately 

1.8^*8. 

The specific energy loss per track is computed 
as a truncated mean from the lowest 80% of the 
individual dE/dx measurements. Various correc- 
tions are applied to remove sources of bias that 
degrade the accuracy of the primary ionization 
measurement. These corrections include the fol- 
lowing: 

• changes in gas pressure and temperature, 
leading to ±9% variation in dE/dx, cor- 
rected by a single overall multiplicative con- 
stant; 

• differences in cell geometry and charge col- 
lection (±8% variation), corrected by a set 
of multiplicative constants for each wire; 

• signal saturation due to space charge build- 
up (±11% variation), corrected by a second- 
order polynomial in the dip angle. A, of the 
form 1/v sin^ A + const: 
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Figure 38. Measurement of dE/dx in the DCH 
as a function of track momenta. The data include 
large samples of beam background triggers, as ev- 
ident from the high rate of protons. The curves 
show the Bethe-Bloch predictions derived from 
selected control samples of particles of different 
masses. 



• non-linearities in the most probable energy 
loss at large dip angles (±2.5% variation), 
corrected with a fourth-order Chebychev 
polynomial as a function of A; and 

• variation of cell charge collection as a func- 
tion entrance angle (±2.5% variation), cor- 
rected using a sixth-order Chebychev poly- 
nomial in the entrance angle. 

The overall gas gain is updated continuously 
based on calibrations derived as part of prompt 
reconstruction of the colliding beam data; the 
remaining corrections are determined once for a 
given HV voltage setting and gas mixture. 

Corrections applied at the single-cell level can 
be large compared to the single-cell dE/dx res- 
olution, but have only a modest impact on the 
average resolution of the ensemble of hits. Global 
corrections applied to all hits on a track arc there- 
fore the most important for the resolution. 

Figure ^ shows the distribution of the cor- 
rected dE/dx measurements as a function of track 
momenta. The superimposed Bethe-Bloch pre- 
dictions for particles of different masses have been 
determined from selected control samples. 
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Figure 39. Difference between the measured and 
expected energy loss dE/dx for from Bhabha 
scattering, measured in the DCH at an operating 
voltage of 1900 V. The curve represents a Gaus- 
sian fit to the data with a resolution of 7.5%. 

The measured dE/dx resolution for Bhabha 
events is shown in Figure |3^. The rms resolu- 
tion achieved to date is typically 7.5%, limited 
by the number of samples and Landau fluctua- 
tions. This value is close to the expected resolu- 
tion of 7%. Further refinements and additional 
corrections are being considered to improve per- 
formance. 

6.9. Conclusions 

The DCH has been performing close to design 
expectations from the start of operations. With 
the exception of a small number of wires that were 
damaged by an unfortunate HV incident during 
the commissioning phase, all cells are fully oper- 
ational. The DCH performance has proven very 
stable over time. The design goal for the intrin- 
sic position and dE/dx resolution have been met. 
Backgrounds are acceptable at present beam cur- 
rents, but there is concern for rising occupancies 
and data acquisition capacity at the high end of 
the planned luminosity upgrades. 

7. Performance of the Charged Particle 
Tracking Systems 

Charged particle tracking has been studied 
with large samples of cosmic ray muons, e"'"e~. 



/i^/i^, and T^r" events, as well as multi-hadrons. 
At this time, these studies are far from complete 
and the results represent the current status. In 
particular, many issues related to the intrinsic 
alignment of the SVT and the DCH, the varia- 
tion with time of the relative alignment of the 
SVT and the DCH, and movement of the beam 
position relative to BABAR remain under study. 

7.1. Track Reconstruction 

The reconstruction of charged particle tracks 
relies on data from both tracking systems, the 
SVT and the DCH. Charged tracks are defined 
by five parameters (c?o, '/'Oj ^Oi tan A) and their 
associated error matrix. These parameters are 
measured at the point of closest approach to the 
z-axis; c?o find zq are the distances of this point 
from the origin of the coordinate system in the 
x-y plane and along the z-axis, respectively. The 
angle 0o is the azimuth of the track, A the dip an- 
gle relative to the transverse plane, and w = 1/pt 
is its curvature, do and lu are signed variables; 
their sign depends on the charge of the track. The 
track finding and the fitting procedures make use 
of Kalman filter algorithm that takes into ac- 
count the detailed distribution of material in the 
detector and the full map of the magnetic field. 

The offline charged particle track reconstruc- 
tion builds on information available from the L3 
trigger and tracking algorithm. It begins with an 
improvement of the event start time to, obtained 
from a fit to the parameters do, 0Oi Sind to based 
on the four-hit track segments in the DCH su- 
pcrlayers. Next, tracks are selected by perform- 
ing helix fits to the hits found by the L3 track 
finding algorithm. A search for additional hits in 
the DCH that may belong on these tracks is per- 
formed, while to is further improved by using only 
hits associated with tracks. Two more sophisti- 
cated tracking procedures are applied which are 
designed to find tracks that either do not pass 
through the entire DCH or do not originate from 
the IP. These algorithms primarily use track seg- 
ments that have not already been assigned to 
other tracks, and thus benefit from a progressively 
cleaner tracking environment with a constantly 
improving Iq. At the end of this process, tracks 
are again fit using a Kalman filter method. 
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Figure 40. The track reconstruction efficiency 
in the DCH at operating vohages of 1900 V and 
1960 V, as a function of a) transverse momentum, 
and b) polar angle. The efficiency is measured in 
multi-hadron events as the fraction of all tracks 
detected in the SVT for which the DCH portion 
is also reconstructed. 

The resulting tracks are then extrapolated into 
the SVT, and SVT track segments are added, pro- 
vided they are consistent with the expected error 
in the extrapolation through the intervening ma- 
terial and inhomogcncous magnetic field. Among 
the possible SVT segments, those with the small- 
est residuals and the largest number of SVT layers 
are retained and a Kalman fit is performed to the 
fuU set of DCH and SVT hits. 

Any remaining SVT hits are then passed to 
two complementary standalone track finding al- 
gorithms. The first reconstructs tracks starting 
with triplets of space points (matched (j) and z 
hits) in layers 1,3, and 5 of the SVT, and adding 
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Figure 41. Reconstruction of low momen- 
tum tracks: the mass difference, AA/ = 
M(K--K+'K+) - M{K-Ti+), both for aU detected 
events (data points) and for events in which the 
low momentum pion is reconstructed both in the 
SVT and DCH (histogram). Backgrounds from 
combinatorics and fake tracks, as well as non- 
resonant data have been subtracted. 

consistent space points from the other layers. A 
minimum of four space points are required to form 
a good track. This algorithm is efficient over a 
wide range of dp ^-nd -^o values. The second al- 
gorithm starts with circle trajectories from cj) hits 
and then adds z hits to form helices. This al- 
gorithm is less sensitive to large combinatorics 
and to missing z information for some of the SVT 
modules. 

Finally, an attempt is made to combine tracks 
that are only found by one of the two tracking 
systems and thus recover tracks scattered in the 
material of the support tube. 

7.2. Tracking EfRciency 

The efficiency for reconstructing tracks in the 
DCH has been measured as a function of trans- 
verse momentum, polar and azimuthal angles in 
multi-track events. These measurements rely on 
specific final states and exploit the fact that the 
track reconstruction can be performed indepen- 
dently in the SVT and the DCH. 
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Figure 42. Monte Carlo studies of low momentum 
tracks in the SVT: a) comparison of data (contri- 
butions from combinatoric background and non- 
SB events have been subtracted) with simulation 
of the transverse momentum spectrum of pions 
from Z)*+ D°7r+ in BB events, and b) effi- 
ciency for slow pion detection derived from simu- 
lated events. 

The absolute DCH tracking efficiency is deter- 
mined as the ratio of the number of reconstructed 
DCH tracks to the number of tracks detected 
in the SVT, with the requirement that they fall 
within the acceptance of the DCH. Such stud- 
ies have been performed for different samples of 
multi-hadron events. Figure ^ shows the re- 
sult of one such study for the two voltage set- 
tings. The measurement errors are dominated by 
the uncertainty in the correction for fake tracks 
in the SVT. At the design voltage of 1960V, 
the efficiency averages 98 ± 1% per track above 
200MeV/c and polar angle > 500mrad. The 
data recorded at 1900 V show a reduction in effi- 
ciency by about 5% for tracks at close to normal 
incidence, indicating that the cells are not fully 
efficient at this voltage. 

The standalone SVT tracking algorithms have 
a high efficiency for tracks with low transverse 
momentum. This feature is important for the de- 
tection oi D* decays. To study this efficiency, 
decays _D*+ D^n^ are selected by recon- 
structing events of the type B D*^X fol- 
lowed by D*+ D°TT+ K-TT+Tr+. For the 



majority of these low momentum pions the mo- 
mentum resolution is limited by multiple scat- 
tering, but the production angle can be deter- 
mined from the signals in innermost layers of 
the SVT. Figure O shows the mass difference 
AAf = M{K-Tr+7^) - M[K-TT+), for the to- 
tal sample and the subsample of events in which 
the slow pion has been reconstructed in both the 
SVT and the DCH. The difference in these two 
distributions demonstrates the contribution from 
SVT standalone tracking, both in terms of the 
gain of signal events and of resolution. The gain 
in efficiency is mostly at very low momenta, and 
the resolution is impacted by multiple scattering 
and limited track length of the slow pions. To 
derive an estimate of the tracking efficiency for 
these low momentum tracks, a detailed Monte 
Carlo simulation was performed. Specifically, the 
pion spectrum was derived from simulation of the 
inclusive D* production in BB events, and the 
Monte Carlo events were selected in the same way 
as the data. A comparison of the detected slow 
pion spectrum with the Monte Carlo prediction is 
presented in Figure ^ Based on this very good 
agreement, the detection efficiency has been de- 
rived from the Monte Carlo simulation. The SVT 
significantly extends the capability of the charged 
particle detection down to transverse momenta of 
~50MeV/c. 

7.3. Track Parameter Resolutions 

The resolution in the five track parameters is 
monitored in OPR using e^e^ and iJ-~^fi~ pair 
events. It is further investigated offline for tracks 
in multi-hadron events and cosmic ray muons. 

Cosmic rays that are recorded during normal 
data-taking offer a simple way of studying the 
track parameter resolution. The upper and lower 
halves of the cosmic ray tracks traversing the 
DCH and the SVT are fit as two separate tracks, 
and the resolution is derived from the difference of 
the measured parameters for the two track halves. 
To assure that the tracks pass close to the beam 
interaction point, cuts are applied on the dg, zg, 
and tan A. The results of this comparison for the 
coordinates of the point of closest approach and 
the angles are shown in Figure ^ for tracks with 
momenta above pt of 3 GeV/c. The distributions 
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Figure 43. Measurements of the differences between the fitted track parameters of the two halves of 
cosmic ray muons, with transverse momenta above 3GeV/c, a) Ado, b) Azq, c) A(/)o, and d) AtanA. 



are symmetric; the non-Gaussian tails are small. 
The distributions for the differences in zo and 
tan A show a clear offset, attributed to residual 
problems with the internal alignment of the SVT. 
Based on the full width at half maximum of these 
distributions the resolutions for single tracks can 
be parametized as 

ad„ = 23//m (t^,-, ~ 0.43 mrad 
a,,=29^lm dtan a = 0.53 • 10~3. 

The dependence of the resolution in c?o a-nd zq 
on the transverse momentum pt is presented in 
Figure The measurement is based on tracks 
in multi-hadron events. The resolution is deter- 
mined from the width of the distribution of the 
difference between the measured parameters, do 
and Zq, and the coordinates of the vertex recon- 
structed from the remaining tracks in the event. 
These distributions peak at zero, but have a tail 
for positive values due to the effect of particle de- 
cays. Consequently, only the negative part of the 
distributions reflects the measurement error and 
is used to extract the resolution. Event shape 
cuts and a cut on the of the vertex fit are ap- 
plied to reduce the effect of weak decays on this 
measurement. The contribution from the vertex 
errors are removed from the measured resolutions 
in quadrature. The do and zq resolutions so mea- 
sured are about 25 /im and 40 /im respectively at 
Pt = 3 GeV/c. These values agree well with ex- 
pectations, and are also in reasonable agreement 
with the results obtained from cosmic rays. 
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Figure 44. Resolution in the parameters do and 
Zq for tracks in multi-hadron events as a function 
of the transverse momentum. The data are cor- 
rected for the effects of particle decays and ver- 
texing errors. 

Figure ^ shows the estimated error in the mea- 
surement of the difference along the z-axis be- 
tween the vertices of the two neutral B mesons, 
one of them is fully reconstructed, the other 
serves as a flavor tag. The rms width of 190 /im 
is dominated by the reconstruction of the tagging 
B vertex, the rms resolution for the fully recon- 
structed B meson is 70 fim. The data meet the 
design expectation 

While the position and angle measurements 
near the IP are dominated by the SVT measure- 
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Figure 45. Distribution of the error on the dif- 
ference Az between the B meson vertices for a 
sample of events in which one B'^ is fuUy recon- 
structed. 



2.0 



1.0 - 



1-2001 
8583A23 




4 8 

Transverse Momentum (GeV/c) 



Figure 46. Resolution in the transverse mo- 
mentum pt determined from cosmic ray muons 
traversing the DCH and SVT. 



ments, the DCH contributes primarily to the pt 
measurement. Figure ^ shows the resolution in 
the transverse momentum derived from cosmic 
muons. The data are well represented by a linear 
function, 

apjpt = (0.13 ± 0.01)% • Pt + (0.45 ± 0.03)%, 

where the transverse momentum pt is measured 
in GeV/ c. These values for the resolution param- 
eters are very close to the initial estimates and 
can be reproduced by Monte Carlo simulations. 
More sophisticated treatment of the DCH time- 
to-distance relations and overall resolution func- 
tion are presently under study. 

Figure ^ shows the mass resolution for J/ip 
mesons reconstructed in the final state, 

averaged over all data currently available. The 
reconstructed peak is centered 0.05% below the 
expected value, this difference is attributed to 
the remaining inaccuracies in the SVT and DCH 
alignment and in the magnetic field parameteriza- 
tion. The observed mass resolution differs by 15% 
for data recorded at the two DCH HV settings, 
it is 13.0 ±0.3 MeV/c^ and 11.4 ± 0.3 MeV/c^ at 
1900 V and 1960 V, respectively. 
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Figure 47. Reconstruction of the decay J/tp 
in selected BB events. 

7.4. Summary 

The two tracking devices, the SVT and DCH, 
have been performing close to design expectations 
from the start of operations. Studies of track res- 
olution at lower momenta and as a function of 
polar and azimuthal angles are still under way. 
Likewise, the position and angular resolution at 
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the entrance to the DIRC or EMC are stiU being 
studied. Such measurements are very sensitive to 
internal ahgnment of the SVT and relative place- 
ment of the SVT and the DCH. A better under- 
standing will not only reduce the mass resolution 
for the reconstruction of exclusive states, it will 
also be important for improvement of the perfor- 
mance of the DIRC. 

8. DIRC 

8.1. Purpose and Design Requirements 

The study of CP-violation requires the abil- 
ity to tag the flavor of one of the B mesons via 
the cascade decay 6 — > c — > s, while fully re- 
constructing the second B decay. The momenta 
of the kaons used for flavor tagging extend up 
to about 2 GeV/c, with most of them below 1 
GeV/c. On the other hand, pions and kaons 
from the rare two-body decays — > tt+tt^ and 
— > K'^Tr~ must be well-separated. They 
have momenta between 1.7 and 4.2 GeV/c with a 
strong momentum-polar angle correlation of the 
tracks (higher momenta occur at more forward 
angles because of the cm. system boost) 

The Particle Identification (PID) system 
should be thin and uniform in terms of radiation 
lengths (to minimize degradation of the calorime- 
ter energy resolution) and small in the radial di- 
mension to reduce the volume, and hence, the 
cost of the calorimeter. Finally, for operation at 
high luminosity, the PID system needs fast sig- 
nal response, and should be able to tolerate high 
backgrounds. 

The PID system being used in BABAR is a new 
kind of ring-imaging Chercnkov detector called 
the DIRC 1^ (the acronym DIRC stands for De- 
tector of Internally Reflected Cherenkov light). It 
is expected to be able to provide n/K separation 
of 4 cr or greater, for all tracks from S-meson 
decays from the pion Cherenkov threshold up to 
4.2 GeV/c. PID below 700 MeV/c rehes primar- 
ily on the dE/dx measurements in the DCH and 
SVT. 

8.2. DIRC Concept 

The DIRC is based on the principle that the 
magnitudes of angles are maintained upon re- 



flection from a flat surface. Figure ^ shows a 
schematic of the DIRC geometry that illustrates 
the principles of light production, transport, and 
imaging. The radiator material of the DIRC is 
synthetic, fused silica in the form of long, thin 
bars with rectangular cross section. These bars 
serve both as radiators and as light pipes for the 
portion of the light trapped in the radiator by 
total internal reflection. Fused, synthetic silica 
(Spectrosil ||5^) is chosen because of its resistance 
to ionizing radiation, its long attenuation length, 
large index of refraction, low chromatic dispersion 
within the wavelength acceptance of the DIRC, 
and because it allows an excellent optical finish 
on the surfaces of the bars p8|| . 

In the following, the variable 9c is used to des- 
ignate the Cherenkov angle, (pc denotes the az- 
imuthal angle of a Cherenkov photon around the 
track direction, and n represents the mean index 
of refraction of fused silica {n — 1.473), with the 
familiar relation cos 0c — {(3 — v/c, v — 

velocity of the particle, c — velocity of light). 

For particles with (3 ~ I, some photons will al- 
ways lie within the total internal reflection limit, 
and will be transported to either one or both ends 
of the bar, depending on the particle incident an- 
gle. To avoid instrumenting both ends of the bar 
with photon detectors, a mirror is placed at the 
forward end, perpendicular to the bar axis, to 
reflect incident photons to the backward, instru- 
mented end. 

Once photons arrive at the instrumented end, 
most of them emerge into a water-filled expan- 
sion region, called the standoff box. A fused silica 
wedge at the exit of the bar reflects photons at 
large angles relative to the bar axis. It thereby 
reduces the size of the required detection sur- 
face and recovers those photons that would oth- 
erwise be lost due to internal reflection at the 
fused silica/ water interface. The photons are 
detected by an array of densely packed photo- 
multiplier tubes (PMTs), each surrounded by re- 
flecting light catcher cones ]59[ | to capture light 
which would otherwise miss the active area of 
the PMT. The PMTs are placed at a distance 
of about 1.2 m from the bar end. The expected 
Cherenkov light pattern at this surface is essen- 
tially a conic section, where the cone opening- 
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Figure 48. Schematics of the DIRC fused sihca 
radiator bar and imaging region. Not shown is a 
6 mrad angle on the bottom surface of the wedge 
(see text). 

angle is the Cherenkov production angle modi- 
fied by refraction at the exit from the fused silica 
window. 

The DIRC is intrinsically a three-dimensional 
imaging device, using the position and arrival 
time of the PMT signals. Photons generated in 
a bar are focused onto the phototube detection 
surface via a "pinhole" defined by the exit aper- 
ture of the bar. In order to associate the photon 
signals with a track traversing a bar, the vector 
pointing from the center of the bar end to the 
center of each PMT is taken as a measure of the 
photon propagation angles ax,ay, and az ■ Since 
the track position and angles are known from the 
tracking system, the three a angles can be used 
to determine the two Cherenkov angles 9c and 0c. 
In addition, the arrival time of the signal provides 
an independent measurement of the propagation 
of the photon, and can be related to the propaga- 
tion angles a. This over-constraint on the angles 
and the signal timing are particularly useful in 
dealing with ambiguities in the signal association 
(see Section 8.6.1) and high background rates. 



8.3. Mechanical Design 

and Physical Description 

The DIRC bars are arranged in a 12-sided 
polygonal barrel. Because of the beam energy 



asymmetry, particles are produced preferentially 
forward in the detector. To minimize interfer- 
ence with other detector systems in the forward 
region, the DIRC photon detector is placed at the 
backward end. 

The principal components of the DIRC are 
shown schematically in Figure The bars 

are placed into 12 hermetically sealed containers, 
called bar boxes, made of very thin aluminum- 
hexcel panels. Each bar box, shown in Figure |50| , 
contains 12 bars, for a total of 144 bars. Within 
a bar box the 12 bars are optically isolated by a 

IbOfim air gap between neighboring bars, en- 
forced by custom shims made from aluminum foil. 

The bars are 17 mm-thick, 35 mm- wide, and 4.9 
m-long. Each bar is assembled from four 1.225 m 
pieces that are glued end-to-end; this length is the 
longest high-quality bar currently obtainable js^ , 

The bars are supported at 600 mm intervals by 
small nylon buttons for optical isolation from the 
bar box. Each bar has a fused silica wedge glued 
to it at the readout end. The wedge, made of the 
same material as the bar, is 91 mm-long with very 
nearly the same width as the bars (33 mm) and a 
trapezoidal profile (27 mm-high at bar end, and 
79 mm at the light exit end) . The bottom of the 
wedge (see Figure Q has a slight (~ 6 mrad) up- 
ward slope to minimize the displacement of the 
downward reflected image due to the finite bar 
thickness. The twelve wedges in a bar box are 
glued to a common 10 mm-thick fused silica win- 
dow, that provides the interface and seal to the 
purified water in the standoff box. 

The mechanical support of the DIRC, shown 
in Figure |4^, is cantilevered from the steel of 
the IFR. The Strong Support Tube (SST) is a 
steel cylinder located inside the end doors of the 
IFR and provides the basic support for the entire 
DIRC. In turn, the SST is supported by a steel 
support gusset that fixes it to the barrel magnet 
steel. It also minimizes the magnetic flux gap 
caused by the DIRC bars extending through the 
flux return, and supports the axial load of the in- 
ner magnetic plug surrounding the beam in this 
region. 

The bar boxes are supported in the active re- 
gion by an aluminum tube, the Central Support 
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Figure 49. Exploded view of the DIRC mechani- 
cal support structure. The steel magnetic shield 
is not shown. 
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Figure 52. Transverse section of the nominal 
DIRC bar box imbedded in the CST. All dimen- 
sions are given in mm. 




Figure 50. Schematics of the DIRC bar box as- 
sembly. 



Tube (CST), attached to the SST via an alu- 
minum transition flange. The CST is a thin, 
double- walled, cylindrical shell, using aircraft- 
type construction with stressed aluminum skins 
and bulkheads having riveted or glued joints. The 
CST also provides the support for the DCH. 

The standoff box is made of stainless steel, 
consisting of a cone, cylinder, and 12 sectors of 
PMTs. It contains about 6,000 liters of purified 
water. Water is used to fill this region because 
it is inexpensive and has an average index of re- 
fraction {n « 1.346) reasonably close to that of 



fused silica, thus minimizing the total internal re- 
flection at silica-water interface. Furthermore, its 
chromaticity index is a close match to that of 
fused silica, effectively eliminating dispersion at 
the silica-water interface. The steel gusset sup- 
ports the standoff box. A steel shield, supple- 
mented by a bucking coil, surrounds the standoff 
box to reduce the field in the PMT region to be- 
low 1 Gauss 

The PMTs at the rear of the standoff box lie 
on a surface that is approximately toroidal. Each 
of the 12 PMT sectors contains 896 PMTs (ETL 
model 9125 ||6^,^) with 29 mm-diameter, in a 
closely packed array inside the water volume. A 
double o-ring water seal is made between the 
PMTs and the wall of the standoff box. The 
PMTs are installed from the inside of the standoff 
box and connected via a feed-through to a base 
mounted outside. The hexagonal light catcher 
cone is mounted in front of the photocathode of 
each PMT, which results in an effective active sur- 
face area light collection fraction of about 90%. 
The geometry of the DIRC is shown in Figures |5^ 
and 

The DIRC occupies 80 mm of radial space in 
the central detector volume including supports 
and construction tolerances, with a total of about 
17% radiation length thickness at normal inci- 
dence. The radiator bars subtend a solid angle 
corresponding to about 94% of the azimuth and 
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Figure 51. Elevation view of the nominal DIRC system geometry. For clarity, the end plug is not shown. 
All dimensions are given in mm. 



83% of the cm. polar angle cosine. 

The distance from the end of the bar to the 
PMTs is ~ 1.17 m, which together with the size 
of the bars and PMTs, gives a geometric contri- 
bution to the single photon Cherenkov angle res- 
olution of ^ 7mrad. This value is slightly larger 
than the rms spread of the photon production 
(dominated by a ~ 5.4mrad chromatic term) and 
transmission dispersions. The overall single pho- 
ton resolution is estimated to be about lOmrad. 

8.3.1. Cherenkov Photon Detection 
Efficiency 

Figure |5^ shows the contribution of various op- 
tical and electronic components of the DIRC to 
the Cherenkov photon detection efhciency as a 
function of wavelength. The data points per- 
tain to a particle entering the center of the bar 
at 90°. A typical design goal for the photon 
transport in the bar was that no single compo- 
nent should contribute more than 10-20% loss of 
detection efficiency. Satisfying this criterion re- 
quired an extremely high internal reflection coef- 
ficient of the bar surfaces (greater than 0.9992 per 
bounce), so that about 80% of the light is main- 



tained after multiple bounces along the bars (365 
bounces in the example of Figure |5^) . The ultra- 
violet cut-off is at ~ 300 nm, determined by the 
epoxy (Epotek 301-2 used to glue the fused 
silica bars together. The dominant contributor 
to the overall detection efficiency is the quantum 
efficiency of the PMT. Taking into account ad- 
ditional wavelength independent factors, such as 
the PMT packing fraction and the geometrical 
efficiency for trapping Cherenkov photons in the 
fused silica bars via total internal reflection, the 
expected number of photoelectrons (iVpe) is ~ 28 
for a = 1 particle entering normal to the sur- 
face at the center of a bar, and increases by over 
a factor of two in the forward and backward di- 
rections. 

8.3.2. DIRC Water System 

The DIRC water system is designed to main- 
tain good transparency at wavelengths as small as 
300 nm. One way to achieve this is to use ultra- 
pure, de-ionized water, close to the theoretical 
limit of 18 Mi7cm resistivity. In addition, the wa- 
ter must be de-gassed and the entire system kept 
free of bacteria. Purified water is aggressive in at- 
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Figure 53. Transmission, reflectivity and quan- 
tum efliciency for various components of the 
DIRC as a function of wavelength for a /3 = 1 
particle at normal incidence to the center of a 
bar 0. 

tacking many materials, and those in contact with 
the water were selected based on known compata- 
bility with purified water. To maintain the neces- 
sary level of water quality, most plumbing compo- 
nents are made of stainless steel or polyvinylidene 
fluoride. 

The system contains an input line with six 
mechanical filters (three 10 /xm, two 5 /im, and 
one 1 /im), a reverse osmosis unit, de- ionization 
beds, a Teflon microtube de-gasser and various 
pumps and valves. To prevent bacteria growth, 
it is equipped with a UV lamp (254 nm wave- 
length) and fllters (two 1 fim, two 0.2 /xm, and 
charcoal filters). Sampling ports are provided to 
check the water quality and to monitor resistivity. 



pH-value, temperature, and flow. A gravity feed 
return system prevents overpressure. The entire 
standoff box water volume can be recirculated up 
to four times a day. 

The operating experience with the water sys- 
tem so far has been very good. The water volume 
is exchanged every ten hours and the resistivity 
of the water is typically 18 MJ7cm in the sup- 
ply line and 8-10 Mf2cm in the return line at a 
temperature of about 23-26° C. The pH- value is 
about 6.5 and 6.6-6.7 in the supply and return 
water, respectively. The water transparency is 
routinely measured using lasers of three differ- 
ent wavelengths. The transmission is better than 
92% per meter at 266 nm and exceeds 98% per 
meter at 325 nm and 442 nm. 

Potential leaks from the water seals between 
the bar boxes and the standoff box are detected 
by a water leak detection system of 20 custom wa- 
ter sensors in and about the bar box slots. Two 
commercial ultrasonic flow sensors are used to 
monitor water flow in two (normally dry) drain 
lines in addition to the 12 humidity sensors on a 
nitrogen gas output line from each bar box (see 
below). Should water be detected, a valve in a 
100 mm diameter drain line is opened, and the 
entire system is drained in about 12 minutes. 

All elements inside the standoff box (PMT, 
plastic PMT housing, gaskets, light catchers) 
were tested at normal and elevated temperatures 
to withstand the highly corrosive action of ultra- 
pure water and to prevent its pollution. For in- 
stance, rhodium-plated mirrors on ULTEM sup- 
port had to be used for the light catchers fs^ . 

8.3.3. DIRC Gas System 

Nitrogen gas from liquid nitrogen boil-off is 
used to prevent moisture from condensing on the 
bars, and used also to detect water leaks. The 
gas flows through each bar box at the rate of 
100-200 cm'^/min, and is monitored for humid- 
ity to ensure that the water seal around the bar 
box remains tight. The gas is filtered through a 
molecular sieve and three mechanical filters to re- 
move particulates (7 /xm, 0.5 /im, and 0.01 fim). 
Dew points of the gas returned from the bar boxes 
are about -40°C. Approximately one third of the 
input nitrogen gas leaks from the bar boxes and 
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keeps the bar box slots in the mechanical support 
structure free of condensation. 

8.4. Electronics 

8.4.1. DIRC PMT Electronics 

The DIRC PMT base contains a single printed 
circuit board, equipped with surface mounted 
components. The operating high voltage (HV) 
of the PMTs is ~ 1.14 kV, with a range between 
0.9 and 1.3 kV. Groups of 16 tubes are selected 
for uniformity of gain to allow their operation at 
a common HV provided from a single distribution 
board. 

The HV is provided by a CAEN SY-527 high 
voltage distribution system. Each of the 12 sec- 
tors receives HV through 56 high voltage chan- 
nels, distributed through a single cable bundle. 
Each voltage can be set between and 1.6 kV. 

8.4.2. DIRC Front-End Electronics 

The DIRC front-end electronics (FEE) is de- 
signed to measure the arrival time of each 
Cherenkov photon detected by the PMT ar- 
ray 1 65 1 to an accuracy that is limited by the in- 
trinsic 1.5 ns transit time spread of the PMTs. 
The design contains a pipeline to deal with the 
LI trigger latency of 12/is, and can handle ran- 
dom background rates of up to 200 kHz/PMT 
with zero dead time. In addition, the pulse height 
spectra can be measured to ensure that the PMTs 
operate on the HV plateau. However, because the 
ADC information is not needed to reconstruct 
events, 64 PMTs are multiplexed onto a single 
ADC for monitoring and calibration. 

The DIRC FEE is mounted on the outside of 
the standoff box and is highly integrated in or- 
der to minimize cable lengths and to retain the 
required single photoelectron sensitivity. Each 
of the 168 DIRC Front-end Boards (DFBs) pro- 
cesses 64 PMT inputs, containing eight custom 
analog chips along with their associated level 
translators, four custom-made TDC ICs, one 8- 
bit flash ADC (FADC), two digitally controlled 
calibration signal generators, multi-event buffers, 
and test hardware. 

The PMT signals are amplified, and pulse- 
shaped by an eight-channel analog IC A 
digital pulse timed with the peak of the input 



pulse is output by a zero-crossing discriminator, 
as well as a pulse shaped by a CR-RC filter with 
80 ns peaking time, which was chosen to allow 
for the ADC multiplexing. The multiplexer se- 
lects the channel to be digitized by the FADC for 
calibration. 

The TDC IC is a 16-channel TDC with 
0.5 ns binning, input buffering, and selective read- 
out of the data in time with the trigger. To cope 
with the LI maximum trigger latency of 12 /xs 
and jitter of 1 /is, the selective readout process ex- 
tracts data in time with the trigger within a pro- 
grammable time window. The acceptance win- 
dow width is programmable between 64 ns and 
2 IIS and is typically set at 600 ns. The twelve 
DIRC Crate Controllers (DCCs) that form the 
interface to the VME front-end crates are con- 
nected to six ROMs via 1.2 Gbit/s optical fibers. 

8.4.3. DAQ Feature Extraction 

Raw data from the DFBs are processed in 
the ROMs by a feature extraction algorithm be- 
fore being transmitted to the segment and event 
builder. This software algorithm reduces the 
data volume by roughly 50% under typical back- 
ground conditions. DFB data that contain er- 
rors are flagged and discarded. The only data 
errors seen to date have been traced to dam- 
aged DFBs that were replaced immediately. Be- 
cause the dataflow system can reliably transmit 
at most 32 kBytes/crate, the feature extraction 
must sometimes truncate data to limit the event 
size. Event data are replaced with a per-DFB 
occupancy summary when a ROM's hit occu- 
pancy exceeds 56%, which occurs about once in 
10* events. An appropriate flag is inserted into 
the feature extraction output whenever trunca- 
tion or deletion occurs. Errors, truncation, and 
feature extraction performance are continuously 
monitored online, and exceptions are either im- 
mediately corrected or logged for future action. 

8.4.4. DIRC Calibration 

The DIRC uses two independent approaches for 
a calibration of the unknown PMT time response 
and the delays introduced by the FEE and the 
fast control system. The first is a conventional 
pulser calibration. The second uses reconstructed 
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tracks from collision data. 

The pulser calibration is performed online us- 
ing a light pulser system which generates precisely 
timed 1 ns duration light pulses from twelve, blue 
LEDs, one per sector. The LEDs are triggered 
by the global fast control calibration strobe com- 
mand sent to the DCCs. The DCC triggers an 
individual LED for each sector upon receipt of 
calibration strobe. Pulses in adjacent sectors are 
staggered by 50 ns to prevent light crosstalk be- 
tween sectors. The pulser is run at roughly 2 kHz 
for the time delay calibration. The LED light is 
transmitted through approximately 47 m-long op- 
tical fibers to diffusers mounted on the inner sur- 
face of the standoff box wall opposite the PMTs. 
This light produces about 10% photoelectron oc- 
cupancy nearly uniformly throughout the stand- 
off box. 

Histograms of TDC times for each PMT are ac- 
cumulated in parallel in the ROMs, and then fit 
to an asymmetric peak function. About 65,000 
light pulses are used to determine the mean time 
delay of each of the PMTs in the standoff box to 
a statistical accuracy of better than 0.1 ns. The 
LED pulser is also used to monitor the photo- 
tube gains using the ADC readout. As with the 
TDC calibration, histograms and fits of the ADC 
spectrum are accumulated and fit in the ROM. 
A calibration run including both TDC and ADC 
information for all PMTs requires a few minutes, 
and is run once per day. Daily calibrations not 
only verify the time delays, but allow the detec- 
tion of hardware failures. 

The data stream calibration uses reconstructed 
tracks from the collision data. For calibration of 
the global time delay, the observed, uncalibrated 
times minus the expected arrival times, Ai^, are 
collected during the online prompt reconstruction 
processing. To calculate individual channel cali- 
brations, Af^ values for each DIRC channel are 
accumulated until statistics equivalent to 100,000 
tracks are collected. The distribution for each 
channel is fit to extract the global time offset cal- 
ibration. 

The data stream and online pulser calibrations 
of the electronic delays, and of the PMT time 
response and gain yield fully consistent results, 
although the data stream results in 15% better 



timing resolution than the pulser calibration. The 
time delay values per channel are typically stable 
to an rms of less than 0.1 ns over more than one 
year of daily calibrations. 

8.4.5. DIRC Environmental Monitoring 
System 

The DIRC environmental monitoring system is 
divided into three parts, corresponding to three 
separate tasks. The first deals with the control 
and monitoring of the HV system for the PMTs. 
The second is devoted to monitoring low voltages 
related to the FEE. The third controls a variety 
of other detector parameter settings. An inter- 
lock system, based on a standard VME module 
(SIAM), is provided. For the purposes of the 
DIRC, three dedicated VME CPUs run the appli- 
cation code. The communication between the HV 
mainframes and the monitoring crate is achieved 
by a CAENET controller (V288). The HV mon- 
itor task controls the step sizes for ramping the 
HV up or down, as well as the communication of 
alarm conditions, and the values and limits for 
the HV and current of each channel. 

The purpose of FEE monitoring is to con- 
trol and monitor parameters related to the 
FEE. For each DIRC sector, a custom multi- 
purpose board, the DCC, equipped with a micro- 
controller | |68t incorporating the appropriate com- 
munication protocol (CANbus), is situated in the 
same crate as the DEB. All monitoring and con- 
trol tasks are implemented on this card. The pa- 
rameters monitored are the low voltages for the 
DFBs and DCCs, the status of the optical link 
(Finisar), the temperature on supply boards, and 
the VME crate status. 

The third part of the monitoring system is 
based on a custom ADC VME board (VSAM), 
used to monitor various type of sensors: magnetic 
field sensors, an ensemble of 12 beam monitoring 
scalers, 16 Csl radiation monitors, the water level 
in the standoff box as well as its pH-value, resis- 
tivity, and temperature. 

8.5. Operational Issues 

The DIRC was successfully commissioned and 
attained performance close to that expected from 
Monte Carlo simulation. The DIRC has been 
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Figure 54. Display of an e+e" event reconstructed in BABAR with two different time cuts. On 

the left, all DIRC PMTs with signals within the ±300 ns trigger window are shown. On the right, only 
those PMTs with signals within 8 ns of the expected Cherenkov photon arrival time are displayed. 



robust and stable, and, indeed, serves also as 
a background detector for PEP-II tuning. Fig- 
ure [5^ shows a typical di-muon event (e+e~ — > 
fj,'^fi~). In addition to the signals caused by the 
Cherenkov light from the two tracks, about 500 
background signals can be seen in the readout 
window of ±300 ns. This background is dom- 
inated by low energy photons from the PEP-II 
machine hitting the standoff box. Some care in 
machine tuning is required to stay under a noise 
limit of about 200 kHz/tube imposed by limited 
DAQ throughput. Lead shielding has been in- 
stalled around the beam line components just 
outside the backward cndcap, and has substan- 
tially reduced this background. 

After about two years of running, approxi- 
mately 99.7% of PMTs and electronic channels 
are operating with nominal performance. 

Some deterioration of the PMT front glass win- 
dows (made of B53 Borosilicate glass) that are im- 
mersed in the ultra-pure water of the standoff box 
has been observed. For most of the tubes, the ob- 
servable effect is typically a slight cloudiness, but 
for about 50 tubes, it is much more pronounced. 
Extensive R&D has demonstrated that the effect 



is associated with a loss of sodium and boron from 
the surface of the glass For most tubes, the 
leaching rate is a few microns per year, and is ex- 
pected to be acceptable for the full projected ten 
year lifetime of the experiment. However, for the 
~ 50 tubes, the incorrect glass was used by the 
manufacturer. That glass does not contain zinc, 
making it much more susceptible to rapid leach- 
ing. This leaching may eventually lead to either 
a loss of performance, or some risk of mechanical 
failure of the face plates for these tubes. Direct 
measurements of the number of Cherenkov pho- 
tons observed in di-muon events as a function of 
time suggest that the total loss of photons from 
all sources is less than 2%/year. 

8.6. Data Analysis and Performance 

Figure ^ shows the pattern of Cherenkov 
photons in a di-muon event, before and after 
reconstruction. The time distribution of real 
Cherenkov photons from a single event is of order 
^ 50 ns wide, and during normal data-taking they 
are accompanied by hundreds of random photons 
in a flat background within the trigger acceptance 
window. Given a track pointing at a particular 
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fused silica bar and a candidate signal in a PMT 
within the optical phase space of that bar, the 
Cherenkov angle is determined up to a 16- fold 
ambiguity: top or bottom, left or right, forward 
or backward, and wedge or no-wedge reflections. 
The goal of the reconstruction program is to as- 
sociate the correct track with the candidate PMT 
signal, with the requirement that the transit time 
of the photon from its creation in the bar to its 
detection at the PMT be consistent with the mea- 
surement error of ~ 1.5 ns. 

8.6.1. Reconstruction 

An unbinncd maximum likelihood formalism is 
used to incorporate all information provided by 
the space and time measurements from the DIRC. 

The emission angle and the arrival time of 
the Cherenkov photons are reconstructed from 
the observed space-time coordinates of the PMT 
signals, transformed into the Cherenkov coordi- 
nate system {9c, (fic, and St) as follows: The 
known spatial position of the bar through which 
the track passed and the PMTs whose signal 
times lie within the readout window of ±300 ns 
from the trigger are used to calculate the three- 
dimensional vector pointing from the center of the 
bar end to the center of each tube. This vector 
is then extrapolated into the radiator bar (using 
Snell's law). This procedure defines, up to the 
16-fold ambiguity described above, the angles dc 
and 4>c of a photon. 

The DIRC time measurement represents the 
third dimension of the photomultiplier hit recon- 
struction. The timing resolution is not competi- 
tive with the position information for Cherenkov 
angle reconstruction, but timing information is 
used to suppress background hits from the beam 
induced background and, more importantly, ex- 
clude other tracks in the same event as the source 
of the photon. Timing information is also used to 
resolve the forward-backward and wedge ambigu- 
ities in the hit-to-track association. 

The relevant observable to distinguish between 
signal and background photons is the difference 
between the measured and expected photon ar- 
rival time, At-y. It is calculated for each photon 
using the track time-of-flight (assuming it to be 
a charged pion), the measured time of the can- 



didate signal in the PMT and the photon prop- 
agation time within the bar and the water filled 
standoff box. The time information and the re- 
quirement of using only physically possible pho- 
ton propagation paths reduce the number of am- 
biguities from 16 to typically 3. Applying the 
time information also substantially improves the 
correct matching of photons with tracks and re- 
duces the number of accelerator induced back- 
ground signals by approximately a factor 40, as 
illustrated in Figure HJ. 

The reconstruction routine currently provides a 
likelihood value for each of the five stable particle 
types (e,/i,7r,K,p) if the track passes through the 
active volume of the DIRC. These likelihoods are 
calculated in an iterative process by maximizing 
the likelihood value for the entire event while test- 
ing different hypotheses for each track. If enough 
photons are found, a fit of 9c and the number of 
observed signal and background photons are cal- 
culated for each track. 

8.6.2. Results 

The parameters of expected DIRC performance 
were derived from extensive studies with a va- 
riety of prototypes, culminating with a full-size 
prototype in a beam at CERN ||7^. The test 
beam results were well-described by Monte Carlo 
simulations of the detector. The performance of 
the full detector is close to expectations, and ad- 
ditional offline work, particularly on geometrical 
alignment, is expected to lead to further improve- 
ments. 

In the absence of correlated systematic errors, 
the resolution (crc,track) on the track Cherenkov 
angle should scale as 

CTCtrack = rivj— ' ^ ) 

where crc,7 is the single photon Cherenkov an- 
gle resolution, and Npc is the number of photons 
detected. Figure ^(a) shows the single photon 
angular resolution obtained from di-muon events. 
There is a broad background of less than 10% 
relative height under the peak, that originates 
mostly from track-associated sources, such as S 
rays, and combinatorial background. The width 
of the peak translates to a resolution of about 
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Figure 55. The difference between (a) the mea- 
sured and expected Cherenkov angle for single 
photons, A9c.-y, and (b) the measured and ex- 
pected photon arrival time, for single muons in 
fX^ IJL~ events. 



10.2 mrad, in good agreement with the expected 
value. The measured time resolution (see Fig- 
ure ^(b)) is 1.7 ns, close to the intrinsic 1.5 ns 
transit time spread of the PMTs. 

The number of photoelectrons shown in Fig- 
ure ^ varies between 20 for small polar angles at 
the center of the barrel and 65 at large polar an- 
gles. This variation is well reproduced by Monte 
Carlo simulation, and can be understood from the 
geometry of the DIRC. The number of Cherenkov 
photons varies with the pathlength of the track 
in the radiator, it is smallest at perpendicular in- 
cidence at the center and increases towards the 
ends of the bars. In addition, the fraction of pho- 
tons trapped by total internal reflection rises with 
larger values of | cos 6'track | • This increase in the 
number of photons for forward going tracks is a 
good match to the increase in momentum and 
thus benefits the DIRC performance. 

With the present alignment, the track 
Cherenkov angle resolution for di-muon events is 
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Figure 56. Number of detected photons versus 
track polar angle for reconstructed tracks in di- 
muon events compared to Monte Carlos simula- 
tion. The mean number of photons in the simu- 
lation has been tuned to match the data. 

shown in Figure The width of the fitted Gaus- 
sian distribution is 2.5 mrad compared to the de- 
sign goal of 2.2 mrad. From the measured sin- 
gle track resolution versus momentum in di-muon 
events and the difference between the expected 
Cherenkov angles of charged pions and kaons, the 
pion-kaon separation power of the DIRC can be 
inferred. As shown in Figure 58, the expected 



separation between kaons and pions at 3 GeV/c 
is about 4.2ct, within 15% of the design goal. 

Figure |5^ shows an example of the use of the 
DIRC for particle identification. The Kir invari- 
ant mass spectra are shown with and without the 
use of the DIRC for kaon identification. The peak 
corresponds to the decay of the particle. 

The efficiency for correctly identifying a 
charged kaon that traverses a radiator bar and 
the probability to wrongly identify a pion as a 
kaon are determined using — > K't:^ decays 
selected kinematically from inclusive D* produc- 
tion and are shown as a function of the track 
momentum in Figure ^ for a particular choice 
of particle selection criteria. The mean kaon se- 
lection efficiency and pion misidcntification are 
96.2 ±0.2% (stat.) and 2.1 ±0.1% (stat.), respec- 
tively. 
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(mrad) 



Figure 57. The difference between the measured 
and expected Cherenkov angle, A9c,track, for sin- 
gle muons in fJ.'^ fi~ events. The curve represents a 
Gaussian distribution fit to the data with a width 
of 2.5 mrad. 
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Figure 59. Invariant Ktt inclusive mass spectrum 
with and without the use of the DIRC for kaon 
identification. The mass peak corresponds to the 
decay of the particle. 
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Figure 58. Expected tt-K separation in 
Tr~^TT~ events versus track momentum inferred 
from the measured Cherenkov angle resolution 
and number of Cherenkov photons per track in 
di-muon events. 



Figure 60. Efficiency and misidentification prob- 
ability for the selection of charged kaons as a 
function of track momentum, determined using 
I?° K'TT^ decays selected kinematically from 
inclusive D* production. 
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8.7. Summary 

The DIRC is a novel ring-imaging Cherenkov 
detector that is well-matched to the hadronic PID 
requirements of BABAR. The DIRC has been ro- 
bust and stable and, two years after installation, 
about 99.7% of all PMTs and electronic channels 
are operating with nominal performance. Addi- 
tional shielding in the standoff box tunnel region 
should reduce the sensitivity to beam-induced 
backgrounds, as should faster FEE, both installed 
during the winter 2000-2001 shutdown. At lumi- 
nosities around 1 x 10'^* cni~^s~^, the TDC IC 
will have to be replaced with a faster version and 
deeper buffering. The design process for this is 
underway. 

The detector performance achieved is rather 
close to that predicted by the Monte Carlo sim- 
ulations. Alignment and further code develop- 
ments are expected to further improve perfor- 
mance. 

9. Electromagnetic Calorimeter 

9.1. Purpose and Design 

The electromagnetic calorimeter (EMC) is de- 
signed to measure electromagnetic showers with 
excellent efficiency, and energy and angular res- 
olution over the energy range from 20MeV to 
9 GeV. This capability allows the detection of 
photons from tt'^ and rj decays as well as from elec- 
tromagnetic and radiative processes. By identify- 
ing electrons, the EMC contributes to the flavor 
tagging of neutral B mesons via semi-leptonic de- 
cays, to the reconstruction of vector mesons like 
J/ip , and to the study of semi-leptonic and rare 
decays of B and D mesons, and t leptons. The 
upper bound of the energy range is set by the need 
to measure QED processes, like e~^e~ e'^e~{^) 
and e+e~ — > 77, for calibration and luminos- 
ity determination. The lower bound is set by 
the need for highly efficient reconstruction of B- 
meson decays containing multiple tt^s and 77^8. 

9.1.1. Requirements 

The measurement of extremely rare decays of 
B mesons containing 7r°s {e.g., B^ tt'^tt^) poses 
the most stringent requirements on energy reso- 
lution, namely of order 1-2%. Below energies of 



2 GeV, the tt" mass resolution is dominated by 
the energy resolution. At higher energies, the an- 
gular resolution becomes dominant, and therefore 
it is required to be of the order of a few mrad. 

Furthermore, the EMC has to be compatible 
with the 1.5 T field of the solenoid and operate re- 
liably over the anticipated ten- year lifetime of the 
experiment. To achieve excellent resolution, sta- 
ble operating conditions have to be maintained. 
Temperatures and the radiation exposure must 
be closely monitored, and precise calibrations of 
the electronics and energy response over the full 
dynamic range must be performed frequently. 

9.1.2. Design Considerations 

The requirements stated above lead to the 
choice of a hermetic, total-absorption calorime- 
ter, composed of a finely segmented array of 
thalhum-doped cesium iodide (CsI(Tl)) crystals. 
The crystals are read out with silicon photodi- 
odes that are matched to the spectrum of scin- 
tillation light. Recent experience at CLEO ||7^ 
has demonstrated the suitability of this choice for 
physics at the T(45') resonance. 

The energy resolution of a homogeneous crystal 
calorimeter can be described empirically in terms 
of a sum of two terms added in quadrature 



(6) 



where E and ctb refer to the energy of a photon 
and its rms error, measured in GeV. The energy 
dependent term a arises primarily from the fluctu- 
ations in photon statistics, but it is also impacted 
by electronic noise of the photon detector and 
electronics. Furthermore, beam-generated back- 
ground will lead to large numbers of additional 
photons that add to the noise. This term is dom- 
inant at low energies. The constant term, b, is 
dominant at higher energies (> IGeV). It arises 
from non- uniformity in light collection, leakage or 
absorption in the material between and in front of 
the crystals, and uncertainties in the calibrations. 
Most of these effects can be influenced by design 
choices, and they are stable with time. Others 
will be impacted by changes in the operating con- 
ditions, like variations in temperature, electronics 
gain, and noise, as well as by radiation damage 
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caused by beam-generated radiation. 

The angular resolution is determined by the 
transverse crystal size and the distance from the 
interaction point. It can also be empirically pa- 
rameterized as a sum of an energy dependent and 
a constant term, 



Table 11 

Properties of CsI(Tl) 



ere 



(7) 



where the energy E is measured in GeV. The de- 
sign of the EMC required a careful optimization 
of a wide range of choices, including the crystal 
material and dimensions, the choice of the photon 
detector and readout electronics, and the design 
of a calibration and monitoring system. These 
choices were made on the basis of extensive stud- 
ies, prototyping and beam tests [Q, and Monte 
Carlo simulation, taking into account limitations 
of space and the impact of other BABAR detector 
systems. 

Under ideal conditions, values for the energy 
resolution parameters a and b close to 1-2% could 
be obtained. A position resolution of a few mm 
will translate into an angular resolution of a few 
mrad; corresponding parameter values are c ~ 
3 mrad and d ^ 1 mrad. 

However in practice, such performance is very 
difficult to achieve in a large system with a small, 
but unavoidable amount of inert material and 
gaps, limitations of electronics, and background 
in multi-particle events, plus contributions from 
beam-generated background. 

Though in CsI(Tl) the intrinsic efficiency for 
the detection of photons is close to 100% down 
to a few MeV, the minimum measurable en- 
ergy in colliding beam data is expected to be 
about 20 MeV, a limit that is largely determined 
by beam- and event-related background and the 
amount of material in front of the calorimeter. 
Because of the sensitivity of the 7r° efficiency to 
the minimum detectable photon energy, it is ex- 
tremely important to keep the amount of material 
in front of the EMC to the lowest possible level. 

9.1.3. CsI(Tl) Crystals 

Thallium-doped Csl meets the needs of BABAR 
in several ways. Its properties are listed in Ta- 
ble |ll|. The high light yield and small Molicre 



Parameter 


Values 


Hacliation Length. 


IOC 

l.oo cm 




3 8 cm 


Density 


4.53 g/cm^ 


Light Yield 


50,000 7/ MeV 


Light Yield Temp. Coeff. 


0.28%/°C 


Peak Emission Amax 


565 nm 


Refractive Index (Amax) 


1.80 


Signal Decay Time 


680 ns (64%) 




3.34 ^s (36%) 



radius allow for excellent energy and angular res- 
olution, while the short radiation length allows 
for shower containment at BABAR energies with a 
relatively compact design. Furthermore, the high 
light yield and the emission spectrum permit effi- 
cient use of silicon photodiodes which operate well 
in high magnetic fields. The transverse size of the 
crystals is chosen to be comparable to the Moliere 
radius achieving the required angular resolution 
at low energies while appropriately limiting the 
total number of crystals (and readout channels). 

9.2. Layout and Assembly 
9.2.1. Overall Layout 

The EMC consists of a cylindrical barrel and 
a conical forward endcap. It has full coverage in 
azimuth and extends in polar angle from 15.8° 
to 141.8° corresponding to a solid-angle cover- 
age of 90% in the cm. system (see Figure |6l] 
and Table ^2]). The barrel contains 5,760 crystals 
arranged in 48 distinct rings with 120 identical 
crystals each. The endcap holds 820 crystals ar- 
ranged in eight rings, adding up to a total of 6,580 
crystals. The crystals have a tapered trapezoidal 
cross section. The length of the crystals increases 
from 29.6 cm in the backward to 32.4 cm in the 
forward direction to limit the effects of shower 
leakage from increasingly higher energy particles. 

To minimize the probability of pre-showering, 
the crystals are supported at the outer radius, 
with only a thin gas seal at the front. The barrel 
and outer five rings of the endcap have less than 
0.3-0.6^0 of material in front of the crystal faces. 
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Figure 61. A longitudinal cross section of the EMC (only the top half is shown) indicating the arrangement 
of the 56 crystal rings. The detector is axially symmetric around the z-axis. All dimensions are given in 



Table 12 

Layout of the EMC, composed of 56 axially sym- 
metric rings, each consisting of Csl crystals of 
identical dimensions. 



Interval 


Length 


# 


Crystals 


(radians) 


(^o) 


Rings 


/Ring 


Barrel 


2.456 - 1.214 


16.0 


27 


120 


1.213-0.902 


16.5 


7 


120 


0.901 - 0.655 


17.0 


7 


120 


0.654 - 0.473 


17.5 


7 


120 


Endcap 


0.469 - 0.398 


17.5 


3 


120 


0.397- 0.327 


17.5 


3 


100 


0.326 - 0.301 


17.5 


1 


80 


0.300 - 0.277 


16.5 


1 


80 



The SVT support structure and electronics, as 
well as the Bl dipole shadow the inner three rings 
of the endcap, resulting in up to 3.0^0 for the 
innermost ring. The principal purpose of the two 
innermost rings is to enhance shower containment 
for particles close to the acceptance limit. 



9.2.2. Crystal Fabrication and Assembly 

The crystals were grown in boules from a melt 
of Csl salt doped with 0.1% thallium [||. They 
were cut from the boules, machined into tapered 
trapezoids (Figure |62|) to a tolerance of ±150 /im, 
and then polished |74|]. The transverse dimen- 
sions of the crystals for each of the 56 rings vary to 
achieve the required hermetic coverage. The typi- 
cal area of the front face is 4.7 x 4.7 cm^, while the 
back face area is typically 6.1 x 6.0 cm^. The crys- 
tals act not only as a total-absorption scintillating 
medium, but also as a light guide to collect light 
at the photodiodes that are mounted on the rear 
surface. At the polished crystal surface light is 
internally reflected, and a small fraction is trans- 
mitted. The transmitted light is recovered in part 
by wrapping the crystal with two layers of diffuse 
white reflector ||7^,]7^, each 165 /im thick. The 
uniformity of light yield along the wrapped crys- 
tal was measured by recording the signal from a 
highly coUimated radioactive source at 20 points 
along the length of the crystal. The light yield 
was required to be uniform to within ±2% in the 
front half of the crystal; the limit increased lin- 
early up to a maximum of ±5% at the rear face. 
Adjustments were made on individual crystals to 
meet these criteria by selectively roughing or pol- 
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Figure 62. A schematic of the wrapped 
CsI(Tl) crystal and the front-end readout package 
mounted on the rear face. Also indicated is the 
tapered, trapezoidal CFC compartment, which is 
open at the front. This drawing is not to scale. 



ishing the crystal surface to reduce or increase its 
reflectivity. 

Following these checks, the crystals were fur- 
ther wrapped in 25 /.tm thick aluminum foil which 
was electrically connected to the metal housing of 
the photodiode-preamplifier assembly to provide 
a Faraday shield. The crystals were covered on 
the outside with a 13 /im-thick layer of mylar to 
assure electrical isolation from the external sup- 
port. 

9.2.3. Photodiodes 

and Preamplifier Assembly 

The photon detector consists of two 2x1 cm^ 
silicon PIN diodes glued to a transparent 1.2 mm- 



thick polysterene substrate that, in turn, is glued 
to the center of the rear face of the crystal by 
an optical epoxy | |77t to maximize light transmis- 
sion . The surrounding area of the crystal face 
is covered by a plastic plate coated with white 
reflective paint |7^. The plate has two 3 mm- 
diameter penetrations for the fibers of the light 
pulser monitoring system. 

As part of the quality control process, the 
1.836 MeV photon line from a ^^Y radioactive 
source was used to measure the light yield of every 
crystal-diode assembly, employing a preamplifier 
with 2 /zs Gaussian shaping. The resulting signal 
distribution had a mean and rms width of 7300 
and 890 photoelectrons/MeV, respectively; none 
of the crystals had a signal of less than 4600 pho- 
toelectrons/MeV 

Each of the diodes is directly connected to a 
low-noise preamplifier. The entire assembly is en- 
closed by an aluminum fixture as shown in Fig- 
ure This fixture is electrically coupled to the 
aluminum foil wrapped around the crystal and 
thermally coupled to the support frame to dissi- 
pate the heat load from the preamplifiers. 

Extensive aging tests were performed to ascer- 
tain that the diodes and the preamplifiers met the 
ten-year lifetime requirements. In addition, daily 
thermal cycles of ±5°C were run for many months 
to assure that the diode-crystal epoxy joint could 
sustain modest temperature variations. 

9.2.4. Crystal Support Structure 

The crystals are inserted into modules that are 
supported individually from an external support 
structure. This structure is built in three sec- 
tions, a cylinder for the barrel and two semi- 
circular structures for the forward endcap. The 
barrel support cylinder carries the load of the bar- 
rel modules plus the forward endcap to the mag- 
net iron through four flexible supports. These 
supports decouple and dampen any acceleration 
induced by movements of the magnet iron during 
a potential earthquake. 

The modules are built from tapered, trape- 
zoidal compartments made from carbon-fiber- 
epoxy composite (CFC) with 300 /im-thick walls 
(Figure ^). Each compartment loosely holds 
a single wrapped and instrumented crystal and 
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Figure 63. The EMC barrel support structure, with details on the modules and electronics crates (not 
to scale). 



thus assures that the forces on the crystal sur- 
faces never exceed its own weight. Each module 
is surrounded by an additional layer of 300 /im 
CFC to provide additional strength. The mod- 
ules are bonded to an aluminum strong-back that 
is mounted on the external support. This scheme 
minimizes inter-crystal materials while exerting 
minimal force on the crystal surfaces; this pre- 
vents deformations and surface degradation that 
could compromise performance. By supporting 
the modules at the back, the material in front of 
the crystals is kept to a minimum. 

The barrel section is divided into 280 sepa- 
rate modules, each holding 21 crystals (7 x 3 in 
6 X (p). After the insertion of the crystals, the 
aluminum readout frames, which also stiffen the 
module, are attached with thermally-conducting 
epoxy to each of the CFC compartments. The en- 
tire lOOkg-module is then bolted and again ther- 
mally epoxied to an aluminum strong-back. The 
strong-back contains alignment features as well 
as channels that couple into the cooling system. 



Each module was installed into the 2.5 cm-thick, 
4 m-long aluminum support cylinder, and subse- 
quently aligned. On each of the thick annular 
end-flanges this cylinder contains access ports for 
digitizing electronics crates with associated cool- 
ing channels, as well as mounting features and 
alignment dowels for the forward endcap. 

The endcap is constructed from 20 identical 
CFC modules (each with 41 crystals), individu- 
ally aligned and bolted to one of two semi-circular 
support structures. The endcap is split vertically 
into two halves to facilitate access to the central 
detector components. 

The entire calorimeter is surrounded by a dou- 
ble Faraday shield composed of two 1 mm-thick 
aluminum sheets so that the diodes and pream- 
plifiers are further shielded from external noise. 
This cage also serves as the environmental bar- 
rier, allowing the slightly hygroscopic crystals to 
reside in a dry, temperature controlled nitrogen 
atmosphere. 
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9.2.5. Cooling System 

The EMC is maintained at constant, accurately 
monitored temperature. Of particular concern 
are the stability of the photodiode leakage current 
which rises exponentially with temperature, and 
the large number of diode-crystal epoxy joints 
that could experience stress due to differential 
thermal expansion. In addition, the light yield 
of CsI(Tl) is weakly temperature dependent. 

The primary heat sources internal to 
the calorimeter are the preamplifiers 
(2 X 50mW/crystal) and the digitizing elec- 
tronics (3kW per end- flange). In the barrel, the 
preamplifier heat is removed by conduction to the 
module strong backs which are directly cooled 
by Fluorinert (polychlorotrifluoro-ethylene) [Q. 
The digitizing electronics are housed in 80 
mini-crates, each in contact with the end-flanges 
of the cylindrical support structure. These crates 
are indirectly cooled by chilled water pumped 
through channels milled into the end-flanges 
close to the inner and outer radii. A separate 
Fluorinert system in the endcap cools both the 
20 mini-crates of digitizing electronics and the 
preamplifiers. 



9.3. Electronics 

The EMC electronics system, shown schemat- 
ically in Figure is required to have negligible 
impact on the energy resolution of electromag- 
netic showers from 20 MeV to 9 GeV, while ac- 
commodating the use of a 6.13 MeV radioactive 
source for calibration. These requirements set a 
limit of less than 250 keV equivalent noise energy 
(ENE) per crystal and define an 18-bit effective 
dynamic range of the digitization scheme. For 
source calibrations, the least significant bit is set 
to 50 keV, while for colliding beam data it is set 
to 200 keV. To reach the required energy reso- 
lution at high energies, the coherent component 
has to be significantly smaller than the incoher- 
ent noise component. In addition, the impact of 
high rates of low energy (<5 MeV) beam- induced 
photon background needs to be minimized. 



9.3.1. Photodiode Readout 
and Preamplifiers 

The ENE is minimized by maximizing the light 
yield and collection, employing a highly efficient 
photon detector, and a low-noise electronic read- 
out. The PIN silicon photodiodes have a 
quantum efficiency of 85% for the CsI(Tl) scin- 
tillation light [Q. At a depletion voltage of 
70 V, their typical dark currents were measured 
to be 4nA for an average capacitance of 85 pF; 
the diodes are operated at a voltage of 50 V. 
The input capacitance to the preamplifier is min- 
imized by connecting the diodes to the preampli- 
fier with a very short cable. The preamplifier is a 
low-noise charge-sensitive amplifier implemented 
as a custom application specific integrated cir- 
cuit (ASIC) |8^. It shapes the signal and acts 
as a band-pass filter to remove high- and low- 
frequency noise components. The optimum shap- 
ing time for the CsI(Tl)-photodiode readout is 
2-3 /is, but a shorter time was chosen to reduce 
the probability of overlap with low-energy pho- 
tons from beam background. The commensurate 
degradation in noise performance is recovered by 
implementing a realtime digital signal-processing 
algorithm following digitization. 

To achieve the required operational reliabil- 
ity ||8^ for the inaccessible front-end readout com- 
ponents, two photodiodes were installed, each 
connected to a preamplifier. In addition, all com- 
ponents were carefully selected and subjected to 
rigorous tests, including a 72-hour burn-in of the 
preamplifiers at 70° C to avoid infant mortality. 
The dual signals are combined in the postam- 
plification/digitization circuits, installed in mini- 
crates at the end-flanges, a location that is acces- 
sible for maintenance. 

9.3.2. Postampliflcation, Digitization 
and Readout 

The two preamplifiers on each crystal, A and B, 
each provide amplification factors of 1 and 32 and 
thus reduce the dynamic range of the signal that 
is transmitted to the mini-crates to 13-bits. A 
custom auto-range encoding (CARE) circuit jsij 
further amplifies the signal to arrive at a total 
gain of 256, 32, 4 or 1 for four energy ranges, 
0-50 MeV, 50-400 MeV, 0.4-3.2 GeV, and 3.2- 
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Figure 64. Schematic diagram of the EMC readout electronics. 



13.0 GeV, respectively. The appropriate range is 
identified by a comparator and the signal is dig- 
itized by a 10-bit, 3.7 MHz ADC. Data from 24 
crystals are multiplexed onto a fiber-optic driver 
and sent serially at a rate of 1.5 Gbytes/s across a 
30 m-long optical fiber to the ROM. In the ROM, 
the continuous data stream is entered into a dig- 
ital pipeline. A correction for pedestal and gain 
is applied to each sample. The pipeline is then 
tapped to extract the input to the calorimeter 
trigger. 

Upon receipt of the LI Accept signal, data 
samples within a time window of ±l/is are se- 
lected for the feature extraction. Up to now, the 
calorimeter feature extraction algorithm performs 
a parabolic fit to the peak of the signal waveform 
to derive its energy and time. In the future, it is 
planned to employ a digital filter prior to the sig- 
nal fit to further reduce noise. For this filter algo- 
rithm, the frequency decomposition of an average 
signal pulse and the typical noise spectrum are 
measured for all channels and subsequently used 
to derive an optimum set of weights that maxi- 
mizes the signal-to-noise ratio. These weights are 
then applied to individual samples to obtain a 
filtered waveform. 

The magnitude of the electronic noise is mea- 
sured as the rms width of the pedestal distribu- 
tion as shown in Figure 65. The observed dis- 



tribution for all channels translates to an ENE 
of 230 keV and 440 keV with and without digital 
filtering; this result is comparable to design ex- 
pectations. Measurements of the auto-correlation 
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Figure 65. The distribution of equivalent noise 
energy (ENE) or all channels of the EMC with 
and without digital filtering. The data were 
recorded in the absence of beams by a random 
trigger. 



function indicate that the coherent noise com- 
ponent is negligible compared to the incoherent 
noise, except for regions where the preamplifiers 
saturate (see below). 

During data-taking, the data acquisition im- 
poses a single-crystal readout threshold in order 
to keep the data volume at an acceptable level. 
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This energy threshold is currently set to 1 MeV 
and during stable colliding beam conditions on 
average 1,000 crystals are read out (measured 
with 600 mA of e~ and 1100 mA of and a ran- 
dom clock trigger), corresponding to an average 
occupancy of 16%. The electronic noise accounts 
for about 10%, while the remaining signals orig- 
inate from beam-generated background (see Sec- 
tion ^) . A typical hadronic event contributes sig- 
nals in 150 crystals. 

9.3.3. Electronics Calibration 
and Linearity 

To measure pedestal offsets, determine the 
overall gain, and to remove non-linearities the 
FEE are calibrated by precision charge injection 
into the preamplifier input. Initially, residual 
non-linearities of up to 12% in limited regions 
near each of the range changes were observed and 
corrected for offline ||8^. These non-linearities 
were traced to oscillations on the ADC cards that 
have since been corrected. The correction re- 
sulted in markedly improved energy resolution 
at high energies. Residual non-linearities (typi- 
cally 2-4%) arise primarily from cross-talk, im- 
pacting both the electronics calibrations and the 
coUiding-beam data. The effect is largest at about 
630 MeV (950 MeV) in a high (low) gain preampli- 
fier channel, inducing a 2 MeV (6 MeV) cross-talk 
signal in an adjacent channel. The implemen- 
tation of an energy dependent correction is ex- 
pected to significantly reduce this small, remain- 
ing effect, and lead to a further improvement of 
the energy resolution. 

9.3.4. Electronics Reliability 

With the exception of minor cable damage 
during installation (leaving two channels inop- 
erative), the system of 13,160 readout channels 
has met its reliability requirements. After the 
replacement of a batch of failing optical-fiber 
drivers, the reliability of the digitizing electronics 
improved substantially, averaging channel losses 
of less than 0.1%. 

9.4. Energy Calibration 

The energy calibration of the EMC proceeds 
in two steps. First, the measured pulse height in 
each crystal has to be translated to the actual en- 



ergy deposited. Second, the energy deposited in 
a shower spreading over several adjacent crystals 
has to be related to the energy of the incident 
photon or electron by correcting for energy loss 
mostly due to leakage at the front and the rear, 
and absorption in the material between and in 
front of the crystals, as well as shower energy not 
associated with the cluster. 

The offline pattern recognition algorithm that 
groups adjacent cry stal s into clusters is described 
in detail in Section 9.6. 



9.4.1. Individual Crystal Calibration 

In spite of the careful selection and tuning of 
the individual crystals, their light yield varies sig- 
nificantly and is generally non-uniform. It also 
changes with time under the impact of beam- 
generated radiation. The absorbed dose is largest 
at the front of the crystal and results in increased 
attenuation of the transmitted scintillation light. 
The light yield must therefore be calibrated at 
different energies, corresponding to different av- 
erage shower penetration, to track the effects of 
the radiation damage. 

The calibration of the deposited energies is 
performed at two energies at opposite ends of 
the dynamic range, and these two measurements 
are combined by a logarithmic interpolation. A 
6.13 MeV radioactive photon source provides 
an absolute calibration at low energy, while at 
higher energies (3-9 GeV) the relation between 
polar angle and energy of from Bhabha events 
is exploited js^ . 

A flux of low-energy neutrons (4xl0^/s) is used 
to irradiate Fluorinert to produce photons of 
6.13 MeV via the reaction ^^F + n -^'^^ N + a, 
167V 0*+P, i^O* + 7. The activated 
^^N has a half-life of 7 seconds and thus does not 
cause radiation damage or long-term activation. 
The fluid is pumped at a rate of 125 £/s from the 
neutron generator to a manifold of thin-walled 
(0.5 mm) aluminum pipes that are moimted im- 
mediately in front of the crystals. At this loca- 
tion, the typical rate of photons is 40 Hz/crystal. 

Figure 66 shows a typical source spectrum that 
was derived from the raw data by employing a 
digital filter algorithm. For a 30-minute exposure, 
a statistical error of 0.35% is obtained, compared 
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Figure 66. A typical pulse- height spectrum 
recorded with the radioactive source to calibrate 
the single-crystal energy scale of the EMC. The 
spectrum shows the primary 6.13 MeV peak and 
two associated escape peaks at 5.62 MeV and 
5.11 MeV. The solid line represents a fit to the 
total spectrum, the dotted lines indicate the con- 
tributions from the three individual photon spec- 
tra. 

to a systematic uncertainty of less than 0.1%. 
This calibration is performed weekly. 

At high energies, single crystal calibration 
is performed with a pure sample of Bhabha 
events 1 88 . As a function of the polar angle of the 



, the deposited cluster energy is constrained to 
equal the prediction of a GEANT-based Monte 
Carlo simulation js^. For a large number of en- 
ergy clusters, a set of simultaneous linear equa- 
tions relates the measured to the expected energy 
and thus permits the determination of a gain con- 
stant for each crystal. In a 12-hour run at a lu- 
minosity of 3 X 10'^^cm~^s~^ some 200 per 
crystal can be accumulated, leading to a statisti- 
cal error of 0.35%. This calibration has been per- 
formed about once per month, and will be fully 
automated in the future. 



9.4.2. Cluster Energy Correction 

The correction for energy loss due to shower 
leakage and absorption is performed as a function 
of cluster energy and polar angle. At low energy 
(E < 0.8 GeV), it is derived from 7r° decays Q. 
The true energy of the photon is expressed as a 
product of the measured deposited energy and a 
correction function which depends on InE and 
cos 9. The algorithm constrains the two-photon 
mass to the nominal tt'^ mass and iterativcly finds 
the coefficients of the correction function. The 
typical corrections are of order 6 ± 1%. The un- 
certainty in the correction is due to systematic 
uncertainties in the background estimation and 
the fitting technique. 

At higher energy (0.8 < E < 9 GeV) the correc- 
tion is estimated from single-photon Monte Carlo 
simulations. A second technique using radiative 
Bhabha events |^l| is being developed. The beam 
energy and the precise track momenta of the e"*" 
and e~, together with the direction of the radia- 
tive photon, are used to fit the photon energy. 
This fitted value is compared to the measured 
photon energy to extract correction coefficients, 
again as a function of InE and cos 9. 

9.5. Monitoring 

9.5.1. Environmental Monitoring 

The temperature is monitored by 256 thermal 
sensors that are distributed over the calorimeter, 
and has been maintained at 20 ± 0.5°C. Dry ni- 
trogen is circulated throughout the detector to 
stabilize the relative humidity at 1 ± 0.5%. 

9.5.2. Light-Pulser System 

The light response of the individual crystals is 
measured daily using a light-pulser system | |9^ , 
|9^ . Spectrally filtered light from a xenon flash 
lamp is transmitted through optical fibers to the 
rear of each crystal. The light pulse is similar in 
spectrum, rise-time and shape to the scintillation 
light in the CsI(Tl) crystals. The pulses are var- 
ied in intensity by neutral-density filters, allowing 
a precise measurement of the linearity of light col- 
lection, conversion to charge, amplification, and 
digitization. The intensity is monitored pulse-to- 
pulse by comparison to a reference system with 
two radioactive sources, ^''^Am and ^''^Gd, that 
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Figure 67. Impact of beam-generated radiation 
on the CsI(Tl) crystals: a) the integrated dose 
measured with RadFETs placed in front of the 
crystals, b) the degradation in light yield mea- 
sured with the radioactive-source calibration sys- 
tem. 

are attached to a small CsI(Tl) crystal that is 
read out by both a photodiode and a photomul- 
tiplier tube. The system is stable to 0.15% over 
a period of one week and has proven to be very 
valuable in diagnosing problems. For example, 
the ability to accurately vary the light intensity 
led to the detection of non-linear response in the 
electronics Ml. 



9.5.3. Radiation Monitoring and Damage 

The radiation exposure is monitored by 56 
and 60 realtime integrating dosimeters (Rad- 
FETs) ||l8| placed in front of the barrel and end- 
cap crystals. In Figure |6^, the accumulated dose 
is compared to the observed loss in scintillation 
light, separately for the endcap, the forward, and 
the backward barrel. The dose appears to follow 
the integrated luminosity, approximately linearly. 
The light loss is greatest in the forward region 
corresponding to the area of highest integrated 
radiation dose. The size of the observed light loss 
is close to expectations, based on extensive irra- 
diation tests. 

9.6. Reconstruction Algorithms 

A typical electromagnetic shower spreads over 
many adjacent crystals, forming a cluster of en- 
ergy deposits. Pattern recognition algorithms 
have been developed to efhciently identify these 
clusters and to differentiate single clusters with 
one energy maximum from merged clusters with 
more than one local energy maximum, referred to 
as a bumps. Furthermore, the algorithms deter- 
mine whether a bump is generated by a charged 
or a neutral particle. 

Clusters are required to contain at least one 
seed crystal with an energy above lOMeV. Sur- 
rounding crystals are considered as part of the 
cluster if their energy exceeds a threshold of 
1 MeV, or if they are contiguous neighbors (in- 
cluding corners) of a crystal with at least 3 MeV. 
The value of the single crystal threshold is set by 
the data acquisition system in order to keep the 
data volume at an acceptable level, given the cur- 
rent level of electronics noise and beam-generated 
background. It is highly desirable to reduce this 
threshold since fluctuations in the effective energy 
loss at the edges of a shower cause a degradation 
in resolution, particularly at low energies. 

Local energy maxima are identified within a 
cluster by requiring that the candidate crys- 
tal have an energy, -BLocaiMax, which exceeds 
the energy of each of its neighbors, and sat- 
isfy the following condition: 0.5{N — 2.5) > 
i^NMaxZ-BLocaiMax, whcrc -^NMax IS the highest en- 
ergy of any of the neighboring N crystals with an 
energy above 2 MeV. 
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Clusters are divided into as many bumps as 
there are local maxima. An iterative algorithm 
is used to determine the energy of the bumps. 
Each crystal is given a weight, Wi, and the bump 
energy is defined as Sbump — J2i ^i^n where the 
sum runs over all crystals in the cluster. For a 
cluster with a single bump, the result is Wi = 1. 
For a cluster with multiple bumps, the crystal 
weight for each bump is calculated as 

exp(-2.5rj/rA/) 



500 



E,, 



J2j Ej exp(-2.5rj/rM) ' 



(8) 



where the index j runs over all crystals in the 
cluster, vm refers to the Moliere radius, and 
is the distance of the ith crystal from the cen- 
troid of the bump. At the outset, all weights are 
set to one. The process is then iterated, whereby 
the centroid position used in calculating is de- 
termined from the weights of the previous itera- 
tion, until the bump centroid position is stable to 
within a tolerance of 1 mm. 

The position of a bump is calculated us- 
ing a center-of-gravity method with logarithmic, 
rather than linear weights |^,^, Wi — 4.0 -I- 
lniJi/i?bump, where only crystals with positive 
weights, i.e., Ei > 0.0184 x_Ebump, are used in the 
calculation. This procedure emphasizes lower- 
energy crystals, while utilizing only those crystals 
that make up the core of the cluster. A system- 
atic bias of the calculated polar angle originates 
from the non-projectivity of the crystals. This 
bias is corrected by a simple offset of — 2.6mrad 
for e > 90° and -H2.6mrad for 6* < 90°. 

A bump is associated with a charged particle 
by projecting a track to the inner face of the 
calorimeter. The distance between the track im- 
pact point and the bump centroid is calculated, 
and if it is consistent with the angle and momen- 
tum of the track, the bump is associated with 
this charged particle. Otherwise, it is assumed to 
originate from a neutral particle. 

On average, 15.8 clusters are detected per 
hadronic event, of which 10.2 are not associated 
with charged particle tracks. At current oper- 
ating conditions, beam-induced background con- 
tributes on average 1.4 neutral clusters with en- 
ergies above 20 MeV. This number is significantly 
smaller than the average number of crystals with 
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Figure 68. The ratio of the EMC measured en- 
ergy to the expected energy for electrons from 
Bhabha scattering of 7.5 GeV/c. The solid line 
indicates a fit using a logarithmic function. 
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Figure 69. The energy resolution for the ECM 
measured for photons and electrons from various 
processes. The solid curve is a fit to Equation ^ 
and the shaded area denotes the rms error of the 
fit. 

energies above 10 MeV (see Section ||). 

9.7. Performance 

9.7.1. Energy Resolution 

At low energy, the energy resolution of the 
EMC is measured directly with the radioactive 
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Figure 70. The angular resolution of the EMC 
for photons from 7r° decays. The solid curve is a 
fit to Equation ^. 

source yielding ge/E = 5.0 ± 0.8% at 6.13 MeV 
(see Figure ^). At high energy, the resolution is 
derived from Bhabha scattering, where the energy 
of the detected shower can be predicted from the 
polar angle of the . The measured resolution is 
cfe/E = 1.9 ± 0.07% at 7.5 GeV (see Figure H). 
Figure ^ shows the energy resolution extracted 
from a variety of processes as a function of en- 
ergy. Below 2 GeV, the mass resolution of tt^ and 
•q mesons decaying into two photons of approx- 
imately equal energy is used to infer the EMC 
energy resolution |0. The decay xa — > J/i^l 
provides a measurement at an average energy of 
about 500 MeV, and measurements at high energy 
are derived from Bhabha scattering. A fit to the 
energy dependence results in 

^^(2.32 ±0.30)% 

Values of these fitted parameters are higher than 
the somewhat optimistic design expectations, but 
they agree with detailed Monte Carlo simulations 
which include the contributions from electronic 
noise and beam background, as well as the impact 
of the material and the energy thresholds. 

9.7.2. Angular Resolution 

The measurement of the angular resolution is 
based on the analysis of 7r° and rj decays to two 
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Figure 71. Invariant mass of two photons in BB 
events. The energies of the photons and the 7r° 
are required to exceed 30 MeV and 300 MeV, re- 
spectively. The solid line is a fit to the data. 

photons of approximately equal energy. The re- 
sult is presented in Figure |7^. The resolution 
varies between about 12 mrad at low energies and 
3 mrad at high energies. A fit to an empirical pa- 
rameterization of the energy dependence results 



<J0 = (Jcj, 

,3.87±0.07 



/E(GeV) 



0.00 ±0.04) mrad. (10) 



These fitted values are slightly better than would 
be expected from detailed Monte Carlo simula- 
tions. 



9.7.3. 



Mass and Width 



Figure 71 shows the two-photon invariant mass 



in BB events. The reconstructed 7r° mass is mea- 
sured to be 135.1 MeV/ and is stable to better 
than 1% over the full photon energy range. The 
width of 6.9 MeV/ agrees well with the predic- 
tion obtained from detailed Monte-Carlo simula- 
tions. In low-occupancy t^t~ events, the width 
is slightly smaller, 6.5 MeV/ c^, for iv^ energies be- 
low 1 GeV. A similar improvement is also ob- 
served in analyses using selected isolated photons 
in hadronic events. 
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9.7.4. Electron Identification 

Electrons are separated from charged hadrons 
primarily on the basis of the shower energy, lat- 
eral shower moments, and track momentum. In 
addition, the dE/dx energy loss in the DCH and 
the DIRC Cherenkov angle are required to be 
consistent with an electron. The most impor- 
tant variable for the discrimination of hadrons 
is the ratio of the shower energy to the track 
momentum (E/p). Figure ^ shows the effi- 
ciency for electron identification and the pion mis- 
identification probability as a function of momen- 
tum for two sets of selection criteria. The elec- 
tron efficiency is measured using radiative Bhab- 



has and 



events. The pion 



misidentification probability is measured for se- 
lected charged pions from K'^ decays and three- 
prong T decays. A tight (very tight) selector re- 
sults in an efficiency plateau at 94.8% (88.1%) 
in the momentum range 0.5 < p < 2 GeV/c. The 
pion misidentification probability is of order 0.3% 
(0.15%) for the tight (very tight) selection crite- 



9.8. Summary 

The EMC is presently performing close to de- 
sign expectations. Improvements in the energy 
resolution are expected from the optimization of 
the feature-extraction algorithms designed to fur- 
ther reduce the electronics noise. Modifications to 
the electronics should allow for more precise cal- 
ibrations. The expected noise reduction should 
permit a lower single-crystal readout threshold. 
However, this decrease in noise might be offset 
by an increase in the beam background that is 
expected for higher luminosities and beam cur- 
rents. 
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Figure 72. The electron efficiency and pion mis- 
identification probability as a function of a) the 
particle momentum and b) the polar angle, mea- 
sured in the laboratory system. 



10. Detector for Muons and Neutral 
Hadrons 

10.1. Physics Requirements and Goals 

The Instrumented Flux Return (IFR) was de- 
signed to identify muons with high efficiency and 
good purity, and to detect neutral hadrons (pri- 
marily K1 and neutrons) over a wide range of 
momenta and angles. Muons are important for 
tagging the fiavor of neutral B mesons via semi- 



leptonic decays, for the reconstruction of vector 
mesons, like the J/ip , and for the study of semi- 
leptonic and rare decays involving leptons of B 
and D mesons and r leptons. detection al- 
lows the study of exclusive B decays, in particular 
CP eigenstates. The IFR can also help in vetoing 
charm decays and improve the reconstruction of 
neutrinos. 

The principal requirements for IFR are large 
solid angle coverage, good efficiency, and high 



75 



background rejection for muons down to mo- 
menta below 1 GeV/c. For neutral hadrons, high 
efhciency and good angular resolution are most 
important. Because this system is very large and 
difScult to access, high reliability and extensive 
monitoring of the detector performance and the 
associated electronics plus the voltage distribu- 
tion are required. 

10.2. Overview and RPC Concept 

The IFR uses the steel flux return of the mag- 
net as a muon filter and hadron absorber. Sin- 
gle gap resistive plate chambers (RPCs) |^] with 
two-coordinate readout have been chosen as de- 
tectors. 

The RPCs are installed in the gaps of the finely 
segmented steel (see Section ^) of the barrel and 
the end doors of the flux return, as illustrated in 
Figure |7^. The steel segmentation has been cho- 
sen on the basis of Monte Carlo studies of muon 
penetration and charged and neutral hadron in- 
teractions. The steel is segmented into 18 plates, 
increasing in thickness from 2 cm for the inner 
nine plates to 10 cm for the outermost plates. The 
nominal gap between the steel plates is 3.5 cm 
in the inner layers of the barrel and 3.2 cm else- 
where. There are 19 RPC layers in the barrel 
and 18 in the endcaps. In addition, two layers of 
cylindrical RPCs are installed between the EMC 
and the magnet cryostat to detect particles exit- 
ing the EMC. 

RPCs detect streamers from ionizing particles 
via capacitive readout strips. They offer several 
advantages: simple, low cost construction and the 
possibility of covering odd shapes with minimal 
dead space. Further benefits are large signals and 
fast response allowing for simple and robust front- 
end electronics and good time resolution, typi- 
cally 1-2 ns. The position resolution depends on 
the segmentation of the readout; a value of a few 
mm is achievable. 

The construction of the planar and cylindrical 
RPCs differ in detail, but they are based on the 
same concept. A cross section of an RPC is shown 
schematically in Figure [tJ. 

The planar RPCs consist of two bakelite (phe- 
nolic polymer) sheets, 2 mm-thick and separated 
by a gap of 2 mm. The gap is enclosed at the 
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Figure 74. Cross section of a planar RPC with the 
schematics of the high voltage (HV) connection. 

edge by a 7 mm wide frame. The gap width is 
kept uniform by polycarbonate spacers (0.8 cm^) 
that are glued to the bakelite, spaced at dis- 
tances of about 10 cm. The bulk resistivity of 
the bakelite sheets has been especially tuned to 
10^^-10^^ Hem. The external surfaces are coated 
with graphite to achieve a surface resistivity of ~ 
100kJ7/square. These two graphite surfaces are 
connected to high voltage 8kV) and ground, 
and protected by an insulating mylar film. The 
bakelite surfaces facing the gap are treated with 
linseed oil. The RPCs are operated in limited 
streamer mode and the signals are read out ca- 
pacitively, on both sides of the gap, by external 
electrodes made of aluminum strips on a mylar 
substrate. 

The cylindrical RPCs have resistive electrodes 
made of a special plastic composed of a conduct- 
ing polymer and ABS plastic. The gap thickness 
and the spacers are identical to the planar RPCs. 
No linseed oil or any other surface treatments 
have been applied. The very thin and flexible 
electrodes are laminated to flberglass boards and 
foam to form a rigid structure. The copper read- 
out strips are attached to the fiberglass boards. 

10.3. RPC Design and Construction 

The IFR detectors cover a total active area of 
about 2,000 m^. There are a total of 806 RPC 
modules, 57 in each of the six barrel sectors, 108 
in each of the four half end doors, and 32 in the 
two cylindrical layers. The size and the shape of 
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Figure 73. Overview of the IFR: Barrel sectors and forward (FW) and backward (BW) end doors; the 
shape of the RFC modules and their dimensions are indicated. 



the modules are matched to the steel dimensions 
with very little dead space. More than 25 differ- 
ent shapes and sizes were built. Because the size 
of a module is limited by the maximum size of 
the material available, i.e., 320xl30cm^ for the 
bakelite sheets, two or three RFC modules are 
joined to form a gap-size chamber. The modules 
of each chamber are connected to the gas system 
in series, while the high voltage is supplied sepa- 
rately to each module. 

In the barrel sectors, the gaps between the steel 
plates extend 375 cm in the z direction and vary 
in width from 180 cm to 320 cm. Three modules 
are needed to cover the whole area of the gap, as 
shown in Figure |7|. Each barrel module has 32 
strips running perpendicular to the beam axis to 
measure the z coordinate and 96 strips in the or- 
thogonal direction extending over three modules 
to measure <j). 

Each of the four half end doors is divided into 
three sections by steel spacers that are needed 
for mechanical strength. Each of these sections is 
covered by two RFC modules that are joined to 
form a larger chamber with horizontal and verti- 



cal readout strips. 

The readout strips are separated from the 
ground aluminum plane by a 4 mm-thick foam 
sheet and form strip lines of 33 fi impedance. The 
strips are connected to the readout electronics at 
one end and terminated with a 2 kJ7 resistor at 
the other. Even and odd numbered strips are 
connected to different front-end cards (FECs), so 
that a failure of a card does not result in a to- 
tal loss of signal, since a particle crossing the gap 
typically generates signals in two or more adja- 
cent strips. 

The cylindrical RFC is divided into four sec- 
tions, each covering a quarter of the circumfer- 
ence. Each of these sections has four sets of two 
single gap RPCs with orthogonal readout strips, 
the inner with helical u-v strips that run paral- 
lel to the diagonals of the module, and the outer 
with strips parallel to and z. Within each sec- 
tion, the strips of the four sets of RFCs in a given 
readout plane are connected to form long strips 
extending over the whole chamber. Details of the 
segmentation and dimensions can be found in Ta- 
ble 
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Table 13 

IFR Readout segmentation. The total number of channels is close to 53,000. 





#of 




# of readout 


strips 


strip length 


strip width 


total # 


section 


sectors 


coordinate 


layers 


layer /sect 


(cm) 


(mm) 


channels 


barrel 


6 




19 


96 


350 


19.7-32.8 


« 11,000 






z 


19 


96 


190-318 


38.5 


« 11,000 


endcap 


4 


y 


18 


6x32 


124-262 


28.3 


13,824 






X 


18 


3x64 


10-180 


38.0 


« 15,000 


cylinder 


4 




1 


128 


370 


16.0 


512 






% 


1 


128 


211 


29.0 


512 






u 


1 


128 


10-422 


29.0 


512 






V 


1 


128 


10-423 


29.0 


512 



Prior to shipment to SLAC, all RPC modules 
were tested with cosmic rays. The single rates, 
dark currents, and efficiency were measured as a 
function of HV. In addition, detailed studies of 
the efficiency, spatial resolution, and strip multi- 
plicity were performed [ pTp^ ]. 

After the assembly of RPC modules into gap- 
size chambers, a new series of cosmic rays tests 
was performed to assure stable and efficient oper- 
ation. Before the installation of the steel flux re- 
turn, the planar chambers were inserted into the 
gaps. The cylindrical chambers were inserted af- 
ter the installation of the solenoid and the EMC. 

For each module, test results and conditions 
are retained in a database, together with records 
of the critical parameters of the components, the 
assembly and cabling. In addition, operational 
data are stored, such as the results of the weekly 
efficiency measurements that are used in the re- 
construction and simulation software. 

10.4. Power and Utilities 

Once the return flux assembly was completed, 
the FECs lH were installed and the low (LV) 
and high voltage (HV), and the gas system were 
connected. There are approximately 3,300 FECs, 
most placed inside the steel gaps, while the re- 
mainder was installed in custom crates mounted 
on the outside of the steel. 

Each FEC is individually connected to the LV 
power distribution. The total power required by 
the entire system is about 8kW at -1-7.0 V and 
2.5 kW at -5.2 V. The LV power is supphed by 



custom-built switching devices with load and line 
regulation of better than 1%. Additional features 
are precision shunts to measure output currents 
and TTL logic to inhibit output. 

The HV power system is a custom adaptation 
by CAEN (lOO). Each HV mainframe can hold 
up to ten cards, each carrying two independent 
10 kV outputs at 1 mA and 2 mA. The RPC mod- 
ules are connected via a distribution box to the 
HV supplies. Each distribution box services six 
RPC modules and up to six distribution boxes are 
daisy-chained to one HV output. Provisions are 
made for monitoring the currents drawn by each 
module. To reduce noise, the RPC ground plane 
is decoupled from the HV power supply ground 
by a 100 kJ7 resistor. 

The RPCs operate with a non-flammable gas 
mixture, typically 56.7% Argon, 38.8% Freon 
134a (1,1,1,2 tetrafluoroethane), and 4.5% isobu- 
tane. This mixture is drawn from a 760 liter 
tank that is maintained at an absolute pressure of 
1500-1600 Torr. The mixing tank is filled on de- 
mand with the three component gases under con- 
trol of mass-flow meters. Samples are extracted 
from the mixing tank periodically and analyzed 
to verify the correct mixture. 

The mixed gas is distributed at a gauge pres- 
sure of approximately 6.5 Torr through a parallel 
manifold system of 12.7 mm-diameter copper tub- 
ing. Each chamber is connected to the manifold 
through several meters of 6 mm-diameter plastic 
tubing (polyamide or Teflon). The flow to each of 
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Figure 75. Block diagram of the IFR electronics. 



these is adjusted individually with a small multi- 
turn metering valve. Protection against overpres- 
sure is provided by an oil bubbler to atmosphere 
in parallel with each chamber, limiting the gauge 
pressure in the chamber to a maximum of about 
1 Torr. Return flow of gas from each chamber 
is monitored by a second oil bubbler which cre- 
ates a back pressure of about 0.2 Torr. The total 
flow through the entire system is approximately 
5 -^/minute and corresponds on average to two gas 
exchanges per day. 

10.5. Electronics 

A block diagram of the IFR electronics system 



101 1 is shown in Figure 75, It includes the FECs, 
the data acquisition, and the trigger. 

The FECs service 16 channels each. They 
shape and discriminate the input signals and set 
a bit for each strip with a signal above a fixed 
threshold. The input stage operates continuously 
and is connected directly to the strips which act 
as transmission lines. A fast OR of all FEC in- 
put signals provides time information and is also 
used for diagnostic purposes. Two types of FECs 
are employed to handle inputs of different polar- 
ity for signals from the opposite sides of the gap. 
Because of the very low occupancy there is no 
provision for buffering during the trigger latency 



Signals from 3,300 FECs are transmitted to 
eight custom IFR front-end crates that are lo- 
cated near the detector. Each front-end crate 
houses up to 16 data handling cards, four trigger 
cards and a crate controller card (ICC) that col- 
lects data from the DAQ cards and forwards them 
to a ROM. There are three kinds of data cards: 
the FIFO boards (IFBs) that buffer strip hits, the 
TDC boards (ITBs) that provide time informa- 
tion, and the calibration boards (ICBs) that in- 
ject test pulses into the FECs. To deliver the data 
and clock signals to all the boards in the front-end 
crate, a custom backplane (PDB) for the standard 
6U Eurocard crate was designed using 9-layer 
strip line technology. Each board is connected 
to the ICC via three point-to-point lines for three 
single-end signals (data-in, data-out and clock), 
all of the same length and impedance (50 J7). 

The IFB reads the digital hit patterns from the 
FECs in less than 2.2 ms, stores the data into FI- 
FOs and transfers the FIFO contents into one of 
the ROMs. Each IFB handles 64 FECs, acting as 
an acquisition master. It receives commands via 
the PDB, and transmits and receives data pat- 
terns from the ROM (via C-Link and ICC). This 
card operates with the system clock frequency of 
59.5 MHz. 

The ICB is used for front-end tests and cali- 
brations. A signal with programmable amplitude 
and width is injected into the FEC input stage. 
To provide timing calibration and to determine 
the correct readout delay, the board is also used 
together with the TDCs. 

The ICC interfaces the crate backplane with 
the C-Link. The physical interface is the Finisar 
transceiver, a low cost and highly reliable data 
link for applications up to 1.5Cbytes/s. 

The TDC boards exploit the excellent time res- 
olution of the RPCs. Each board has 96 ECL 
differential input channels for the fast OR sig- 
nals from the FECs. Time digitization is achieved 
by three custom TDCs, designed at CERN pO^ . 
Upon receipt of a LI Accept, data are selected and 
stored untfl readout by the ROM. The 59.5 MHz 
clock signal is synchronized with the data and 
distributed to the 16 boards. High performance 
drivers provide a reliable clock distribution with 
a jitter of less than 0.5 ns. 
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10.6. Slow Controls 

and Online Monitoring 

The IFR is a system with a large number of 
components and electronics distributed all over 
the BABAR detector. To assure safe and stable 
operation, an extensive monitoring and control 
system was installed. The IFR Online Detector 
Control (lODC) monitors the performance of the 
RPCs by measuring the singles counting rate and 
the dark current of every module. It also controls 
and monitors the operation of the electronics, the 
DAQ and trigger, as well as the LV, the HV, and 
the gas system. The total number of hardware 
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channels is close to 2,500 [103 



The system has been easy to operate. HV trips 
are rare. Temperature monitoring in the steel 
structure and the electronics crates has proven 
very useful for the diagnosis of operational prob- 
lems. The occupancy is extremely low every- 
where, except in layer 18 of the forward end door 
which lacks adequate shielding from machine- 
generated background. On average, there are 
about 100-150 strip signals per event. 

10.7. Efficiency Measurements 
and Performance 

The efficiency of the RPCs is evaluated both for 
normal collision data and for cosmic ray muons 
recorded with the IFR trigger. Every week, cos- 
mic ray data are recorded at different voltage set- 
tings and the efficiency is measured chamber-by- 
chamber as a function of the applied voltage. The 
absolute efficiency at the nominal working voltage 
(typically 7.6 kV) is stored in the database for use 
in the event reconstruction software. 

To calculate the efficiency in a given chamber, 
nearby hits in a given layer and hits in different 
layers are combined to form clusters. Two differ- 
ent algorithms are used. The first is based solely 
on the IFR information and uses data recorded 
with a dedicated IFT trigger; the second matches 
the IFR clusters with the tracks reconstructed in 
the DCH. Both these algorithms start from one- 
dimensional IFR clusters defined as a group of ad- 
jacent hits in one of the two readout coordinates. 
The cluster position is defined as the ccntroid of 
the strips in the cluster. In the first algorithm, 
two-dimensional clusters are formed by joining 
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Figure 76. Distribution of the efficiency for all 
RFC modules measured with cosmic rays in June 
1999. Some 50 modules were not operational at 
that time. 

one-dimensional clusters (of the same readout 
coordinate) in different layers. In each sector, 
two-dimensional clusters in different coordinates 
are combined into three-dimensional clusters pro- 
vided there are fewer than three layers missing in 
one of the two coordinates. The second algorithm 
extrapolates charged tracks reconstructed by the 
DCH. IFR clusters which are less than 12 cm 
from the extrapolated track are combined to form 
three-dimensional or two-dimensional clusters. A 
detailed discussion o f the clustering algorithm can 
be found elsewhere |104]. 

The residual distributions from straight line fits 
to two-dimensional clusters typically have an rms 
width of less than 1 cm. An RFC is considered 
efficient if a signal is detected at a distance of less 
than 10 cm from the fitted straight line in either 
of the two readout planes. Following the instal- 
lation and commissioning of the IFR system, all 
RFC modules were tested with cosmic rays and 
their efficiency was measured. The results are 
presented in Figure ^ Of the active RFC mod- 
ules, 75% exceed an efficiency of 90%. 

Early tests indicated that the RFC dark cur- 
rent was very temperature dependent, specifi- 
cally, the current increases 14-20% per °C. Be- 
cause the IR experimental hall does not have tem- 
perature regulation this presents a serious prob- 
lem. The FECs that are installed in the steel 
gaps dissipate 3 W each, generating a total power 
of 3.3 kW in the barrel and 1.3 kW in the forward 
end door. 
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Figure 77. History of the temperature and dark 
current in the RPC modules since January 2000. 
a) temperature in the IR-2 haU and in the back- 
ward end door; b) total dark current in the 216 
modules of the backward end door. 



During the first summer of operation, the daily 
average temperature in the IR hall was 28°C 
and the maximum hall temperature frequently 
exceeded 31°C. The temperature inside the steel 
rose to more than 37° C and the dark currents 
in many modules exceeded the capabilities of the 
HV system and some RPCs had to be temporarily 
disconnected. 

To overcome this problem, water cooling was 
installed on the barrel and end door steel, remov- 
ing salOkW of heat and stabilizing the tempera- 
ture at 20-21°C in the barrel, 22°C in the back- 
ward and 24° C in the forward end doors. Fig- 
ure ^ shows the history of temperature in the 
hall and temperature and total dark current in 
the backward end door. While the current closely 
follows the temperature variations, the range of 
change is now limited to a few degrees. 

During operation at high temperatures, a large 
fraction of the RPCs (>50%) showed not only 
very high dark currents, but also some reduc- 
tion in efficiency compared to earlier measure- 
ments [105|. After the cooling was installed and 
the RPCs were reconnected, some of them contin- 
ued to deteriorate while others remained stable, 
some of them (> 30%) at full efficiency, (see Fig- 
ure ^). Detailed studies revealed large regions of 
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Figure 78. Efficiency history for 12 months start- 
ing in June 1999 for RPC modules showing dif- 
ferent performance: a) highly efficient and stable; 
b) continuous slow decrease in efficiency; c) more 
recent, faster decrease in efficiency. 

very low efficiency in these modules, but no clear 
pattern was identified. 

The cause of the efficiency loss remains under 
investigation. Several possible causes have been 
excluded as the primary source of the problem, 
such as a change in the bakelite bulk resistivity, 
loosened spacers, gas flow, or gas composition. 
A number of prototype RPCs developed similar 
efficiency problems after being operated above a 
temperature 36°C for a period of two weeks. In 
some of these modules, evidence was found that 
the linseed oil had failed to cure and had accumu- 
lated at various spots under the influence of the 
electric field. 

10.8. Muon Identification 

While muon identification relies almost entirely 
on the IFR, other detector systems provide com- 
plementary information. Charged particles are 
reconstructed in the SVT and DCH and muon 
candidates are required to meet the criteria for 
minimum ionizing particles in the EMC. Charged 
tracks that are reconstructed in the tracking sys- 
tems are extrapolated to the IFR taking into ac- 
count the non-uniform magnetic field, multiple 



81 



scattering, and the average energy loss. The pro- 
jected intersections with the RPC planes are com- 
puted and for each readout plane all clusters de- 
tected within a predefined distance from the pre- 
dicted intersection are associated with the track. 

A number of variables are defined for each IFR 
cluster associated with a charged track to discrim- 
inate muons from charged hadrons: 1) the to- 
tal number of interaction lengths traversed from 
the IP to the last RPC layer with an associated 
cluster, 2) the difference between this measured 
number of interaction lengths and the number of 
interaction lengths predicted for a muon of the 
same momentum and angle, 3) the average num- 
ber and the rms of the distribution of RPC strips 
per layer, 4) the for the geometric match be- 
tween the projected track and the centroids of 
clusters in different RPC layers, and 5) the of 
a polynomial fit to the two-dimensional IFR clus- 
ters. Selection criteria based on these variables 
are applied to identify muons. 

The performance of muon identification has 
been tested on samples of muons from /i/zee and 
/i/i7 final states and pions from three-prong r de- 
cays and Kg tt+tt" decays. The selection of 
these control samples is based on kinematic vari- 
ables, and not on variables used for muon identifi- 
cation. As illustrated in Figure |79|, a muon detec- 
tion efficiency of close to 90% has been achieved 
in the momentum range of 1.5 < p < 3.0GeV/c 
with a fake rate for pions of about 6-8%. Decays 
in flight contribute about 2% to the pion misiden- 
tification probability. The hadron misidcntifica- 
tion can be reduced by a factor of about two by 
tighter selection criteria which lower the muon 
detection efficiency to about 80%. 

10.9. K1 and Neutral Hadron Detection 

KI^s and other neutral hadrons interact in the 
steel of the IFR and can be identified as clus- 
ters that are not associated with a charged track. 
Monte Carlo simulations predict that about 64% 
of K^s above a momentum of 1 GeV/c produce a 
cluster in the cylindrical RPC, or a cluster with 
hits in two or more planar RPC layers. 

Unassociated clusters that have an angular sep- 
aration of < 0.3 rad are combined into a com- 
posite cluster, joining clusters that originate from 



1.0 



0.8 



o 0.6 
92 

£ 0.4 



0.2 



0.0 



1 

" a) 


1 1 


• 
• 




• 




• 
• 

• 


± 
± 

□ 7T 


• 












#1 1 1 1 



0.5 



0.4 



- 0.3 



0.2 



1 2 
Momentum (GeV/c) 



0.1 



0.0 




3-2001 
8583A44 



40 80 1 20 
Polar Angle (degrees) 



160 



Figure 79. Muon efficiency (left scale) and pion 
misidcntification probability (right scale) as a 
function of a) the laboratory track momentum, 
and b) the polar angle (for 1.5 < p < 3.0GeV/c 
momentum), obtained with loose selection crite- 
ria. 



showers that spread into adjacent sectors of the 
barrel, several sections of the end doors and/or 
the cylindrical RPC. This procedure also com- 
bines multiple clusters from large fiuctuations in 
the hadronic showers. The direction of the neu- 
tral hadron is determined from the event vertex 
and the centroid of the neutral cluster. No infor- 
mation on the energy of the cluster can be ob- 
tained. 
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Figure 80. Angular difference, cosA^, between 
the direction of the missing momentum and the 
closest neutral IFR cluster for a sample of (f) 
mesons produced in the reaction e^e~ 4>"/ with 

Since a significant fraction of hadrons interact 
before reaching the IFR, information from the 
EMC and the cylindrical RPCs is combined with 
the IFR cluster information. Neutral showers in 
the EMC are associated with the neutral hadrons 
detected in the IFR, based on a match in produc- 
tion angles. For a good match, a probability 
of > 1% is required. 

An estimate of the angular resolution of the 
neutral hadron cluster can be derived from a sam- 
ple of K'^'s produced in the reaction e+e^ 
07 — > K'^K'^^. The K'l direction is inferred from 
the missing momentum computed from the mea- 
sured particles in the final state. The data in 
Figure ^ indicate that the angular resolution of 
the derived from the IFR cluster information 
is of the order of 60mrad. For K^'s interacting 
in the EMC, the resolution is better by about a 
factor of two. 

For multi-hadron events with a reconstructed 
J/ip decay. Figure |8l| shows the angular differ- 
ence, between the missing momentum and 
the direction of the nearest neutral hadron clus- 
ter. The observed peak demonstrates clearly that 
the missing momentum can be associated with a 
neutral hadron, assumed to be a K^. The de- 
tection efficiency increases roughly linearly with 
momentum; it varies between 20% and 40% in the 
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Figure 81. Difference between the direction of 
the reconstructed neutral hadron cluster and the 
missing transverse momentum in events with a 
reconstructed J/ip decay. The Monte Carlo simu- 
lation is normalized to the luminosity of the data; 
the background is obtained using neutral hadrons 
and the missing momentum from different events. 

momentum range from 1 GeV/c to 4 GeV/ c (EMC 
and IFR combined). 

10.10. Summary and Outlook 

The IFR is the largest RFC system built to 
date. It provides efficient muon identification and 
allows for the detection of if °'s interacting in the 
steel and the calorimeter. During the first year of 
operation, a large fraction of the RFC modules 
have suffered significant losses in efficiency. This 
effect appears to be correlated with high temper- 
atures, but the full extent of the problem and its 
cause remain under study. Thanks to the large 
number of RFC layers, this problem has not yet 
impacted the overall performance severely. But 
present extrapolations, even after installation of 
water cooling on the steel, indicate a severe prob- 
lem for the future operation. Recently, 24 end 
door modules have been replaced by new RPCs 
with a substantially thinner coating of linseed oil 
and improved treatment of the bakelite surfaces. 
Results with these new RPCs and other tests will 
need to be evaluated before decisions on future 
improvements of the IFR can be made. Further- 
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more, it is planned to reduce the contamination 
from hadron decays and punch through by in- 
creasing the absorber thickness, i.e., adding more 
steel on the outside and replacing a few of the 
RPCs with absorber plates. 

11. Trigger 

11.1. Trigger Requirements 

The basic requirement for the trigger system is 
the selection of events of interest (see Table 14) 
with a high, stable, and well-understood efficiency 
while rejecting background events and keeping 
the total event rate under 120 Hz. At design lumi- 
nosity, beam-induced background rates are typi- 
cally about 20 kHz each for one or more tracks 
in the drift chamber with pt > 120MeV/c or at 
least one EMC cluster with E > 100 MeV. Ef- 
ficiency, diagnostic, and background studies re- 
quire prescaled samples of special event types, 
such as those failing the trigger selection criteria, 
and random beam crossings. 

The total trigger efficiency is required to exceed 
99% for all BB events and at least 95% for contin- 
uum events. Less stringent requirements apply to 
other event types, e.g., t^t^ events should have 
a 90-95% trigger efficiency, depending on the spe- 
cific decay channels. 

The trigger system must be robust and flexible 
in order to function even under extreme back- 
ground situations. It must also be able to operate 
in an environment with dead or noisy electronics 
channels. The trigger should contribute no more 
than 1% to dead time. 

11.2. Trigger Overview 

The trigger is implemented as a two-level hier- 
archy, the Level 1 (LI) in hardware followed by 
the Level 3 (L3) in software. It is designed to ac- 
commodate up to ten times the initially projected 
m PEP-II background rates at design luminos- 
ity and to degrade slowly for backgrounds above 
that level. Redundancy is built into the system 
to measure and monitor trigger efficiencies. 

During normal operation, the LI is configured 
to have an output rate of typically 1 kHz. Triggers 
are produced within a fixed latency window of 11- 
12 fj,s after the e+e^collision, and delivered to the 



Fast Control and Timing System (FCTS). Data 
used to form the trigger decision are preserved 
with each event for efficiency studies. 

The L3 receives the output from LI, performs a 
second stage rate reduction for the main physics 
sources, and identifies and flags the special cat- 
egories of events needed for luminosity determi- 
nation, diagnostic, and calibration purposes. At 
design luminosity, the L3 filter acceptance for 
physics is ^90 Hz, while ^30 Hz contain the other 
special event categories. The L3 algorithms com- 
ply with the same software conventions and stan- 
dards used in all other BABAR software, thereby 
simplifying its design, testing, and maintenance. 

11.3. Level 1 Trigger System 

The LI trigger decision is based on charged 
tracks in the DCH above a preset transverse mo- 
mentum, showers in the EMC, and tracks de- 
tected in the IFR. Trigger data are processed 
by three specialized hardware processors. As de- 
scribed below, the drift chamber trigger (DCT) 
and electromagnetic calorimeter trigger (EMT) 
both satisfy all trigger requirements indepen- 
dently with high efficiency, and thereby provide 
a high degree of redundancy, which enables the 
measurement of trigger efficiency. The instru- 
mented flux return trigger (IFT) is used for trig- 
gering fJ.~^ fJ.~ and cosmic rays, mostly for diagnos- 
tic purposes. 

The overall structure of the LI trigger is illus- 

Table 14 

Cross sections, production and trigger rates for 
the principal physics processes at 10.58 GeV for a 
luminosity of 3 x 10^'^ cm~^s~^. The e+e^ cross 
section refers to events with either the e^, e~, or 
both inside the EMC detection volume. 
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Figure 82. Simplified LI trigger schematic. Indi- 
cated on the figure are the number of components 
(in square brackets), and the transmission rates 
between components in terms of total signal bits. 

trated in Figure |8^. Each of the three LI trigger 
processors generates trigger primitives, summary 
data on the position and energy of particles, that 
are sent to the global trigger (GLT) every 134 ns. 
The DCT and EMT primitives sent to the GLT 
are (p-maps. An individual ^-map consists of an 
n-bit word representing a particular pattern of 
trigger objects as distributed in fixed-width (j) re- 
gions from to 27r. A trigger object is a quantity 
indicating the presence of a particle, such as a 
drift chamber track or a calorimeter energy de- 
posit. The IFT primitive is a three-bit pattern 
representing the hit topology in the IFR. The 
meaning of the various trigger primitive inputs 



to the GLT are summarized in Table 15 



The GLT processes all trigger primitives to 
form specific triggers and then delivers them to 
the FCTS. The FCTS can optionally mask or 
prescale any of these triggers. If a valid trigger 
remains, a LI Accept is issued to initiate event 
readout. The trigger definition logic, masks, and 
prescale values are all configurable on a per run 



basis. 

The LI hardware is housed in five 9U VME 
crates. The LI trigger operates in a continuous 
sampling mode, generating trigger information at 
regular, fixed time intervals. The DCH front-end 
electronics (FEEs) and the EMC untriggered per- 
sonality cards (UPCs) send raw data to the DCT 
and EMT about 2 /is after the e'^e'~ collision. 
The DCT and EMT event processing times are 
4-5 /iS, followed by another ~3 /iS of processing 
in the GLT to issue a LI trigger. The LI trigger 
takes approximately 1 /is to propagate through 
the FCTS and the readout modules (ROMs) to 
initiate event readout. These steps are all ac- 
complished within the 12.8 /zs FEE buffer capac- 
ity limit. 

The DCT, EMT and GLT each maintain a four- 
event buffer to hold information resulting from 
the various stages of the LI trigger. These data 
are read out by the normal data acquisition sys- 
tem. 



11.3.1. Level 1 Drift Chamber Trigger 

The input data to the DCT consist of one bit 
for each of the 7104 DCH cells. These bits convey 
time information derived from the sense wire sig- 
nal for that cell. The DCT output primitives are 
candidate tracks encoded in terms of three 16-bit 
0-maps as listed in Table |5[ 

The DCT algor ithm s are executed in three 
types of modules |106|. First, track segments, 
their <j) positions and drift time estimates are 
found using a set of 24 Track Segment Finder 
(TSF) modules |T07| . These data are then 
passed to the Binary Link Tracker (BLT) mod- 
ule ]10^ , where segments are linked into complete 
tracks. In parallel, the (p information for segments 
found in axial superlayers is transmitted to eight 
transverse momentum discriminator (PTD) mod- 
ules |10£], which search for tracks above a set pt 
threshold. 

Each of the three DCT modules (TSF, BLT, 
and PTD) rches heavily on multiple FPGA's fll^ 
which perform the control and algorithmic func- 
tions. All cabling is handled by a small (6U) back- 
of-crate interface behind each main board. 
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Table 15 

Trigger primitives for the DCT and EMT. Most energy thresholds are adjustable; those listed are typical 
values. 
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Track Segment Finder 

The TSF modules are responsible for find- 
ing track segments in 1776 overlapping eight-cell 
pivot groups. A pivot group is a contiguous set 
of cells that span all four layers within a super- 
layer. The pivot group shape is such that only 
reasonably straight tracks originating from the 
interaction point can produce a valid segment. 
Figure ^ shows the arrangement of cells within 
a pivot group. Cell 4 is called the pivot cell; 
the TSF algorithm is optimized to find track seg- 
ments that pivot about this cell. 

The DCH signals are sampled every 269 ns. 
The passage of a single particle through the DCH 
will produce ionization that drifts to the sense 
wires in typically no more than four of these 
clock ticks. Each cell is associated with a two-bit 
counter that is incremented at every clock tick for 
which a signal is present. In this way, a short time 
history of each cell is preserved. For each clock 
tick, the collection of two-bit counters for each 
pivot group forms a 16-bit value used to address 
a look-up-table. This look-up-table contains two- 
bit weights indicating whether there is no accept- 
able segment, a low-quality segment, a three-layer 
segment (allowing for cell inefRciencies) , or a four- 
layer segment. When an acceptable segment is 
found, that pivot group is examined to determine 
which of three subsequent clock ticks produce the 
highest weight or best pattern. 

The look-up-table also contains position and 
time information which, along with a summary of 
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Figure 83. Track Segment Finder pivot group. 

cell occupancies, forms the basis of data sent to 
the BLT and PTD. The TSF algorithm uses the 
time- variation of the look-up-table weights to re- 
fine both the event time and its uncertainty, thus 
enabling it to output results to the BLT every 
134 ns. 

The position resolution as measured from the 
data after calibration, is typically ~600 fim for a 
four-layer segment and ~900 /xm for a three-layer 
segment. For tracks originating from the IP, the 
efficiency for finding TSF segments is 97%, and 
the efficiency for high-quality three-layer or four- 
layer TSF segments is 94%. 

Binary Link Tracker 

The BLT receives segment hit information from 
all 24 TSF's at a rate of 320 bits every 134 ns 
and links them into complete tracks. The seg- 
ment hits are mapped onto the DCH geometry in 
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terms of 320 supercells, 32 sectors in and ten 
radial superlayers. Each bit indicates whether a 
segment is found in that superceU or not. The 
BLT input data are combined using a logical OR 
with a programmable mask pattern. The mask- 
ing allows the system to activate track segments 
corresponding to dead or highly inefficient cells 
to prevent efficiency degradation. The linking al- 
gorithm uses an extensio n of a method developed 
for the CLEO-11 trigger [pTT| . It starts from the 
innermost superlayer, Al, and moves radially out- 
ward. 

Tracks that reach the outer layer of the DCH 
(superlayer AIO) are classified as type A. Tracks 
that reach the middle layer (superlayer U5) are 
classified as type B. An A track is found if there 
is a segment in at least eight superlayers and if 
the segments in two consecutive superlayers fall 
azimuthally within three to five supercells of each 
other (depending on the superlayer type). This 
allows for track curvature and dip angle varia- 
tions. The data are compressed and output to 
the GLT in the form of two 16-bit 0-maps, one 
each for A and B tracks. 

PT Discriminator 

The eight PTD modules receive (f) information 
of high quality track segments in the axial super- 
layers (Al, A4, A7 and AlO), and determine if the 
segments are consistent with a track pt greater 
than a configurable minimum value. An envelope 
for tracks above the minimum pt is defined using 
the IP, and a track segment position in one of the 
seed superlayers, A7 or AIO. A high pt candidate, 
denoted as A', is identified when sufficient track 
segments with accurate cf) information from the 
other axial superlayers lie within this envelope. 

Each PTD module searches for seed segments 
in superlayers A7 and AIO, and within a 45- 
degree azimuthal wedge of the DCH. This search 
region spans eight supercells, and the processing 
for each supercell is performed by its own process- 
ing engine on the PTD. The principal components 
in each engine are an algorithmic processor and 
look-up-tables containing the limits for each in- 
dividual seed position. The contents of the look- 
up-tables specify the allowed track segment posi- 
tions for each of the three other axial superlayers 
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Figure 84. DCT track efficiency versus trans- 
verse momentum for A, B, and A' tracks. The 
A' threshold is set to 800 MeV/c. 

and consequently define the effective pt discrim- 
ination threshold. The resulting pt threshold for 
the PTD A' tracks is shown in Figure |8j together 
with the BLT A, B track efficiency. 

11.3.2. Level 1 Calorimeter Trigger 

For trigger purposes, the EMC is divided into 
280 towers, 7 x 40 (6* x 0). Each of the barrel's 
240 towers is composed of 24 crystals in a 8 x 
3 (6* X 0) array. The endcap is divided into 40 
towers, each forming a wedge in cj) containing 19- 
22 crystals. For each tower, all crystal energies 
above a threshold of 20 MeV are summed and sent 
to the EMT every 269 ns. 

The conversion of the tower data into the GLT 
(/)-maps is performed by ten Trigger Processor 
Boards (TPBs). The TPBs determine energies 
in the 40 (j) sectors, summing over various ranges 
of 6, compare these energies against thresholds 
for each of the trigger primitives (see Table p^), 
estimate the time of energy deposition, correct for 
timing jitter, and then transmit the result to the 
GLT. 

Each TPB receives data from 28 towers, corre- 
sponding to an array of 7 x 4 in x 0, or four 0- 
sectors. Each of the 40 0-sectors is summed inde- 
pendently. To identify energy deposits that span 
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two adjacent ^-sectors, the energy of each sector 
is also made available to the summing circuit for 
a single adjoining sector in such a way that all 
possible pairs of adjacent 0-sectors are summed. 
These energy sums are compared against thresh- 
olds to form trigger objects. Each sum is also sent 
to an eight-tap finite impulse response (FIR) dig- 
ital filter which is used to estimate the energy 
deposition time. A look-up-table is used to make 
an energy-dependent estimate of the timing jit- 
ter which, along with the FIR output, is used to 
time the transmission of any trigger objects to the 
GLT. Pairs of 0-sectors are ORed to form 20-bit 
(^-maps for the M, G, E, and X primitives, while 
for the Y primitive, groups of four are ORed to 
form a 10-bit 0-map. The complete algorithm is 
implemented in one FPGA |112| | for each 0-sector, 
with four identical components per TPB. Further 
details of the EMT system can be found in [113|. 

The basic performance of the EMT can be ex- 
pressed in terms of the efficiency and timing jit- 
ter of the trigger primitives. The efficiency of 
the primitives can be measured by the number of 
times a trigger bit is set for a specific energy re- 
constructed offline in events from a random trig- 
ger. Figure ^ shows this efficiency for energies 
near the M threshold. The efficiency changes 
from 10% to 90% in the range of 110 to 145 MeV, 
and reaches 99% at 180 MeV, close to the average 
energy deposition of a minimum ionizing particle 
at normal incidence. 

The EMT time jitter is measured by comparing 
the time centroid of (/)-strip M hits in events 
with the DCH track start time, to. The difference 
has an rms width of 90 ns with >99.9% of the 
matching M hits within a ±500 ns window. 



11.3.3. Level 1 IFR Trigger 

The IFT is used for triggering on /i^/i^ and 
cosmic rays. For the purposes of the trigger, the 
IFR is divided into ten sectors, namely the six 
barrel sextants and the four half end doors. The 
inputs to the IFT are the Fast OR signals of all 
(j) readout strips in eight selected layers in each 
sector. 

A majority logic algorithm defines trigger ob- 
jects for every sector in which at least four of the 
eight trigger layers have hits within a time win- 
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Figure 85. EMT M efficiency vs. EMC cluster 
energy for an M threshold setting of 120 MeV. 



Table 16 

IFR trigger pattern (U) definition, where fi refers 
to a signal within a sector. 

U Trigger condition 

1 > 2/^ topologies other than U = 5 — 7 

2 1 /i in backward endcap 

3 1 /i in forward endcap 

4 1 /X in barrel 

5 2 back-back /i's in barrel -f 1 forward fi 

6 1 /X in barrel -1-1 forward ^ 

7 2 back-back /z's in barrel 

dow of 134 ns. The IFR trigger synchronization 
module processes the trigger objects from the ten 
sectors and generates the three-bit trigger word 
(U) encoding seven exclusive trigger conditions, 
as defined in Table |l^. The trigger U > 5, for 
example, covers all /.t^/i~ topologies of interest. 

The efficiency of the IFT has been evaluated us- 
ing cosmic rays triggered by the DCT and cross- 
ing the detector close to the IP. For these events, 
98% were triggered by the IFT as events with 
at least one track, and 73% as events with two 
tracks, inside the geometrical region of the IFR. 
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Most of the IFT inefficiency is concentrated at 
the boundaries between sectors. 



plemented as an array of f6 memory chips with 
8 Mbytes of configuration data. 



11.3.4. Global Trigger 

The GLT receives the eight trigger primitives 
in the form of i/i-maps as hsted in TableJTs] along 
with information from the IFT (Table |16|) to form 
specific triggers that are then passed to the FCTS 
for the final trigger decision. Due to the different 
latencies associated with the production of these 
primitives, the GLT forms a time alignment of 
these input data using configurable delays. 

The GLT then forms some additional combined 
0-maps from the DCT and EMT data. These 
maps include matched objects such as BM for B 
tracks matched to an M cluster in 0, back-to-back 
objects, B* and M*, which require a pair of (p 
bits separated by a configurable angle of typically 
- 120°, and an EM* object for back-to-back EM 
pairs. 

All 16 (p-maps are then used to address indi- 
vidual GLT look-up-tables which return three-bit 
counts of trigger objects contained within those 
maps, e.g., the number of B tracks or number of 
M clusters. To count as distinct trigger objects, 
the map bits are typically required to have a sep- 
aration of more than one (j) bin. The resulting 16 
counts plus the IFT hit pattern are then tested in 
logical operations. The permissible operations in- 
clude: always-pass; or a comparison (>, —, or <) 
with a configurable selection parameter. A trig- 
ger line is then set as the logical AND of these 17 
operations. This process is performed for each of 
the 24 trigger lines. 

The GLT derives the LI trigger time from the 
centroid of the timing distribution of the highest 
priority trigger, binned in the 134 ns interval and 
spanning about 1 /is. Other trigger signals com- 
patible with this time are retained and cached. 
The average time is calculated to the nearest 67 ns 
and the 24-bit GLT output signal is sent to the 
FCTS every 67 ns. The achieved timing resolu- 
tion for hadronic events has an rms width of 52 ns; 
and 99% of the events are within 77 ns. 

The GLT hardware consists of a single 9U 
VME module. Most of the logic, including di- 
agnostic and DAQ memories, are implemented in 
FPGA's [[llOt . The look-up-table section is im- 



11.4. Level 1 Trigger Performance 
and Operational Experience 

The LI trigger configuration consists of DCT- 
only, EMT-only, mixed and prescaled triggers, 
aimed not only for maximum efficiency and back- 
ground suppression, but also for the convenience 
of trigger efficiency determination. 

Although most triggers target a specific physics 
source, they often also select other processes. For 
example, two-track triggers are not only efficient 
for Bhabha, and t^t~ events, but are also 

useful for selecting jet-like hadronic events and 
some rare B decays. 

The efficiencies and rates of selected LI trig- 
gers for various physics processes are listed in 
Table ^ Although triggering on generic BB 
events is relatively easy, it is essential to en- 
sure high efficiencies for the important rare low- 
multiplicity B decays. For this reason, efficiencies 
for B^ 7r°7r° and B^ t^V are also listed in 
Table 0. 

The efficiencies listed for the hadronic events 
are absolute and include acceptance losses based 
on Monte Carlo simulation, and local inefficiency 
effects. The efficiencies for r-pair events are for 
fiducial events, i.e., events with two or more 
tracks with pt > 120MeV/c and polar angle 9 to 
reach at least DCH superlayer U5. The Bhabha 
and /i-pair efficiencies are determined from the 
data, for events with two high momentum parti- 
cles, which are back-to-back in the cm. system, 
and within the EMC fiducial volume. The data in 
Table ^ demonstrate that the DCT and the com- 
bined EMT/IFT provide fully efficient, indepen- 
dent triggers for most physics processes, although 
independent triggers for fJ-^ fJ." and t~^t~ are not 
individually fully efficient. The efficiencies pre- 
dicted by the Monte Carlo simulation are gener- 
ally in good agreement with data when tested us- 
ing events passing typical analysis selections and 
based on orthogonal triggers. Prescaled triggers 
with a very open acceptance of physics events, 
such as (B>2 & A>1) or (M>2) are also used to 
measure the trigger efficiencies. 



The trigger rates listed in Table 17 are for a 
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Table 17 

Level 1 Trigger efficiencies (%) and rates (Hz) at a luminosity of 2.2 x 10^'^ cm~^s~"'^ for selected triggers 
applied to various physics processes. The symbols refer to the counts for each object. 
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Rate 


A>3 & B* >1 


97.1 


66.4 


81.8 


88.9 


81.1 








17.7 


180 


A>1 & B* >1 & A' >1 


95.0 


63.0 


83.2 


89.2 


85.2 


98.6 


99.1 


79.9 


410 


Combined DCT (ORed) 


99.1 


79.7 


92.2 


95.3 


90.6 


98.9 


99.1 


80.6 


560 


M>'? ^ M* >1 


99.7 


98.6 


93.7 


98.5 


94.7 






53.7 


160 


EM* >1 


71.4 


94.9 


55.5 


77.1 


79.5 


97.8 




65.8 


150 


Combined EMT (ORed) 


99.8 


99.2 


95.5 


98.8 


95.6 


99.2 




77.6 


340 


B>3 & A>2 & M>2 


99.4 


81.2 


90.3 


94.8 


87.8 






19.7 


170 


M* >1 & A>1 & A' >1 


95.1 


68.8 


83.7 


90.1 


87.0 


97.8 


95.9 


78.2 


250 


E>1 & B>2 & A>1 


72.1 


92.4 


60.2 


77.7 


79.2 


99.3 




72.8 


140 


M* >1 & U>5 (M-pair) 














60.3 




70 


Combined Level 1 triggers 


>99.9 


99.8 


99.7 


99.9 


98.2 


>99.9 


99.6 


94.5 


970 



typical run with HER (LER) currents at 650 mA 
(1350 mA) and a luminosity of 2.2 x 10^^ cm^^s"^ 
These rates are stable to within 20% for the same 
PEP-II configuration, but they are impacted by 
changes in vacuum conditions, beam currents, 
and orbits. There are occasional background 
spikes which can double the LI rate. However, 
due to the 2 kHz capability of the data acquisi- 
tion, these spikes do not induce significant dead 
time. 

For a typical LI rate of IkHz, Bhabha and 
annihilation physics events contribute ~130Hz. 
There are also 100 Hz of cosmic ray and 20 Hz of 
random beam crossing triggers. The remaining 
triggers are due to lost particles interacting with 
the beam pipe or other components. The distri- 
bution of single track zq values as reconstructed 
by L3 for all LI triggers is shown in Figure ^ 
The most prominent peaks aX z — ±20 cm cor- 
respond to a flange of the beam pipe. The peak 
at zq = —55 cm corresponds to a step in the syn- 
chrotron mask. 

The LI trigger hardware operation has been 
very stable. For the first one and half years of op- 
eration, there have been only four hardware fail- 
ures in the LI system, mainly auxiliary or com- 
munication boards. Occasional adjustments to 
the EMT tower mask were used to temporarily 



■p, 8000 - 



4-H 

o 

d 6000 - 

2: 



4000 - 



2000 - 




-80 -40 40 80 

L3 Track z„ (cm) 



Figure 86. Single track zq for all LI tracks, re- 
constructed by L3. 

suppress noisy channels in the EMC electronics. 

11.5. Level 3 Trigger System 

The L3 trigger software comprises event recon- 
struction and classification, a set of event selec- 
tion filters, and monitoring. This software runs 
on the online computer farm. The filters have ac- 
cess to the complete event data for making their 
decision, including the output of the LI trigger 
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processors and FCTS trigger scalers. L3 operates 
by refining and augmenting the selection meth- 
ods used in LI. For example, better DCH track- 
ing (vertex resolution) and EMC clustering filters 
allow for greater rejection of beam backgrounds 
and Bhabha events. 

The L3 system runs within the Online Event 
Processing (OEP) framework (see Section [l^ ). 
OEP delivers events to L3, then prescales and 
logs those which pass the L3 selection criteria. 

To provide optimum flexibility under different 
running conditions, L3 is designed according to a 
general logic model that can be configured to sup- 
port an unlimited variety of event selection mech- 
anisms. This provides for a number of different, 
independent classification tests, called scripts, 
that are executed independently, together with 
a mechanism for combining these tests into the 
final set of classification decisions. 

The L3 trigger has three phases. In the first 
phase, events are classified by defining L3 input 
lines, which are based on a logical OR of any num- 
ber of the 32 FCTS output lines. Any number of 
L3 input lines may be defined. 

The second phase comprises a number of 
scripts. Each script executes if its single L3 input 
line is true and subsequently produces a single 
pass-fail output fiag. Internally, a script may ex- 
ecute one or both of the DCH or EMC algorithms, 
followed by one or more filters. The algorithms 
construct quantities of interest, while the filters 
determine whether or not those quantities satisfy 
the specific selection criteria. 

In the final phase, the L3 output lines are 
formed. Each output line is defined as the log- 
ical OR of selected script flags. L3 can treat 
script flags as vetoes, thereby rejecting, for exam- 
ple, carefully selected Bhabha events which might 
otherwise satisfy the selection criteria. 

L3 utilizes the standard event data analysis 
framework and depends crucially on several of its 
aspects. Any code in the form of modules can be 
included and configured at run time. A sequence 
of these software modules compose a script. The 
same instance of a module may be included in 
multiple scripts yet it is executed only once, thus 
avoiding significant additional CPU overhead. 



11.5.1. Level 3 Drift Chamber 
Tracking Algorithm 

Many events which pass LI but must be re- 
jected by L3 are beam-induced charged particle 
background that are produced in material close 
to the IP. LI does not currently have sufficient 
tracking resolution to identify these background 
tracks. The DCH-based algorithm, L3Dch, per- 
forms fast pattern recognition (track finding) and 
track fitting, which determines the five helix track 
parameters for tracks with pt above 250MeV/c. 
To speed up the process of pattern recognition, 
L3Dch starts with the track segments from the 
TSF system and improves the resolution by mak- 
ing use of the actual DCH information. 

For those TSF segments that have a simple so- 
lution to the left-right ambiguity, a track to is 
determined. The to values for each segment in an 
event are binned and the mean produced from the 
values in the most populated bin is used as the 
estimated event tg. All events which pass LI typi- 
cally have enough segments to form a Iq estimate. 
The measured rms resolution on this estimate is 
1.8 ns for Bhabha events and 3.8 ns for hadronic 
events. 

The pattern recognition for L3Dch is done with 
a look-up-table. For this track table, the DCH is 
divided into 120 0-sectors, corresponding to the 
number of cells in the innermost layers. The track 
table is populated with the hit patterns of Monte 
Carlo generated tracks with a pt above 250 MeV/c 
and originating within 2 cm of the IP in the x-y 
plane, and within 10 cm in z. The pattern recog- 
nition algorithm searches the table entries looking 
for matches to segments found by the TSFs. The 
matched set of segments for a given track is then 
passed to the track fitting algorithm. The track 
table allows for up to two missing DCH TSF seg- 
ments per track. 

The track fitting algorithm is provided with 
both the track segments found in pattern recog- 
nition and the individual hits within those seg- 
ments. From this information the five helix pa- 
rameters are fitted. The fit is then iterated, 
adding segments close to the initially fitted track, 
and dropping hits with large residuals. The final 
fit does not demand that the track originate from 
the IP. 
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Figure 87. Transverse and longitudinal miss dis- 
tances between the two tracks in Bhabha events. 

The two-track miss distances for Bhabha events 
are plotted in Figure [s^. The resolutions for indi- 
vidual tracks are 0.80 mm and 6.1 mm for do and 
zq, respectively. Similarly, the 1/pt difference be- 
tween the two tracks in /i-pair events yields a pt 
resolution of (5pt /pt ^ 0.019 -pt, withpt in GeV/c. 

11.5.2. Level 3 Calorimeter 
Clustering Algorithm 

The all-neutral trigger for L3 is based on infor- 
mation from the EMC. In addition, calorimeter 
information is a vital complement to the DCH 
data for the identification of Bhabha events. 

The L3 EMC-based trigger, L3Emc, identifies 
energy clusters with a sensitivity sufficient for 
finding minimum ionizing particles. EMC data 
are processed in two steps: first, lists of crystals 
with significant energy deposits are formed; and 
second, clusters are identified. The EMC typi- 
cally sends data for ~1400 crystals (of 6580 total). 
The majority of these are caused by electronics 
noise and beam-induced background. For each 
crystal, these data include the peak energy and 
time of the crystal waveform. To filter out noise, 
LSEmc rejects individual crystal signals below an 
energy threshold of 20 MeV or which lie outside 
a 1.3 /xs time window around the event time. For 
the remaining crystals, raw energies and times are 
converted into physical units and added to the 
L3Emc crystal list. Clusters are formed using an 
optimized look-up-table technique requiring only 
a single pass over the crystal list. Clusters with a 
total energy above 100 MeV are retained, and the 
energy weighted centroid and average time, the 
number of crystals, and a lateral moment describ- 



ing the shower shape for particle identification are 
calculated. 

11.5.3. Level 3 Filters 

Based on the L3 tracks and clusters, a variety 
of filters perform event classification and back- 
ground reduction. The logging decision is pri- 
marily made by two orthogonal filters, one based 
exclusively on DCH data and the other based only 
on EMC data. 

The drift chamber filters select events with one 
tight (high pt) track or two loose tracks originat- 
ing from the IP, respectively. To account for the 
fact that the IP is not exactly at the origin, track 
selection is based on its x~y closest approach dis- 
tance to the IP, , and Zq^, the corresponding 
z coordinate for that point. The IP position is 
a fixed location close to the average beam po- 
sition over many months. The high pt track is 
required to have a transverse momentum of pt > 
600 MeV/c and to satisfy a vertex condition de- 
fined as Idp^l < 1.0cm, and \zq^ — zip\ < 7.0cm. 
Two tracks are accepted with pt>250MeV/c and 
a somewhat looser vertex condition defined as 
14^1 < 1.5cm, \z^P - zip\ < 10.0cm. 

Two calorimeter cluster filters select events 
with cither high energy deposits or high cluster 
multiplicity. Each filter also requires a high effec- 
tive mass calculated from the cluster energy sums 
and the energy weighted centroid positions of all 
clusters in the event assuming massless particles. 
The first filter requires at least two clusters of 
EcM > 350 MeV (cm. system energy) and event 
mass greater than 1.5 GeV; the second filter re- 
quires at least four clusters, and an event mass 
greater than 1.5 GeV. 

At current luminosities, the output of both the 
DCH and EMC filters is dominated by Bhabha 
events, which need to be rejected. This is ac- 
complished by a Bhabha veto filter that selects 
one-prong (with only a positron in the backward 
part of the detector) and two-prong events (with 
both e+ and e~ detected). Stringent criteria on 
EMC energy deposits are imposed, relying on the 
track momenta and on E/p. The two-prong veto 
requires either colinearity between the tracks in 
the cm. system or an acolinearity that is consis- 
tent with initial state radiation (ISR). 
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For purposes of calibration and offline luminos- 
ity measurements, Bhabha, radiative Bhabha, 77 
final state, and cosmic ray events are flagged. The 
output rate of flagged Bhabha events is adjusted 
to generate an approximately flat distribution of 
events in polar angle. Radiative Bhabha events 
are identified by selecting two-prong events with 
missing energy and requiring an EMC cluster in 
the direction of the missing momentum. Events 
with two high energy clusters, back-to-back in the 
cm. system select the e+e^ 77 process. The 
cosmic ray selection is DCH-based and requires 
two back-to-back tracks in the laboratory frame 
with nearly equal impact parameters and curva- 
ture. A significant background from ISR Bhabha 
events faking this topology is removed using the 
same kinematic constraints used in the two-prong 
veto. 

The online luminosity monitoring and energy 
scale monitoring are performed in L3. A track- 
based lepton-pair selection with a well known ef- 
ficiency monitors the luminosity. Hadronic filters 
for selection of continuum and _B_B-enriched sam- 
ples monitor the energy scale. The latter two cat- 
egories are distinguished by an event shape selec- 
tion using a ratio of Fox- Wolfram moments ]114| . 
The ratio of the SB-enriched sample to the lu- 
minosity is a sensitive measure of relative posi- 
tion on the T{4:S) peak and thereby monitors the 
beam energies. 

11.6. Level 3 Performance 

and Operational Experience 

The L3 trigger efficiency for Monte Carlo sim- 
ulated events are tabulated in Table ^ for events 
passing Level 1. High efficiencies are indepen- 
dently achieved for the DCH and EMC based fil- 
ters applied to simulated hadronic events. The 
comparison between data and Monte Carlo L3 
trigger pass fractions for the various filters also 
show good agreement when requiring tracking, 
and EMC based hadronic event selections in turn. 

An example of the event display used for online 
trigger monitoring is shown in Figure |8^. L3 re- 
constructed tracks and EMC clusters are shown 
together with the LI and L3 trigger line states 
for the event. The left column lists the LI trig- 
ger lines and their states: on (1); off (0); or on 



Table 19 

Composition of the L3 output at a luminosity of 
2.6x10^3 cm-2s-i. 



Event type Rate (Hz) 

Hadrons, rr, and /i/i 16 

Other QED, 2-photon events 13 

Unidentified Bhabha backgrounds 18 

Beam-induced backgrounds 26 

Total physics accept 73 



Calibration Bhabhas (e+e^) 30 

77, Radiative Bhabhas (e+e^7) 10 

Random triggers and cosmic rays 2 

L1,L3 pass through diagnostics 7 

Total calibration/diagnostics 49 



but ignored due to prescale factor (-1). The right 
column shows the same information for the L3 
trigger lines. 

For a typical run on the T(4S') peak with an 
average luminosity of 2.6x 10^'^ cm~^s~^, the L3 
event composition is tabulated in Table The 
desired physics events contribute 13% of the total 
output while the calibration and diagnostic sam- 
ples comprise 40%. 

The L3 executable currently takes an average 
processing time of 8.5 ms per event per farm com- 
puter. A Level 1 input rate of 2700 Hz saturates 
the Level 3 processors, well above the 2 kHz de- 
sign requirement. At this input rate the L3 pro- 
cess consumes ~72% of the CPU time, the rest is 
spent in OEP, including the network event builder 
and in the operating system kernel. 

11.7. Summary and Outlook 

Both the LI and L3 trigger systems have met 
their original design goals at a luminosity of 
3 X lO'^'^ cm~^s^^. The triggering efficiencies for 
BB events generally meet the 99% design goal for 
both LI and L3. The orthogonal triggers based on 
DCH-only and EMC-only information have suc- 
cessfully delivered stable and measurable overall 
trigger efficiency. The current system also pro- 
vides a solid foundation for an upgrade path to 
luminosities of 10'^'* cm^^s^^ or more. 

Short-term LI trigger improvements will pri- 
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Figure 88. A Level 3 event display. The small circles and small crosses in the DCH volume are 
DCH hits and TSF segment hit wires respectively. The filled EMC crystals represent energy deposit 
(full crystal depth — 2 GeV) from Level 3 EMC clusters while the small triangles just inside the EMC 
indicate the location of the cluster centroid. 



marily come from further background rejection, 
afforded by algorithm refinements and upgrades 
of the DCT. This is essential for reducing the load 
on the DAQ and L3. The new PTD algorithm wiU 
effectively narrow the track do acceptance win- 
dow, while a new BLT algorithm will narrow the 
track zq acceptance. 

For the longer term future, a major DCT up- 
grade is planned. By adding the stereo layer in- 
formation, a Zq resolution of 4 cm is expected, al- 
lowing for an efficient rejection of beam-induced 
background beyond z = ±20 cm. 



Future improvements for L3 will also empha- 
size background rejection. Improvements in the 
L3 IP track filter are expected to further reduce 
beam-induced background to about one third of 
current levels. The physics filter algorithms will 
be tuned and improved, primarily for rejecting 
Bhabha, QED, and two-photon events. Improve- 
ments in the L3 tracking algorithms are expected 
to lower the pt thresholds below 250MeV/c. A 
moderate CPU upgrade for the L3 online farm 
will be sufficient to keep up with luminosities of 
~ lO^-^cm-^s-i. 



94 



Table 18 

L3 trigger efficiency (%) for various physics processes, derived from Monte Carlo simulation. 



L3 Trigger 


^BB 






^CC 


^uds 




1 track filter 


oy.y 


oy.y 


oO.O 




OO.Z 




2 track filter 


98.9 


84.1 


94.5 
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12. The Online Computing System 

12.1. Overview 

The BABAR online computing system comprises 
the data acquisition chain from the common FEE, 
through the embedded processors in the data ac- 
quisition system and the L3 trigger, to the logging 
of event data. It also includes those components 
required for detector and data acquisition control 
and monitoring, immediate data quality monitor- 
ing, and online calibration. 

12.1.1. Design Requirements 

The data acquisition chain was designed to 
meet the following basic performance require- 
ments. It must support a LI trigger accept rate 
of up to 2 kHz, with an average event size of 
~32 kbytes and a maximum output (L3 trigger 
accept) rate of 120 Hz. While performing these 
functions it should not contribute more than 
a time-averaged 3% to deadtime during normal 
data acquisition. 

The online system is also required to be capa- 
ble of performing data acquisition simultaneously 
on independent partitions — sets of detector sys- 
tem components — to support calibrations and di- 
agnostics. 

Normal detector operation, data acquisition 
and routine calibrations are performed efficiently 
and under the control of a simple user interface 
with facilities for detecting, diagnosing, and re- 
covering from common error conditions. 

Following standard practice, the event data ac- 



quired by the system are subjected to monitoring. 
Such monitoring is configurable by experts and 
designed to detect anomalies in the detector sys- 
tems which, if present, are reported to operators 
for rapid assessment and, if necessary, corrective 
action. 

Environmental conditions of the detector, such 
as the state of low and high voltage power, the pu- 
rity of gas supplies, and the operating conditions 
of the accelerator, such as beam himinosity and 
currents, are measured and recorded in a fash- 
ion that permits the association with the event 
data logged. Conditions relevant to data quality 
are monitored for consistency with specified stan- 
dards. Operators are alerted if these are not met. 
Data-taking is inhibited or otherwise flagged if 
conditions are incompatible with maintaining the 
quality of the data. 

Operational configurations, calibration results, 
active software version numbers, and routine mes- 
sages and error messages are also recorded. Dur- 
ing data analysis or problem diagnosis, these data 
help in reconstructing the detailed operating con- 
ditions. 

12.1.2. System Components 

The online computing system is designed as a 
set of subsystems using elements of a common 
software infrastructure running on a dedicated 
collection of hardware. 

The major subsystems arc: 

• Online Dataflow (ODF) — responsible for 
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communication with and control of the de- 
tector systems' front-end electronics, and 
the acquisition and building of event data 
from them; 

• Online Event Processing (OEP) — 
responsible for processing of complete 
events, including L3 (software) triggering, 
data quality monitoring, and the final 
stages of calibrations; 

• Logging Manager (LM) — responsible for re- 
ceiving selected events sent from OEP and 
writing them to disk files for use as input to 
the Online Prompt Reconstruction process- 
ing; 

• Online Detector Control (ODC)— 
responsible for the control and monitoring 
of environmental conditions of the detector 
systems; 

• Online Run Control (ORC) — ties together 
all the other components, and is responsi- 
ble for sequencing their operations, inter- 
locking them as appropriate, and providing 
a graphical user interface (GUI) for opera- 
tor control. 

Each of these components, as well as a selection 
of the common tools which tie them together are 
described below. 

The entire system is coded primarily in the 
C-|— I- language, with some use of Java for graph- 
ical user interfaces. Object-oriented analysis and 
design techniques have been used throughout. 
This has been an important feature, enhancing 
development speed, maintainability, and extensi- 
bility. 

12.1.3. Hardware Infrastructure 

The hardware infrastructure for the online sys- 
tem is shown schematically in Figure |8^. 

The data from the FEEs of the various detector 
systems are routed via optical fiber links to a set 
of 157 custom VME Readout Modules (ROMs). 
These ROMs are grouped by detector system and 
housed in 23 data acquisition VME crates that 
are controlled by the ODF software. One ROM 
in each crate aggregates the data and forwards 



them for event building to 32 commercial Unix 
workstations |115| which are part of the online 
farm. Other farm machines perform data moni- 
toring and calibrations. The crates and farm ma- 
chines communicate via full-duplex lOOMbits/s 
Ethernet, linke d by a network switch — the event 
builder switch [ 116 |. The ROMs are supported by 
a boot server providing core and sy stem -specific 
code and configuration information | 117 ]. 

The thirty-two online farm machines host the 
OEP and L3 trigger software. The events ac- 
cepted by the trigger are logged via TCP /IP to a 
logging server |117] and written to a disk buffer 
for later reconstruction and archival storage. Var- 
ious data quality monitoring processes run on 
farm machines not used for data acquisition. 

Several additional file servers hold the online 
databases and production software releases. A 
further set of application servers host the central 
functions of the various online subsystems. Op- 
erator displays are s upported by a group of ten 
console servers |118|. 

An additional set of 15 VME crates, each with 
an embedded processor, contain the data acquisi- 
tion hardware for the detector control subsystem. 

All VME crates, the online farm, and all the 
application and console servers are connected via 
a switched 100 Mbits/s Ethernet network distinct 
from that used for event building, with 1 Gbits/s 
fiber Ethernet used for the file servers and inter- 
switch links. 

12.1.4. User interaction 

Operator control of the online system is 
achieved primarily through a custom Motif GUI 
for run control and an extensive hierarchy of dis- 
plays for detector control, including control pan- 
els, strip charts and an alarm handler. An elec- 
tronic logbook is made available through a Web 
browser interface. These and other GUIs are or- 
ganized across seventeen displays for the use of 
the experiment's operators. This operator envi- 
ronment provides for basic control of data acquisi- 
tion, the overall state of the detector, and certain 
calibration tasks. 

Each detector system has developed a set of 
specialized calibration and diagnostic applica- 
tions using the tools provided in the online sys- 
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Figure 89. Physical infrastructure of the BABAR onhne system, including VME crates, computers, and 
networking equipment. 



tern. A subset of these calibrations has been spec- 
ified to be run once per day, during a ten-minute 
scheduled beam-off period. The run control logic, 
combined with the capability for creating parti- 
tions, allows calibrations for all detector systems 
to be run in parallel and provides the operator 
with basic feedback on the success or failure of 
each. 

12.2. Online Dataflow 

ODF handles data acquisition and processing 
from the detector systems' FEE through the de- 
livery of complete events to the online farm ]119| . 
The ODF subsystem receives the LI trigger out- 
puts, filters and distributes them to the FEE, 
reads back the resulting data and assembles them 
into events. It provides interfaces for control of 
data acquisition, processing and calibration of de- 
tector system data, and FEE configuration. Mul- 
tiple independent partitions of the detector may 
be operated simultaneously. 

Event data acquisition proceeds from a trig- 
ger decision formed i n the Fast Control and Tim- 
ing System (FCTS) |120| based on inputs from 



the LI trigger. LI Accepts are distributed, in 
the full detector configuration, to the 133 ROMs 
connected via optical fibers to the detector sys- 
tem FEE. These ROMs read and process the data 
from the FEE. One to ten such ROMs from a sin- 
gle detector system are located in each of the data 
acquisition VME crates. ODF builds complete 
events from these ROMs, first collecting the data 
in each crate into an additional dedicated ROM, 
and then collecting the data from the 23 of these, 
across the event builder network switch, into the 
online farm. 

The operation of the system is controlled by 
ODF software running on one of the application 
servers, under the direction of run control. A sin- 
gle ROM in the VME crate containing the central 
FCTS hardware supports the software interface 
to ODF. The distribution of ROMs by detector 
system is shown in Table ^ The numbers of 
ROMs is shown as a sum of those connected di- 
rectly to the detector FEE, and those used for 
event building. 

AU of the ROM CPUs boot via NFS over the 
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Table 20 

VME crates and ROMs used by ODF 
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event building network from the boot server de- 
scribed above. About 1.5 Mbytes of core ODF 
code plus another ~4 Mbytes of detector-specific 
code are loaded into each ROM. This, along with 
the booting process, takes about 40 seconds. 

The ODF software allows all the components 
of this heterogeneous system to be represented in 
a uniform object-oriented application framework. 
These components are organized into five levels 
which map closely onto the physical structure. 

For each component at each level, its operation 
is abstracted as a finite state machine. The com- 
plete set of these machines is kept coherent by 
passing messages and data regarding state tran- 
sitions along the chain of levels. The basic flow 
of control and data is shown in Figure 1^. The 
mapping of levels to components is as follows: 

Control The Unix-based process controlling 
the operation of each partition and the source 
of all state transitions except for LI Accept. It 
transmits state transition messages over the net- 
work to the source level, waiting for acknowledge- 
ment of their successful processing by all levels. 

Source The FCTS hardware and the software 
running in the ROM located in the FCTS VME 
crate. For each partition, its source level receives 
control level transitions and LI trigger outputs 
and distributes them via the FCTS hardware to 
all ROMs in the VME crates included in the par- 
tition. LI triggers are modeled in the subsystem 



as an additional, idempotcnt state transition, LI 
Accept, and are treated uniformly with the others 
wherever possible. 

Segment The ROMs connected to the detec- 
tor FEE, with their ODF and detector system- 
specific software. Each segment level ROM re- 
ceives state transition messages from the source 
level and runs appropriate core and detector 
system-specific tasks in response. These tasks in- 
clude the acquisition of raw data from the FEE 
in response to LI Accepts, and feature extraction. 
Output data resulting from this processing is at- 
tached to the transition messages, which are then 
forwarded over the VME backplane to the frag- 
ment level ROM in each crate. 

Fragment The per-crate event builder ROMs 
and software. The single fragment level ROM 
in each crate aggregates the messages from the 
crate's segment level ROMs — the first stage of 
event building — and forwards the combined mes- 
sage to one of the event level Unix nodes. 

Event The processes on the online farm nodes 
receiving complete events and handing them over 
to OEP for filtering and logging. The ODF event 
level code aggregates messages, with their at- 
tached data, from all the crates in a partition — 
the second and final stage of event building. The 
resulting data may be further processed by user 
code in the event level, but are normally just 
passed on to OEP. The control level is noti- 
fied of the completion of processing of all tran- 
sitions other than LI Accept. Both the fragment 
and event level event builders use a data-driven 
"push" model, with a back pressure mechanism 
to signal when they are unable to accept more 
data. 

Test stands of varying complexity are sup- 
ported. The simplest possible consists of a single 
Unix machine which runs both control and event 
level code, with two FCTS modules and a sin- 
gle ROM, running source, segment, and fragment 
level code, in one VME crate. Configuration is 
detected at run-time, so the same code that runs 
in the full system can also run in test stand sys- 
tems. 
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Figure 90. Schematic of the ODF levels, their mapping onto physical components, and the flow of control 
signals and data between them. 



12.2.1. Control and Source Levels 

The control level sends state transition mes- 
sages for a partition over the net work , using 
the User Datagram Protocol, UDP [121|, to the 
source level in the single ROM inside the FCTS 
crate. In the source level, the transition message 
is sent over VME to an FCTS module which for- 
wards it as a 104-bit 59.5 MHz serial word to all 
VME crates in the relevant partition. This serial 
word contains a 56-bit event time stamp (count- 
ing at 59.5 MHz), a 32-bit transition-specific word 
and additional control bits. LI Accept transitions 
and calibration sequences, however, originate in 
the source level and the same mechanism is used 
to transmit them through the system. 

The FCTS hardware receives the 24 LI trigger 
output lines and eight additional external trigger 
lines. The FCTS crate is a 9U VME crate, with a 
custom P3 backplane on which all the trigger lines 
are bussed. For each partition, an FCTS module 
receives these lines. It is configurable with a bit 
mask specifying the trigger lines enabled for its 
partition, and an optional prescale factor for each 
line. A trigger decision is formed for the partition 
by taking the logical OR of the enabled prescaled 
lines. Twelve of these modules are installed in the 
full system, thus setting its maximum number of 
partitions. A detector system can belong to only 
one partition at a time. 

The FCTS crate receives two timing signals 
from the accelerator: a 476 MHz clock tied to 
the RF structure of PEP-II and a 136 kHz fidu- 
cial that counts at the beam revolution fre- 
quency. The former is divided by eight to create 
a 59.5 MHz system clock. The fiducial is used to 
start timing counters and to check the synchro- 
nization of the clocks. 

There are two types of deadtime in the ODF 



subsystem. The first arises from the minimum 
2.7 /xs spacing between LI Accept transitions. 
This restriction simplifies the logic design of the 
FEE readout, because each signal in the silicon 
tracker and drift chamber is thus associated with 
only one LI Accept. The FCTS hardware enforces 
this minimum separation between transitions, in- 
troducing an irreducible, yet minimal dead time 
of 0.54% at 2 kHz. 

The second type of deadtime arises when all 
FEE buffers are full and thus unable to accept 
another event. In a time required to be less than 
the inter-command spacing, each VME crate in 
a partition may send back a full signal indicat- 
ing that it is no longer able to process further 
LI Accept transitions. The FCTS hardware de- 
tects these signals and disables triggering until 
the FEE are once again prepared to accept data. 

An actual LI Accept signal is only generated 
from a partition's trigger decision when neither 
form of dead time is asserted. 

12.2.2. Segment and Fragment Levels 

The segment and fragment levels reside in the 
23 detector system VME crates. These are stan- 
dard 9U crates with a custom P3 backplane. 

The 104-bit serial transition messages that 
leave the source level are received by a FCTS 
module in each VME crate in a partition. This 
module in turn forwards these messages to the 
ROMs in the crate over the custom backplane, 
along with the 59.5 MHz system clock. 

A ROM consists of four components (see Fig- 
ure pl|), a commercial single-board computer 
(SBC) 122 and three custom boards. The cus- 
tom boards include: a controller card for receiv- 
ing FCTS commands and supporting FEE reads 
and writes; a personality card that transmits com- 
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Table 21 

Typical event sizes from detector systems 



FULL 
S9.S MHz Clock 
FCTS Transition Messages 



Figure 91. A ROM with a triggered personality 
card (TPC) 

mands to and receives data from the FEE; and a 
PCI mezzanine card with a 33 MHz Intel i960 I/O 
processor. The SBCs run the VxWorks [| oper- 
ating system with custom code written in CH — h 
and assembly language. 

There are two styles of personality cards in the 
system: triggered (TPC) and untriggered (UPC). 
UPCs are used only in the EMC system. UPCs 
accept data continuously from the FEE into a 
buffer pipeline, at a rate of 3.7 MHz. From these 
samples EMC trigger information is derived and 
sent over a dedicated serial link to the trigger 
hardware, providing it with a continuous data 
stream. An LI Accept causes up to 256 samples 
of the raw data stream to be saved to an inter- 
mediate memory on the UPC. 

A TPC (used in all other systems) reads out 
FEE data only when an LI Accept signal is re- 
ceived, again saving it into an intermediate mem- 
ory. Each detector reads out data in a time win- 
dow around the trigger signal, large enough to 
allow for trigger jitter and detector time resolu- 
tion. For instance, this window is about 500 ns 
wide for the SVT. The actual event time within 
this window is estimated in the L3 trigger soft- 
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ware and then refined offline in the course of full 
event reconstruction. 

FEE commands are sent and data received 
by the personality cards over uni-directional 

All 



1.2Gbits/s serial optical fiber links |123]. 
FEEs provide zero suppression in hardware ex- 
cept in the EMC and IFR. Data are transferred 
from the personality card to the SBC memory us- 
ing the i960 as a direct memory access (DMA) 
engine. This DMA runs at nearly the ideal 
133Mbytes/s rate of the PCI bus. 

The FEEs for various systems are able to buffer 
data for three to five LI Accept transitions. The 
ROM keeps track of the buffer occupancy and 
sends, when necessary, a full signal as previously 
described. The full condition is removed when 
event reading by the ROM frees sufficient buffer 
space. This mechanism handles back pressure 
from any stage of the data acquisition through 
to data logging by OEP. 

The ODE application framework provides uni- 
form software entry points for the insertion of user 
code at each level of the system. This capability 
is used primarily at the segment level, for FEE 
configuration and feature extraction. Table |2^ 
presents typical data contributions from each de- 
tector system and the trigger. 

Data from the segment level ROMs in a crate 
are gathered by the fragment level ROM using 
a chained sequence of DMA operations. The 
maximum throughput of the fragment level event 
builder is about 31 Mbits/s. 
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In calibrations, ODF may be operated in a 
mode in which LI Accept data are not transferred 
out of the segment level ROMs. This allows for 
calibration data accumulation at high rates inside 
the ROMs, not limited by the throughput of the 
event builders or any downstream software. Com- 
pleted calibration results are computed, read out, 
and written to a database. 

12.2.3. Event Level 

For each LI Accept transition passing through 
the ODF subsystem, all fragment ROM data are 
sent to one of the farm machines. The destination 
is chosen by a deterministic calculation based on 
the LI AccepVs 56-bit time stamp, available from 
the FCTS in each ROM. This technique produces 
a uniform quasi-random distribution and intro- 
duces no detectable inefficiency. Events sent to a 
farm machine still busy with a previous event are 
held in a buffer to await processing. 

All fragment data for an event are sent over 
the switched lOOMbits/s Ethernet event building 
network to the selected farm machine. The con- 
nectionless UDP was chosen as the data trans- 



port protocol 1 124 1, allowing a flow control mech- 
anism to be tailored specifically to this applica- 
tion. Dropped packets are minimized by the net- 
work's purely point-to-point, full duplex switched 
architecture, and by careful tuning of the buffer- 
ing in the network switch and other parameters. 
The rare instance of packet loss is detected by the 
event builder and the resulting incomplete event 
is flagged. 

The event level provides the standard software 
entry points for user code. During normal oper- 
ation, these are used only to transfer events via 
shared memory to the OEP subsystem for L3 trig- 
gering, monitoring, and logging. 

12.2.4. System Monitoring 

It is critical that the clocks of the FEEs stay 
synchronized with the rest of the system. Each 
FEE module maintains a time counter which is 
compared to the time stamp of each LI Accept in 
order to ensure that the system remains synchro- 
nized. If it becomes unsynchronizcd, a special 
synch command can be sent through the FCTS, 
causing all systems to reset their clocks. 



To ensure that the data from the correct event 
is retrieved from the FEEs, a five-bit counter is in- 
cremented and sent from the FCTS to the FEEs 
with each LI Accept. These bits are stored in 
the FEEs along with the data and are compared 
on read-back. If they disagree, a special clear- 
readout command is sent which resynchronizes 
ROM buffer pointers with FEE buffer pointers. 

All transitions, including LI Accept^ are logged 
in a 4 kbytes-deep by 20 byte- wide FIFO as they 
pass through the FCTS crate. The transition 
type, the event time stamp, a bit list of the trigger 
lines contributing to the decision, and the current 
FULL bit list from all VME crates are recorded in 
this FIFO. There are also scalers which record 
delivered and accepted luminosity, deadtime due 
to the 2.7 /is minimum inter-command spacing, 
deadtime caused by VME crates being full and 
triggers on each line. These FIFOs and scalers are 
read out by the FCTS ROM, which then trans- 
mits the data to monitoring programs that calcu- 
late quantities such as luminosity, deadtime and 
trigger rates. The UDP multicast protocol ]125| | is 
used to allow efficient simultaneous transmission 
of data to multiple clients. 

To provide diagnostics, a system which multi- 
casts additional performance information on de- 
mand from each CPU, typically at 1 Hz, is used. 
This information is currently received by a single 
client on one of the Unix application servers and 
archived. It can be retrieved subsequently to in- 
vestigate any unusual behaviour observed in the 
system. 

12.3. Online Event Processing 

The online event processing (OEP) subsys- 
tem provides a framework for the processing of 
complete events delivered from the ODF event 
builder |126|. The L3 software trigger operates 
in this framework, along with event-based data 
quality monitoring and the final stages of online 
calibrations. Figure ^ shows the basic flow of 
data in the OEP subsystem. 

The OEP subsystem serves as an adapter be- 
tween the ODF event builder interface and the ap- 
plication framework originally developed for the 
offline computing system. Raw data delivered 
from the ODF subsystem arc put into an object- 
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oriented form and made available through the 
standard event data analysis interface. 

The use of this technique permits the L3 trig- 
ger and most of the data quality monitoring soft- 
ware to be written and debugged within the of- 
fline environment. This software is decomposed 
into small, reusable units — modules, pluggable 
software components in the framework — many of 
which are shared among multiple applications. 

The OEP interfaces allow user applications to 
append new data blocks to the original raw data 
from ODF. The results of L3 event analysis are 
stored in this manner so that the trigger decision 
and the tracks and calorimeter clusters on which 
it is based may be used in later processing, such as 
reconstruction and trigger performance studies. 

Histograms and other monitoring data are ac- 
cumulated across the farm. A distributed his- 



tograming package (DHP) |127| was developed to 
provide networked clients with a single view of 
histograms and time history data. This data is 
summed across all nodes via CORBA-based com- 



munication protocols |12S,129|. 

The fast monitoring system provides auto- 
mated comparisons of monitoring data against 
defined references. Statistical comparisons of live 
histograms, or the results of fits to reference his- 
tograms, analytic spectra, or nominal values of 
fit parameters may be performed at configurable 
time intervals. Comparison failures, tagged with 
configurable severity levels based on the confi- 
dence levels of the comparisons, are displayed to 
operators and logged in the common occurrence 
database, described below. 



The Java Analysis Studio (JAS) package [ 130 | 
previously developed at SLAC was enhanced with 
the ability to serve as a DHP client. It is used for 
viewing of monitoring data. This feature was im- 
plemented by devising a Java server that adapts 
the DHP protocol to the native JAS data proto- 
col. 

In addition to the primary triggering and mon- 
itoring functions carried out on 32 online farm 
machines, OEP provides a "trickle stream" pro- 
tocol that allows networked clients to subscribe to 
a sampling of the event data. This scheme pro- 
vides support for event displays and additional 
detailed data quality monitoring ( "Fast Monitor- 



ing" ) . A dedicated operator console supports the 
JAS-based data quality monitoring system. This 
console is used to display histograms from Fast 
Monitoring and the L3 trigger processes, along 
with any error conditions detected by the auto- 
matic histogram analysis facility. 

12.4. Data logging 

Events selected by the L3 trigger algorithms 
in OEP are retained for subsequent full recon- 
struction. The events are sent from the 32 OEP 
nodes via TCP /IP to a single multithreaded pro- 
cess, the Logging Manager (LM), running on the 
logging server. The LM writes these data to 
RAID storage arrays in a format specific to OEP. 
Data from all 32 nodes are combined into a sin- 
gle file for each data-taking run (typically two to 
three hours of data acquisition, resulting in files 
of about 15-20 Gbytes in size). 

Completed data files are copied to the S LAC 
High Performance Storage System (HPSS) 131 
system for archiving to tape. Within eight hours 
of data acquisition these files are retrieved from 
HPSS for event reconstruction. The data files are 
also retrievable for other tasks such as detector 
system hardware diagnostics and offline tests of 
the L3 trigger algorithms. 

12.5. Detector Control 
12.5.1. Design Principles 

The Experimental Physics and Industrial Con- 
trol System, EPICS §, was selected to provide 
the basis for the ODC subsystem. This provides 
direct connection to the electrical signals of the 
power supplies and other hardware, with suffi- 
cient monitoring and control to allow commission- 
ing, fault diagnosis, and testing. A summary of 
monitoring and control points is presented in Ta- 
ble n 

Beyond the writing of custom drivers, only mi- 
nor additions or changes were required to EPICS. 
EPICS and the additional BABAR-specific soft- 
ware are written in the C language. 

Detector-wide standard hardware was adopted 
to ease development and maintenance. The stan- 
dard ODC crate is a 6U VME chassis contain- 



ing a single board computer |132| serving as an 
EPICS input/output controller (IOC). Fifteen 
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Figure 92. Flow of data in the OEP subsystem: ODF event level (EL) and L3 trigger processes on 
each OEP node; the Logging Manager (LM) on the logging server; the DHP "requestor" process that 
combines histograms from all 32 L3 processes; one instance of a Fast Monitoring (FM) process with DHP 
histograms; the Java server that makes DHP histograms available to JAS clients; two such clients, and 
one event display for the L3 trigger. OEP-specific data transport protocols are identified. 



Table 22 

ODC Distribution by system of approximately 12,000 recorded monitor channels 

SVT DCH DIRC EMC IFR Central 



Radiation Dose 


12 


8 


22 


116 






Data Rate 






12 




1612 


20 


Temperature 


208 


87 


36 


506 


146 


100 


Humidity 4 


5 


12 


7 






2 


Magnetic Field — 




8 








50 


Position 


30 


22 










Gas System 


4 


115 


12 


1 


32 




Fluid System 


20 


2 


18 


18 


20 


90 


Liquid Source 








12 




3 


HV System 


2080 


1299 


672 


24 


1574 




LV System (non-VME) 




62 


1080 


442 


94 




VME Crates 


5 


2 


12 


12 


8 


1500 


CAN Micro Controller 


48 


40 


3 


47 


45 




Finisar Monitor 


28 


12 


12 









System Totals 



2439 1654 1899 1185 3531 1765 
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such crates are used in the experiment. EPICS 
is fully distributed. For example, each IOC sup- 
plies its own naming service, notify-by-exception 
semantics, and processing. The lOCs boot from 
a dedicated server. 

Analog data are either digitized by modules 
within the crates or, more commonly, on digitizer 
boards located directly on the detector. In the 
latter case, the CANbus standard |132] is used 
for the transport of signals to and from the de- 
tector. A c ustom "general monitoring board" 
(GMB) [13^ was developed to interface CANbus 
to the on-detector electronics. The GMB contains 
a microcontroller, an ADC, multiplexors, and op- 
erational amplifiers. It can digitize up to 32 sig- 
nals. 

12.5.2. User Interface 

The operator view of this part of the control 
system is via screens controlled by the EPICS dis- 
play manager (DM). Dedicated control and dis- 
play panels were developed using DM for each of 
the detector systems, using common color rules 
to show the status of devices. A top-level panel 
for GDC summarizes the status of all systems and 
provides access to specialized panels. 

The EPICS alarm handler with some BABAR- 
specific modifications is used to provide oper- 
ators with audible and color-coded alarms and 
warnings in a hierarchical view of all the systems 
and components. Conditions directly relevant to 
personnel or detector safety are further enforced 
by hardware interlocks, the status of which are 
themselves displayed in a set of uniform EPICS 
screens, in the alarm handler, and on an alarm 
annunciator panel. 

12.5.3. Interfaces to Other 
BABAR Softv^fare 

A custom C++ layer above EPICS consisting 
of Component Proxies and Archivers provides for 
device-oriented state management and archival 
data collection. This is ODG's interface to the 
rest of the online system. 

The 27 component proxies (CPs), running on 
a Unix application server, each define a logical 
component representing some aspect of a detec- 
tor system or the experiment's central support 



systems, aggregated from the < lO'"^ individual 
EPICS records. 

The CPs present a simple finite state ma- 
chine model as their interface to Run Control. 
The most important actions available are Config- 
ure, on which the CP accesses the configuration 
database, retrieving set points for its component's 
channels, and Begin Run, which puts the CP into 
the Running state, in which setpoint changes are 
prohibited and readbacks are required to match 
settings. While in the Running state, the CP 
maintains a Runnahle flag which reflects that re- 
quirement and allows Run Control to ensure that 
data acquisition is performed only under satisfac- 
tory conditions. 

The CP's other principal function is to provide 
an interface for the rest of the online system to 
the ambient data collected by ODC on the state 
of the detector hardware and its environment. It 
is the task of the archiver processes, each paired 
with a CP, to collect the ambient data, aggregate 
them and write out histories approximately ev- 
ery hour to the ambient database. These recorded 
data are associated with times so that they may 
be correlated with the time stamps of the event 
data. Data from the archiver processes or from 
the database may be viewed with a custom graph- 
ical browser. 

Ambient data typically vary only within a nar- 
row noise range or dead-band. The storage of 
unnecessary data is avoided by recording only 
those monitored quantities which move outside of 
a per-channel dead-band range or across an alarm 
threshold. 

12.5.4. Integration With the Accelerator 

Close integration between the BABAR detector 
and the PEP-II accelerator is essential for safe 
and efficient data collection. Data from the accel- 
erator control system are transferred via EPICS 
channel access to BABAR for display and storage, 
managed by a dedicated CP. In turn, background 
signals from the detector are made available to 
PEP-II to aid in injection and tuning, minimiz- 
ing backgrounds, and optimizing integrated lu- 
minosity. An important component of this com- 
munication is the "injection request" handshake. 
When the PEP-II operator requests a significant 
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change in the beam conditions, such as injection, 
the request can only procede following confirma- 
tion from BABAR. This procedure complements 
the safety interlocks based on radiation dose mon- 
itors. 

12.5.5. Operational Experience 

The ODC subsystem has been operational since 
the initial cosmic-ray commissioning of the detec- 
tor and the beginning of data-taking with col- 
liding beams. The core EPICS infrastructure 
has proven to be very robust. The large size 
of the subsystem, with its 15 lOCs and < 10^ 
records, produces heavy but manageable traffic 
on the experiment's network. The rate of data 
into the ambient database averages 4.6 Mbytes/hr 
or 110 Mbytes/day 

12.6. Run Control 

The ORG subsystem is implemented as an ap- 
plication of SMI-I--I-, a toolkit for designing dis- 
tributed control systems [135|. Using this soft- 



ware, the BABAR experiment is modeled as a col- 
lection of objects behaving as finite state ma- 
chines. These objects represent both real entities, 
such as the ODF subsystem or the drift chamber 
high voltage controller, and abstract subsystems 
such as the "calibrator," a supervisor for the co- 
ordination of online components during detector 
calibration. The behavior of the objects are de- 
scribed in a specialized language (SML) which is 
interpreted by a generic logic engine to implement 
the control system. 

The SML descriptions of the objects which 
make up the experiment simply specify their own 
states and transitions as well as the connections 
between the states of different objects. Ob- 
jects perform actions on state transitions, which 
may include explicitly commanding transitions in 
other objects; objects may also be programmed 
to monitor and automatically respond to changes 
of state in other objects. Anticipated error con- 
ditions in components of the online system are 
reflected in their state models, allowing many er- 
rors to be handled automatically by the system. 
To reduce complexity, logically related objects are 
grouped together into a hierarchy of cooperating 
domains. 



The system is highly automated; user input is 
generally required only to initialize the system, 
start and stop runs, and handle unusual error 
conditions. The user communicates with ORG 
via a configurable Motif-based GUI included in 
SM1++. 

The states and behavior of ORG objects repre- 
senting external systems are provided by a spe- 
cial class of intermediate software processes called 
proxies. A proxy monitors its system, provides 
an abstraction of it to ORG, and receives state 
transition commands. These commands are in- 
terpreted and applied to the underlying hardware 
or software components, implementing the tran- 
sitions' actions. The control level of an ODF par- 
tition is an example of such a proxy. 

Gommunication between the various proxies 
and the ORG engines is provided by DIM |136| ], 
a fault tolerant "publish and subscribe" commu- 
nications package based on TGP/IP sockets, al- 
lowing ORG to be distributed transparently over 
a network. 

Essential to the operation of the online sys- 
tem is the notion of the Runnahle status of its 
various ODG and data acquisition components, 
indicating that they are in a state suitable for 
production-quality data-taking. The ORG logic 
interlocks data-taking to the logical AND of all 
components' Runnable status. Whenever this 
condition is not satisfied, data-taking may not 
start and any existing run will be paused with 
an alert sent to the operator. 

12.7. Common Software Infrastructure 
12.7.1. Databases 

Five major databases are used by the online 
system: 

1. Configuration Database: This database, im- 
plemented using the commercial object-oriented 
database management system Objectivity Q, al- 
lows the creation of hierarchical associations of 
system-specific configuration data with a single 
numeric configuration key. This key is distributed 
to all online components, which can then use it 
to retrieve from the database all the configuration 
information they require. Convenient mnemonics 
are associated with the keys for currently rele- 
vant configurations, and may be selected for use 
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via the ORG GUI |l37 |. 

2. Conditions Database: The Gonditions 
Database is used to record cahbration and ahgn- 
ment constants, and the configuration keys in 
force during data-taking runs. It has the addi- 
tional feature that the data for a given time in- 
terval may be updated as they are refined in the 
cour se of improved understanding of the appara- 
tus ll38| . 

The Gonfiguration and Gonditions Databases 
are both made available for reconstruction and 
physics analysis. 

3. Ambient Database: The Ambient Database 
is used principally by the ODG subsystem to 
record detector parameters and environmental 
data at the time they are measured [137|. 

Both the Ambient and Gonditions databases, 
are implemented using Objectivity, and are based 
on the notion of time histories of various data 
associated with the experiment. The history for 
each item is divided into intervals over which a 
specific value is consistent. 

4. Occurrence (Error) Log: Informational and 
error messages generated in the online system 
are routed through the GMLOG system |139| to 
a central database, from which they are avail- 
able for operators' realtime viewing or historical 
browsing, using a graphical tool, as well as for 
subscription by online components which may re- 
quire notification of certain occurrences. 

5. Electronic Logbook: An Oracle-based [ 140 | 



12.7.2. Softv^fare Release Control 

and Configuration Management 

All of the online software is maintained in the 
common BABAR code repository, based on the 
freely available Goncurrent Versions System soft- 



logbook is used to maintain the history of the 
data-taking, organized by runs. It contains infor- 
mation on beam parameters — instantaneous and 
integrated luminosity, currents, and energies — 
as well as records of data acquisition parame- 
ters such as trigger rates, data volumes, and dead 
times, and the detector configuration used for a 
run. The logbook also contains text comments 
and graphics added by the operations staff. 

A number of other databases are used in the 
online system for various tasks such as indexing 
logged data files, the repair history of online hard- 
ware and spares, and software problem reports. 



ware, GVS p41 |. 

The online's Unix and VxWorks applications 
are built and maintained with an extensi on o f the 
standard BABAR software release tools |142]. At 



the start of every data-taking run, the identities 
of the current production software release and any 
installed patches are recorded; thus it is possible 
at a later date to reconstruct the versions of online 
software used to acquire data. 

12.8. Summary and Outlook 

The online system has exceeded its data acqui- 
sition performance goals. It is capable of acquir- 
ing colliding beam events, with an average size 
of 28 kbytes, at a ~ 2500 Hz LI trigger rate and 
reducing this rate in L3 to the required ~ 120 Hz 
limit. This provides comfortable margins, since 
under normal beam conditions the LI trigger rate 
is 800-1000 Hz. 

The system is capable of logging data at much 
higher rates; the nominal 120 Hz figure represents 
a compromise between data volume and its con- 
sequential load on downstream processing and 
archival storage, and trigger efficiency for low 
multiplicity final states. 

During normal data-taking, the online system 
routinely achieves an efficiency of over 98%, tak- 
ing both data acquisition livetime and the sys- 
tem's overall reliability into account. 

There are several hardware options for enhanc- 
ing ODF capacity. Gurrently most ROMs receive 
more than one fiber from the FEE. These fibers 
could be distributed over more ROMs to add pro- 
cessing power. There are also commercial up- 
grade paths for the ROMs' SBG boards available. 
Grates can be split (up to a maximum total of 32) 
to create more VME event building bandwidth, as 
well as more fragment level GPU power and net- 
work bandwidth. Gigabit Ethernet connections 
could also be installed to improve the network 
event builder's bandwidth. 

Various software upgrade options are being 
investigated, including optimizing the VxWorks 
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network drivers and grouping sets of events to- 
gether in order to reduce the impact of per-event 
overhead. 

Current background projections indicate that 
fragment level CPU, segment level memory bus 
bandwidth, and network event building band- 
width are the most likely bottlenecks for future 
running. 

Increases in the LI trigger rate or in the back- 
ground occupancy and complexity of events are 
expected to necessitate additional capacity for 
OEP, principally for L3 triggering. The online 
farm machines could be replaced with faster mod- 
els. More machines could be added, at the ex- 
pense of increases in coherent loading on various 
servers and of additional management complex- 
ity. 

No significant capacity upgrades to the data 
logging subsystem or to ODC are foreseen at this 
time. 

13. Conclusions 

During the first year of operation, the BABAR 
detector has performed close to expectation with 
a high degree of reliability. In parallel, the PEP-II 
storage rings gradually increased its performance 
and towards the end of the first year of data- 
taking routinely delivered close to design luminos- 
ity. In fact, the best performance surpassed the 
design goals, both in terms of instantaneous as 
well as integrated luminosity per day and month. 
Of the total luminosity of 23.5 fb~^ delivered by 
PEP-II during the first ten months of the year 
2000, BABAR logged more than 92%. The data 
are fully processed with a delay of only a few 
hours. They are of very high quality and have 
been extensively used for physics analysis. 

A large variety of improvements to the 
event reconstruction and detailed simulation are 
presently being pursued. They include improve- 
ments in many aspect of the calibration and re- 
construction procedures and software, for exam- 
ple, the calibration and noise suppression in the 
EMC, and the development of techniques for pre- 
cision alignment of the SVT and DCH. The latter 
effort will not only benefit the overall efficiency 
and precision of the track reconstruction, it will 



also improve the matching of tracks with signals 
in the DIRC and EMC. Detailed studies and the 
full integration of all available information per- 
tinent to the identification of charged and neu- 
tral particles are expected to result in better un- 
derstanding and improved performance of various 
techniques. 

Beyond routine maintenance, minor upgrades 
and a few replacements of faulty components are 
currently planned. They include the replacement 
of SVT modules that are expected to fail in the 
next few years due to radiation damage, plus a 
few others that cannot be correctly read out due 
to broken connections. A large fraction of the 
RPCs are showing gradually increasing losses in 
efficiency and plans are being developed for the 
replacement of the RPC modules over the next 
few years. Furthermore, 20-25 cm of absorber will 
be added to the fiux return to reduce the hadron 
misidentification rates. 

With the expected increase in luminosity, 
machine-induced backgrounds will rise. Measures 
are being prepared to reduce the sources and the 
impact of such backgrounds on BABAR. Apart 
from the addition of shielding against shower de- 
bris, upgrades to the DCH power supply system 
and to the DIRC electronics are presently under 
way. Most important are upgrades to the trigger, 
both at levels LI and L3. Specifically, the DCH 
stereo layer information will be added to allow for 
a more efficient suppression of background tracks 
from outside the luminous region of the beam. 
The L3 processing will be refined so as to reject 
both backgrounds and high rate QED processes 
with higher efficiency. In addition, data acquisi- 
tion and processing capacity will be expanded to 
meet the demands of higher luminosity. 

In summary, the BABAR detector is perform- 
ing very well under current conditions and is well 
suited to record data at significantly higher than 
design luminosity. 
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Abstract 

BABAR, the detector for the SLAC PEP-II asymmetric e+e^ B Factory operating at the T(4S') resonance, 
was designed to allow comprehensive studies of CP-violation in i?-meson decays. Charged particle tracks 
are measured in a multi-layer silicon vertex tracker surrounded by a cylindrical wire drift chamber. 
Electromagentic showers from electrons and photons are detected in an array of Csl crystals located just 
inside the solenoidal coil of a superconducting magnet. Muons and neutral hadrons are identified by arrays 
of resistive plate chambers inserted into gaps in the steel flux return of the magnet. Charged hadrons are 
identified by dE/dx measurements in the tracking detectors and in a ring-imaging Cherenkov detector 
surrounding the drift chamber. The trigger, data acquisition and data-monitoring systems, VME- and 
network-based, are controlled by custom-designed online software. Details of the layout and performance 
of the detector components and their associated electronics and software are presented. 
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